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ABSTRACT 

Several  experiments  are  presented  in  this  thesis  which  examine 
methods  to  measure  and  monitor  fluid  flow  from  hydrothermal  vent  fields. 

Simultaneous  velocity,  temperature,  and  conductivity  data  were 
collected  in  the  convective  flow  emanating  from  a  hydrothermal  vent  field 
located  at  10°56'N,  103°41'W  on  the  East  Pacific  rise.  The  horizontal 
profiles  obtained  indicate  that  the  flow  field  approaches  an  ideal  plume 
in  the  temperature  and  velocity  distribution.  Such  parameters  as  total 
heat  flow  and  maximum  plume  height  can  be  estimated  using  either  the 
velocity  or  the  temperature  information.  The  results  of  these 
independent  calculations  are  in  close  agreement,  yielding  a  total  heat 
flow  from  this  vent  site  of  3.7  +  0.8  MW  and  a  maximum  height  of  150+10 
m.  The  nonlinear  effects  of  large  temperature  variations  on  heat 
capacity  and  volume  changes  slightly  alter  the  calculations  applied  to 
obtain  these  values. 

In  Guaymas  Basin,  a  twelve  day  time  series  of  temperature  data  was 
collected  from  a  point  three  centimeters  above  a  diffuse  hydrothermal 
flow  area.  Using  concurrent  tidal  gauge  data  from  the  town  of  Guaymas  it 
is  shown  that  the  effects  of  tidal  currents  can  be  strong  enough  to 
dominate  the  time  variability  of  a  temperature  signal  at  a  fixed  point  in 
hydrothermal  flow  and  are  a  plausible  explanation  for  the  variations  seen 
in  the  Guaymas  Basin  temperature  data. 

Theoretical  examination  of  hot,  turbulent,  buoyant  jets  exiting  from 
hydrothermal  chimneys  revealed  acoustic  source  mechanisms  capable  of 
producing  sound  at  levels  higher  than  ambient  ocean  noise.  Pressure 
levels  and  frequency  generated  by  hydrothermal  jets  are  dependent  on 
chimney  dimensions,  fluid  velocity  and  temperature  and  therefore  can  be 
used  to  monitor  changes  in  these  parameters  over  time. 

A  laboratory  study  of  low  Mach  number  jet  noise  and  amplification  by 
flow  inhomogeneities  confirmed  theoretical  predictions  for  homogeneous 
jet  noise  power  and  frequency.  The  increase  in  power  due  to  convected 
flow  inhomogeneities,  however,  was  lower  in  the  near  field  than  expected. 

Indirect  evidence  of  hydrothermal  sound  fields  (Reidesel  et  al., 
1982;  Bibee  and  Jacobson,  1986)  showing  anomalous  high  power  and  low 
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frequency  noise  associated  with  vents  is  due  to  processes  other  than  jet 
noise . 

On  Axial  Seamount,  Juan  de  Fuca  Ridge,  high  quality  acoustic  noise 
measurements  were  obtained  by  two  hydrophones  located  3  m  and  40  m  from 
an  active  hydrothermal  vent,  in  an  effort  to  determine  the  feasibility  of 
monitoring  hydrothermal  vent  activity  through  flow  noise  generation. 

Most  of  the  noise  field  could  be  attributed  to  ambient  ocean  noise 
sources  of  microseisms,  distant  shipping  and  weather,  punctuated  by  local 
ships  and  biological  sources.  Water/rock  interface  waves  of  local 
origin,  were  detected  which  showed  high  pressure  amplitudes  near  the 
seafloor  and,  decaying  with  vertical  distance,  produced  low  pressures  at 
40  m  above  the  bottom. 

Detection  of  vent  signals  was  hampered  by  unexpected  spatial  non- 
stationarity  due  to  shadowing  effects  of  the  caldera  wall.  No  continuous 
vent  signals  were  deemed  significant  based  on  a  criterion  of  90% 
probability  of  detection  and  5%  probability  of  false  alarm.  However,  a 
small  signal  near  40  Hz,  with  a  power  level  of  1x10'  Pa  /Hz  was  no^ced 
on  two  records  taken  near  the  Inferno  black  smoker.  The  frequency  of 
this  signal  is  consistent  with  predictions  and  the  power  level  suggests 
the  occurrence  of  jet  noise  amplification  due  to  convected  density 
inhomogeneities . 

Ambient  noise  from  the  TAG  (Trans -Atlantic  Geotraverse)  hydrothermal 
area  on  the  Mid-Atlantic  Ridge  near  26°N,  in  the  frequency  band  1-30  Hz 
at  a  range  of  0. 75-^4  km  from  the  site  of  an  extremely  active  high 
temperature  hydrothermal  vent  field  (Rona ,  1986)  was  examined.  The 
ambient  noise  field  exhibits  great  temporal  and  spatial  variations 
attributed  in  part  to  typical  ocean  noise  sources  such  as  distant 
shipping  and  microseisms.  Power  spectral  levels  as  measured  at  each  of 
six  ocean  bottom  hydrophones  (OBH)  were  used  to  estimate  the  location  of 
point  sources  of  sound  in  the  area,  if  any. 

The  hydrothermal  vent  did  not  produce  enough  sound  to  be  located  as 
a  point  source  using  data  from  the  OBH  array.  The  only  consistently 
identifiable  point  source  found  with  the  data  set  was  generating  sound  in 
a  0.8 -3. 5  Hz  bandwidth  and  located  outside  the  median  valley.  It  appears 
to  be  harmonic  tremor  associated  with  the  tip  of  a  ridge  on  the  western 
side  of  the  spreading  axis  and  may  be  volcanic  in  origin. 
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INTRODUCTION 


The  interaction  of  a  complex  set  of  geological,  biological,  chemical 
and  fluid  dynamical  processes  create  the  life  supporting  environment 
called  earth.  Except  for  the  influx  of  solar  radiation,  the  earth  is 
basically  a  self-contained,  internally  cycling  system  upon  which  a  great 
diversity  of  life  depends.  The  different  cycles  occur  on  all  scales  and 
levels  of  Interaction.  Slow,  large  scale  geologic  processes  such  as 
continental  formation  and  drift,  seafloor  creation,  and  mountain  building 
interact  with  similar  scale  processes  in  other  fields:  species  evolution 
in  biology,  mantle  convection  in  fluid  dynamics,  to  create  long  time 
period  trends  In  the  environment.  Medium  scale  variations  seen  in  ocean 
chemical  composition,  hydrothermal  activity,  sea  level,  biological 
extinction  and  global  climate  are  clearly  cyclic  and  interdependent.  And 
on  the  small  scale,  volcanic  eruptions,  earthquakes,  storms,  pollution, 
and  species  blooms  and  extinctions  interact  to  produce  a  challenging, 
rapidly  changing  habitat.  The  ability  to  predict  and  control  changes  in 
this  environment  is  necessary  for  our  survival  on  such  a  planet, 
especially  as  we  begin  to  participate,  through  technology  and  increased 
population,  in  the  cycles  and  processes  themselves. 

The  ambitious  task  of  understanding  the  earth  as  a  system  must  be 
approached  by  breaking  down  the  problem  into  a  set  of  studies  of 
individual  processes  which  can  then  be  related  to  the  interdisciplinary, 
global  view.  In  this  way,  a  large  number  of  scientists  from  many  fields 
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contribute  to  a  level  of  understanding  necessary  to  prevent  adverse 
changes  and  precipitate  beneficial  ones  in  this  self-contained  system. 

In  the  context  of  earth  system  science,  an  important,  individual 
process  is  the  circulation  of  seawater  through  oceanic  lithosphere,  a 
process  which  transfers,  heat  and  chemicals  from  the  mantle  to  the  oceans 
and  enables  them  to  react  with  and  influence  the  biosphere.  Over  30 
active  hydrothermal  vent  sites  have  been  found  to  date,  both  in  the 
Atlantic  and  Pacific  oceans.  It  is  becoming  increasingly  apparent  that 
hydrothermal  systems  permeate  a  large  portion  of  the  Mid-Ocean  Ridge 
system,  and  may  significantly  contribute  to  the  chemical  composition  of 
the  ocean  and  atmosphere. 

Hydrothermal  circulation  begins  as  cold  seawater  seeps  into  the 
exposed,  fractured  basalt  found  along  mid- ocean  ridges.  The  shallow 
magma  chambers  and  intruded  lavas  characteristic  of  these  ocean  floor 
spreading  centers  provide  heat  to  drive  the  seawater  in  a  convection 
system  within  the  crust.  Unreacted  seawater  enters  a  broad  downflow  zone 
which  extends  kilometers  along  and  away  from  the  unsedimented  spreading 
axis.  As  the  cool  water  descends  and  approaches  the  heat  source,  its 
temperature  and  reactability  increase  and  it  rapidly  leaches  minerals 
from  the  surrounding  rock.  Upon  descending  to  the  top  of  a  magma 
chamber,  2-4  km  into  the  crust,  the  water,  now  rich  in  minerals,  turns 
and  rises  through  narrow  upflow  zones  a  few  tens  of  meters  wide  and  exits 
onto  the  sea  floor.  The  surface  expression  of  these  hydrothermal  vents 
is  a  spectacular  manifestation  of  the  importance  of  the  heat  and  mass 
load  of  the  fluid.  The  precipitation  of  minerals  that  occurs  when  the 
hot  vent  fluid  mixes  with  cold,  ambient  seawater  forms  tall,  15-20  m  high 
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sulfide  chimneys  out  of  which  spews  '->50°C  black,  particle  laden  smoke. 
Fluid  *_r  c  has  mixed  in  the  subsurface  with  ambient  seawater  and  has 
lost  its  mineral  lead  issues  from  cracks  and  fissures  around  these 
chimneys  as  cooler,  clear  fluid.  The  hydrogen  sulfide  present  in  vent 
water  provides  the  primary  energy  source  for  a  suite  of  unique  organisms 
which  colonize  hydrothermal  vents.  Massive  ore  deposits  of  iron,  copper 
and  zinc  sulfide  are  formed  by  continuous  deposition  of  the  mineral 
content  of  vent  fluid.  And  the  plume  of  warm,  particle  rich  fluid  that 
rises  above  a  vent  field  can  be  detected  many  kilometers  away  and  may 
drive  local  ocean  circulation. 

Studies  of  hydrothermal  circulation  have  relied  on  instantaneous 
measurements  of  surface  characteristics  like  morphology,  biota, 
temperature,  chemical  composition,  and  flow  velocity,  which  are  then  can 
be  used  to  model  subsurface  chemical  and  flow  processes.  Even  these 
measurements  are  difficult  to  obtain,  requiring  the  use  of  a  submersible 
for  mosc  sampling.  The  high  pressure,  hot  water,  particle  precipitation 
and  corrosive  fluids  make  long-term  measurements  even  more  scarce.  In 
designing  a  method  to  determine  simply  the  temperature  field,  for 
example,  one  is  faced  with  numerous  difficulties  starting  with  the 
complex  and  variable  morphology  of  vent  fields.  Flow  ranges  in  character 
from  10°C  water  slowly  seeping  through  cracks,  fissures,  masses  of  tube 
worms,  and  bacterial  mats,  to  hot,  350°C  water  jetting  out  of  narrow,  3 
cm  diameter  sulfide  chimneys  at  2-3  m/s.  The  sensor  must  be  able  to 
withstand  high  temperature  and  also  cover  an  area  of  extremely  variable 
flow  and  temperature  Long-term  measurements  are  even  more  problematic, 
as  they  are  complicated  by  the  fact  that  chimneys  grow,  fall  and 
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effectively  change  position,  that  many  materials  degrade  under  vent 
conditions,  and  that  variable  cross-currents  move  the  plume  of  hot  water 
back  and  forth  over  the  vent  field. 

The  full  impact  of  hydrothermal  fluxes  of  heat  and  chemicals  on 
biogeochemical  cycles  can  only  be  assessed  by  determining  both  the  total 
spatial  extent  and  temporal  variability  of  hydrothermal  vents,  neither  of 
which  are  presently  well  known.  The  overall  purpose  of  the  work 
contained  in  this  volume  was  to  examine  heat  and  fluid  fluxes  from 
hydrothermal  vents,  and  to  develop  methods  with  which  to  measure  these, 
both  instantaneously  and  over  the  long-term. 

The  first  study  was  directed  at  measuring  instantaneous  heat  flux 
from  a  black  smoker  vent  field.  Velocity  and  temperature  measurements 
were  taken  in  the  convective  flow  emanating  from  a  hydrothermal  vent 
field  and  used  with  plume  theory  to  estimate  total  heat  flux  from  the 
site.  The  close  agreement  between  heat  flux  calculations  using  velocity 
and  independent  calculations  using  temperature  justified  the  validity  of 
this  method  for  estimating  instantaneous  hydrothermal  vent  heat  flux. 

The  second  study  examined  cross-current  effects  on  temperature 
measurements  made  in  diffuse  hydrothermal  flow.  A  twelve  day  time  series 
of  temperature  from  a  few  centimeters  above  an  area  of  diffuse  flow 
showed  significant  temporal  variations  which  seemed  to  be  correlated  with 
tide.  It  was  shown  through  a  boundary  layer  flow  model  that  the  effects 
of  tidal  currents  can  be  strong  enough  to  dominate  the  time  variability 
of  a  temperature  signal  at  a  fixed  point  in  hydrothermal  flow.  The 
results  show  the  necessity  of  recording  simultaneous  bottom  currents  when 
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trying  to  use  temperature  measurements  to  determine  variations  in 
hydrothermal  flow  caused  by  subsurface  processes. 

The  three  other  studies  which  make  up  this  work  were  concentrated  on 
determining  the  feasibility  of  monitoring  high  temperature,  high  velocity 
fluid  flux  through  passive  acoustic  measurements  of  hydrothermally 
generated  sound.  They  include  a  theoretical  examination  of  sound  source 
mechanisms,  a  laboratory  experiment  to  determine  the  acoustic  behavior  of 
low  Mach  number,  an  examination  of  noise  recorded  near  hydrothermal  vents 
by  microearthquake  surveys,  and  a  field  experiment  designed  to  measure 
the  sound  produced  at  a  hydrothermal  vent  jet. 

Theoretical  examination  of  hot,  turbulent,  buoyant  jets  exiting  from 
hydrothermal  chimneys  revealed  acoustic  source  mechanisms  capable  of 
producing  sound  at  levels  higher  than  ambient  ocean  noise.  Pressure 
levels  and  frequency  generated  by  hydrothermal  jets  are  dependent  on 
chimney  dimensions,  fluid  velocity  and  temperature  and  therefore  can  be 
used  to  monitor  changes  in  these  parameters  over  time. 

A  laboratory  study  of  low  Mach  number  jet  noise  and  amplification  by 
flow  inhomogeneities  confirmed  theoretical  predictions  for  homogeneous 
jet  noise  power  and  frequency.  The  increase  in  power  due  to  convected 
flow  inhomogeneities,  however,  was  lower  in  the  near  field  than  expected. 

Indirect  evidence  of  hydrothermal  sound  fields  showing  anomalous 
high  power  and  low  frequency  noise  associated  with  vents  is  due  to 
processes  other  than  jet  noise. 

Direct  measurements  of  hydrothermal  vent  sound  fields  show  that 
frequency  and  power  levels  are  comparable  to  ambient  ocean  noise  and  are 
consistent  with  jet  quadrupole  near-field  sound  which  has  been  amplified 
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by  the  dipole  behavior  of  convected  flow  inhomogeneities.  This  near¬ 
field  amplification  is  not  as  great  as  that  predicted  for  the  far-field 
but  is  consistent  with  theoretical  considerations  for  near-field  dipole 
and  quadrupole  behavior. 

Hydrothermal  systems  are  extremely  variable,  complex,  inaccessible 
and  difficult  to  characterize.  Yet,  they  embody  a  process  through  which 
significant  quantities  of  heat  and  chemicals  enter  the  ocean,  massive  ore 
deposits  are  formed,  and  on  which  unique  vent  communities  depend.  The 
value  of  studying  them  outweighs  the  difficulties  inherent  in  such 
studies.  The  work  presented  here  has  provided  a  method  with  which  to 
measure  instantaneous  heat  flux  from  black  smoker  vent  fields.  The 
proportion  of  heat  loss  from  diffuse  flow  to  that  of  total  flux  remains 
unknown,  however.  The  importance  of  measuring  cross-currents  in  long¬ 
term  monitoring  has  been  illuminated.  The  method  of  examining  long-term 
variability  in  flow  through  passive  acoustic  monitoring  has  been  shown  to 
be  theoretically  possible,  but  given  the  levels  of  ambient  ocean  noise, 
not  a  simple  tool  useful  at  all  vent  sites. 

The  development  of  hydrothermal  vent  monitoring  techniques  is  still 
an  active  field  where  many  important  problems  remain  unsolved. 
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Simultaneous  velocity,  temperature,  and  conductivity  data  were  taken  in  the  convective  (Tow  emanating 
from  a  hydrothermal  vent  field  located  at  lO^SSTM.  I03“4I'W  on  the  East  Pacilic  Rise  The  horizontal 
profiles  obtained  indicate  that  the  (low  field  approaches  an  ideal  plume  in  the  temperature  and  velocity 
distribution.  Such  parameters  as  total  heat  (Tow  and  maximum  plume  height  can  be  estimated  using  cither 
the  velocity  or  the  temperature  information.  The  results  of  these  independent  calculations  are  in  close 
agreement,  yielding  a  total  heat  flow  from  this  vent  site  of  3  7  ±0  8  MW  and  a  maximum  height  of  150  ±  10 
m.  The  nonlinear  effects  of  large  temperature  variations  on  heat  capacity  and  volume  changes  slightly  alter 


the  calculations  applied  to  obtain  these  values. 

Introduction 

Near  midocean  ridges  (MOR's),  cold  seawater  seeps  into 
unsedimented  crust  and  is  heated  by  the  shallow  magma 
chambers  and  intruded  lavas  of  newly  formed  oceanic  lithos¬ 
phere  in  the  process  called  hydrothermal  circulation.  The  hot 
water,  rich  in  minerals,  rises  buoyantly  through  the  crust  and 
exits  through  hydrothermal  vents  located  on  the  ridge  axis. 
These  relatively  sparsely  distributed  and  inaccessible  vents  are 
geologically,  chemically,  and  biologically  active  areas  currently 
under  intense  study. 

Active  hydrothermal  vents  have  been  found  in  the  Pacific 
Ocean  on  medium  to  fast  spreading  ridges  near  the  Galapagos 
at  86°W,  along  the  East  Pacific  Rise  (EPR)  at  20°S,  ll°N, 
I3°N,  and  2I°N,  in  theGuaymas  Basin,  and  on  the  Gorda  and 
Juan  de  Fuca  ridges  between  40®  and  50°N  [Green  el  at,  1981; 
Baker  el  al.,  1985;  Francheieau  and  Ballard,  1983;  Crane  el  at, 
1985;  Rona  ei  at,  1985].  Numerous  extinct  vent  sulfide  chim¬ 
neys  have  also  been  found  in  these  areas.  In  the  Atlantic  Ocean, 
vents  have  been  photographed  at  23®N  and  26®N  on  the  Mid- 
Atlantic  Ridge  [Rona  el  at,  1984  and  R.  Detrick,  personal 
communication.  1985], 

The  processes  of  hydrothermal  convection  can  be  separated 
into  three  distinct  scales  of  venting.  The  smallest  Kale  is  that  of 
How  from  a  single  opening  on  a  hydrothermal  sulfide  edifice. 
Vents  on  this  scale  are  often  called  smokers  (due  to  the  particu¬ 
lates  in  the  venting  fluid)  and  typically  have  openings  2-8  cm  in 
diameter  [Converse  el  at.  1984).  The  next  Kale  is  that  of  a  vent 
field  or  vent  area.  This  includes  all  the  hot  smokers  (30-4 20°C) 
on  the  sulfide  structure  as  well  as  the  cooler  (J-30°C)  diffuse 
flow  from  cracks  and  crevices  in  the  vicinity  of  the  main  sulfide 
edifice.  A  typical  vent  field  is  100  x  100  m  and  is  usually  defined 
by  the  lateral  extent  of  thermal  anomalies  and  the  associated 
animal  communities  [Francheieau  and  Ballard,  1983].  The  larg¬ 
est  Kale  of  hydrothermal  venting  is  that  of  an  active  ridge 
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segment.  On  this  scale  of  tens  of  kilometers,  hydrothermal 
deposits  and  venting  sites  appear  on  regularly  spaced  (30-80 
km)  topographic  highs  [ Francheieau  and  Ballard.  1983; 
Schouten  ei  at.  .  1 985].  These  highs  are  thought  to  overlie  sites 
of  geologically  recent  (thousands  of  years)  magma  injection 
[Schouten  el  at.  1985).  This  scale  could  include  10-100  km  of 
midocean  ridge  crest  containing  20-30  active  vent  areas. 

Attempts  to  determine  the  thermal  output  of  hydrothermal 
vents  have  been  limited  in  accuracy  by  either  the  very  small  scale 
(e  g.,  a  single  vent  orifice)  or  the  very  large  scale  (e  g.,  a  mido¬ 
cean  ridge  segment)  of  the  measurements.  Macdonald  ei  al. 
[1980]  used  vent  exit  velocities  and  an  orifice  diameter  esti¬ 
mated  from  video  tapes  coupled  with  assumed  exit  tempera¬ 
tures  to  calculate  a  heat  flow  of  60  MW  from  a  single  smoker  at 
the  National  Geographic  vent  site  at  2I°N  on  the  East  Pacific 
Rise.  Converse  ei  al.  [1984]  measured  flow  rates  and  tempera¬ 
tures  at  2I®N  by  inserting  a  small  vane-type  flowmeter  into 
individual  smokers,  estimating  orifice  diameters,  and  using 
buoyant  plume  theory  to  obtain  values  of  heat  flow  ranging 
from  0.5  to  10  MW  for  individual  smokers.  These  heat  flow 
estimates  depend  on  measurements  made  near  the  plume  source 
and  are  thus  highly  sensitive  to  both  the  location  of  velocity  and 
temperature  measurements  and  visual  estimates  of  the  vent  exit 
geometry.  In  addition,  subsequent  estimates  of  total  heat  flow 
from  the  vent  field  obtained  by  multiplying  an  individual  smok¬ 
er’s  heat  flow  by  the  number  of  observed  orifices  can  neither 
account  for  the  variations  of  individual  smoker  size  nor  the 
contribution  to  heat  flow  from  diffuse  (low. 

On  (be  large  scale.  Crane  el  at  [1985]  estimated  heat  loss 
from  ridge  segments  by  applying  buoyant  plume  theory  and  a 
gross  heat  balance  to  temperature  measurements  obtained  by  a 
ship-towed  sensor  array.  Their  estimates  varied  by  a  factor  of 
more  than  100.  between  I3.000and  1,560.000  MW.  for  1 4 1  km 
of  the  Juan  de  Fuca  Ridge,  largely  as  a  result  of  the  sensitivity  of 
heat  flow  calculations  to  the  resolution  of  horizontal  and  verti¬ 
cal  positioning  of  temperature  measurements  and  to  uncertain¬ 
ties  in  the  magnitude  of  regional  ocean  currents  and 
hydrography. 


I  ll  II  F  FI  Al  HMIXOIMI  KMAI  VfNI  1*1  l<ME  Kl  ow 


\ 


\ 


<  > 


oppronmoftf  boyndory  of  v#nf  fiold 
boi«d  on  (r«qu«ncf  of  biological  activity 

logit,  dip  di'«<*iOft 
tifUlft 

Iqvo  lok«,  direction  of  collopi* 
boioll  »K««l  flow 


Legend 


pillow  lovO 
(  bulboui ) 

pillow  lovo 


m 

\ 


Pillow  lava  With  ftligb*  dusting 
of  t«dim«nt  on  ond  b«iw««n 

pillow  lo.o  w 1 1 K  il.ghl  lob'* 
dulling  ol  todimont  on  ond 
bftwoon  piliowi 

ini«ni«  conctnlfolion  of 
»#»pulid  wotms 
inl»m#  conconlrotion  of 

tub*  worm*  Ond  mulioll 

mottiv*  lull'd* 


lolui  blocks  common 
ol  boia  of  scorp 


\\ 

\\ 
l  \ 

ilNVzQNE 

r.dga  of  .  j  T\  Z 

durupiad  i  \  \ 

bo  Foil  I  I 

A  \\  'S§ifc:->3» 

\\ZONEflp.y:  ^rr* 

voiloy  Moor  cooiom*  I  !\  * 
ir-oo»h,tow  profit*  j  \  U  1  {. 

lobotr  flows  ond  ---  I  )  \  f  l  || 

sh««t  flows;  !  j  <J"  ►  11 

no  sadimanl  covar  -J- »  I  y 

, 

1  2*3*m-K»gh  rim 

Q  5  10  20  3p  4Q 

m 


Cron  Section  A- A  looking  north 
(  no  vorlicol  •  ■oggarolion  ) 


Fig.  I.  Map  of  vent  field  il  !0*55.73'N.  103*40  6'Wihowingsulfide  Jlructureonedge  of  fault  scnrp(from 
MvConachv  n  al.  1986].  Zone  I  depicts  the  location  of  warm  water  seeps;  lone  2  shows  the  location  of  the 
main  sulfide  edifice  with  black  smokers. 


Our  objectives  were  to  obtain  more  accurate  heat  flow  esti¬ 
mates  and  to  characterize  the  fluid  How  from  a  vent  field  by 
measuring  the  three-dimensional  structure  of  the  convective 
flow  above  a  vent  field  at  a  scale  (meters)  between  that  of 
MacDonald  and  Converse  et  al.  (centimeters)  and  Crane  et  al. 
(kilometers).  Measurements  made  al  the  scale  of  a  vent  field  can 
avoid  problems  due  to  sensitivity  to  vent  exit  quantity  and 
geometry,  and  depend  minimally  on  estimates  of  oceanic  condi¬ 
tions  such  as  ambient  current  and  stratification.  Precise  map¬ 
ping  of  temperature  anomalies  also  permits  the  quantification 
of  the  relative  contributions  of  high-temperature  smokers  and 
lower  temperature  diffuse  flow  to  the  total  heat  flow.  To  collect 
the  appropriate  data,  the  DSV  Alvin  and  R/V  Atlantis  II 
visited  I  l°N  on  the  EPR  in  May  1984  to  survey  chemically  and 
physically  the  water  above  a  vent  field  composed  of  both  black 
smokers  and  the  more  diffuse  warm  vents,  temperature,  pres¬ 


sure,  conductivity,  and  velocity  were  measured  at  various 
heights  above  this  vent  field  and  used  to  determine  the  physical 
structure  of  the  flow  and  to  calculate  maximum  plume  height 
and  total  advective  heat  flow. 

Sue  Dfscripiion 

The  target  vent  field,  located  at  IO°55.73'N.  I03°40.6'W  at  a 
depth  of  2531  m  (see  Figure  I),  contains  an  actively  venting 
sulfide  edifice  surrounded  by  patches  of  warm  venting.  The  65  x 
45  m  vent  field  lies  along  the  eastern  edge  of  the  rift  valley 
against  a  20  m-high  north-south  trending  fault  scarp  [McCon- 
achv  et  al..  1984,  1986].  The  vent  area  is  comprised  of  pillow 
lava,  basalt  sheet  flows,  lava  lakes  with  lava  spires,  and  talus 
blocks  of  sulfide  near  the  smokers.  The  chimney  structure, 
consisting  of  both  dead  and  broken  spires  and  several  active 
regions,  reaches  a  height  of  approximately  8  m  and  is  3-4  m  wide 
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Fig  2.  Diagram  of  50  m  cable  array  deployed  from  the  DSR/V  Alvin 
with  a  water  sampler,  transmissamelcr.  current  meter  and  nine 
thermistors. 

at  the  base"  Hot.  gray-black  water  flowed  from  the  base  from 
three  main  and  three  or  four  lesser  openings  2-8  cm  in  diameter, 
all  within  I  m  of  the  seafloor.  Shimmering,  milky  white  water 
seeped  out  several  meters  from  the  base  through  beds  of  vesli- 
mentiferan  tube  worms. 

Equipment 

The  temperature,  conductivity,  and  velocity  fields  above  the 
hydrothermal  vent  field  were  sampled  with  an  array  of  sensors 
deployed  vertically  from  the  DR  V  Alvin  (Figure  2).  The  instru¬ 
ments  used  in  this  work  included  a  vertically  mounted  seven- 
thermistor  chain  surrounding  a  self-contained  flowmeter, 
which  recorded  three-component  current  velocity,  tempera¬ 
ture.  conductivity,  pressure,  tilt,  and  compass  heading.  The 
flowmeter  was  positioned  24  m  above  the  sampling  basket 
attached  to  base  of  the  submarine;  the  thermistors  were 
mounted  at  40.5,  38.2,  35.6,  34.2,  30.1,  25.4,  and  22.5  m  above 
the  base. 

The  flowmeter,  built  by  J.  Dean  at  the  Woods  Hole  Oceano¬ 
graphic  Institution,  operates  in  the  range  of  i-100  cm/s(  Weller 
el  at.,  1985].  Its  four  fans,  each  22  cm  in  diameter,  detect  two 
vertical  (for  redundancy)  and  two  horizontal  components  of 
velocity.  A  tilt  meter  and  a  compass  enable  the  measured  veloci¬ 
ties  to  be  transposed  into  an  earth-based  reference  frame  despite 
rotations  of  the  flowmeter.  The  thermistor  response  time  is  7  s 
with  a  resolution  of  0  005°C.  Tilt,  temperature,  and  pressure 
are  each  averaged  over  5  s,  conductivity  over  3.6  s,  and  the 
current  velocities  over  the  entire  14  s  sampling  interval. 

The  thermistor  chain,  consisting  of  seven  thermistors  and  a 
self-contained,  internally  recording  package,  was  deployed  on 
dive  1386.  A  data  sample  with  a  resolution  of  Q.OOI°C  was 
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I  lie  Ml  m  cable,  with  ils  attached  I  low  meter,  thermistor  chain, 
ami  I lol :il ion.  was  connected  hv  a  quick  tclcnsc  mechanism  to 
the  basket  ol  Alvui  prior  to  launch  During  descent,  tcmpcia- 
ture.  conductivity,  and  pressure  were  recorded  to  obtain  the 
ambient  stratification  of  the  water  column  Following  descent, 
each  dive  consisted  of  several  stations  in  the  vicinity  of  the 
active  main  sullide  edifice  At  each  station  the  submarine 
remained  stationary  on  the  bottom  for  5-30  min  to  allow  the 
instrument  array  to  stabilize  in  a  vertical  position  and  to  record 
the  time  variability  of  the  turbulent  fluctuations  in  the  plume 
The  submarine  moved  between  stations  while  remaining  as 
close  to  the  bottom  as  possible,  usually  within  1-2  m  I  his 
resulted  in  flowmeter  positioning  of  24  ±  2  m  above  the  seafloor 
during  most  of  the  dives  I  arger  excursions  can  be  seen  by 
examining  the  pressure  records  All  stations  were  conducted  in 
the  vent  field  at  distances  of  1-35  m  from  (he  chimney  Naviga¬ 
tion  was  by  dead  reckoning  Video  and  still  pi'-"*"'’  were  taken 
of  the  chimney  structure  and  environs 

I  hree  dives  were  completed  with  a  fully  operational,  continu¬ 
ously  recording  flowmeter  ( 1 384.  1385,  1386).  and  onlhclasiof 
these  the  thermistor  chain  internally  recorded  temperature  at 
seven  points  along  the  vertical  array. 

Data  Reduciion 

The  experiment  resulted  in  a  data  set  that  includes  seven 
simultaneous  temperature  profiles  of  the  plume  from  one  dive 
and  three  sets  of  temperature,  conductivity,  pressure,  three- 
component  velocity,  tilt,  and  compass  taken  at  24  m  above  the 
bottom  from  each  of  three  dives 

Salinity  was  calculated  via  a  salinity  algorithm [Fofimntfand 
Millard.  1983]  using  the  temperature,  conductivity,  and  pres¬ 
sure  records.  7  he  few-second  time  lag  between  measurements  of 
temperature  and  conductivity  was  sufficient  to  cause  errors  in 
salinity  calculations  in  the  rapidly  varying  plume  environment 
However,  the  chloride  content  of  a  water  sample  from  this  vent 
taken  on  dive  1377.  several  weeks  earlier,  was  found  to  be  only 
3Fr  higher  than  that  of  the  ambient  water  (.1  Edmond,  personal 
communication.  1985)  Accordingly,  it  was  assumed  in  further 
data  analysis  that  the  plume  salinity  was  equal  to  the  ambient 
value  of  34.651°/,*,  obtained  during  the  dive  series  1384-1386 
Plume  density  was  calculated  employing  an  equation  of  state 
for  seawater.  [Fofonoff  and  Millard.  1983]  and  assuming  a 
constant  ambient  salinity  and  pressure  (2535  x  I0J  Pa). 

Corrections  are  applied  in  order  to  recover  the  true  vertical 
plume  velocities  in  cases  w  hen  the  instrument  array  was  moving 
up  or  down  in  response  to  the  submarine  motion.  I  hese  vertical 
motions  are  removed  from  the  measured  vertical  velocities  to 
obtain  true  vertical  plume  velocities  by  subtracting  the  rate  of 
pressure  change  of  the  instrument  from  the  measured  velocities. 
This  procedure  introduces  an  unwanted  but  tractable  amount 
of  noise  (±4  cm/s)  into  the  resulting  vertical  plume  velocities 
due  to  (he  limited  resolution  of  the  pressure  sensor. 

Both  the  vertical  and  horizontal  velocities  have  been  pm  in  an 
earth-based  reference  frame  by  using  the  compass  and  tilt 
records  to  remove  the  effects  of  spinning  and  tipping  1  hese 
horizontal  velocity  records  were  used  to  calculate  approximate 
instrument  trajectories  because  an  operational  bottom  naviga¬ 
tion  system  was  unavailable  during  the  dives  The  trajectories 
were  reconstructed  by  estimating  the  ambient  current  from  the 
long-term  stations  and  subtracting  this  from  the  horizontal 
vclocilv  data  I  he  resultant  velocities  were  integrated  to  obtain 


recorded  from  each  thermistor  every  24  s. 
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the  relative  positions  of  the  instrument  during  a  given  dive  By 
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Fig  3.  Approximate  dive  Crack  from  dive  1384  shown  with  annotated 
temperatures  and  ambient  current  vector  (0  24  cm/s  W,  3  4  cm/s  N). 
Track  obtained  from  horizontal  current  records  (accurate  to  10m)  The 
approximate  location  of  the  main  sulfide  chimney  is  outlined  with  a 
circle. 

combining  temperature  records  with  these  positions  we  were 
able  to  make  an  approximate,  two-dimensional  map  of  the 
measured  temperature  anomalies  24  m  off  the  bottom  (see 
Figure  3).  The  major  source  of  error  in  position  (on  the  order  of 


10  ml  results  from  the  assumption  of  constant  ambient  current, 
which  ranged  from  3  to  4  cm/s  toward  the  north  to  northwest 
during  Ihe  dive 


Ri  sm  t s  and  Dim  ussion 

In  order  to  characterize  the  fluid  flow  from  a  hydrothermal 
vent  field  and  improve  the  accuracy  of  quantifying  advected 
heat  flow  it  is  important  to  determine  the  extent  to  which  the 
advective  flow  is  plumelikc  and  how  it  may  be  complicated  by 
multiple  sources,  diffuse  flow,  ambient  stratification,  or  shear 
flow  caused  by  cross  currents. 

To  determine  Ihe  form  of  the  conceded  flow,  it  is  noted  that 
when  the  instruments  intersect  the  vent  field  at  a  height  of  24  m 
above  Ihe  bottom,  the  data  show  constant  ambient  values  of 
I  83°C  punctuated  by  sharp  peaks  in  temperature,  velocity,  and 
conductivity  (see  Figures  4-7).  Ihe  data  obtained  from  the 
thermistors  at  seven  different  heights  also  show  this  pattern 
with  similarly  shaped  anomalous  peaks  in  temperature.  This 
indicates  a  plumelike  structure,  characterized  by  positive, 
peaked  anomalies,  where  the  overall  structure  of  each  anomaly 
is  present  at  more  than  one  height  (Figure  7)  with  a  general 
decrease  in  amplitude  as  height  increases 

The  four  complicating  factors  mentioned  above  are  all  pres¬ 
ent  at  this  vent  site  and  must  initially  be  addressed  todetermine 
ihe  applicability  of  simple  plume  theory  at  this  site,  and  later 
will  be  tested  more  rigorously  with  the  specific  values  calculated 
from  the  theory.  (I)  Multiple  sources:  Photos  show  three 
strongly  venting  orifices  and  three  or  four  weaker  ones  located 
at  various  levels  on  the  main  sulfide  deposit  within  3  m  horizon¬ 
tally  of  each  other.  (2)  Diffuse  flow:  Several  pockets  of  cooler 
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Fig.  4.  Dive  1384  vertical  velocity,  temperature,  conducts  ity.  pressure. and  density  profiles  through  the 
vent  Held  The  ascent  at  the  end  of  ihe  di\ccan  be  seen  at  21  5  hours  The  spikes  in  the  velocity  result  from 
the  corrections  applied  to  compensate  for  the  vertical  motions  of  the  submarine  (±4  cm/s). 
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Fig.  5.  Dive  1 385  vertical  velocity,  temperature,  conductivity.  pressure.  and  density  data  profile*:  through 
the  vent  field  The  descent  at  the  beginning  of  the  dive  can  be  seen  at  17  hours.  ^  he  pressure  drop  indicates 
the  occurrence  of  a  submarine  reconnaissance  excursion  off  the  seafloor. 


diffuse  flow  are  visible  wilhin  20  m  of  ihe  base  of  the  sulfide 
edifice.  (3)  Ambient  stratification:  Ambient  stratification  moni¬ 
tored  during  ascent  and  descent  indicated  that  the  ambient 
density  of  I039.3330kg/mjat  2535  m  remained  constant  uptoa 
height  of  75  m  off  the  bottom,  after  which  it  dropped  steadily  (in 
50  m  averages)  to  1039.3185  kg/m1  at  2235  m.  (4)  Shear  flow: 
The  ambient  current  remained  fairly  constant  during  the  3  days 
at  3-4  cm/s  toward  the  north  to  northwest. 

If  each  of  the  six  or  seven  orifices  was  emitting  a  single  plume, 
the  diameters  of  each  would  be  approximately  5  m  at  a  height  of 
24  m  above  the  base  of  the  sulfide  edifice.  The  plume  center- 
lines.  however,  are  no  more  than  3  m  apart,  so  by  24  m.  the 
plumes  will  have  merged  to  some  extent.  1  he  temperature 
sampling  interval  of  14  s  coupled  with  a  bottom  speed  ranging 
from  5  to  50  cm/s  resulted  in  a  typical  sample  spacing  of  from 
0.7  to  7.0  m.  This  makes  it  difficult  to  determine  whether  unique 
signatures  of  the  plumes  from  individual  orifices  exist  at  24  m. 
No  evidence  of  a  distributed,  diffuse,  convective  flow  is 
observed  in  the  recorded  data.  The  effects  of  the  diffuse  flow  are 
either  below  the  limits  of  detection  or  have  been  incorporated 
through  entrainment  into  the  main  plume.  On  the  basis  of  these 
measurements,  we  hypothesise  that  at  this  vent  site  the  convec¬ 
tive  heat  transport  is  primarily  in  the  form  of  a  single  plume 
created  by  the  merging  of  multiple  individual  smokers  Further¬ 
more.  although  the  plume  was  measurably  deflected  by  a  lateral 
current  (no  anomalies  were  found  when  the  submarine  was  on 
the  south  side  of  the  smokers),  it  is  hypothesised  that  this 
current  has  no  significant  effect  on  calculations  yielding  heat 
flow  and  maximum  plume  height  Finally,  there  is  very  little 
ambient  density  stratification  over  the  height  range  of  measure¬ 
ments  (50  m).  This  characterization  of  the  vent  field  (low  pro¬ 


vides  the  basis  for  heat  flow  estimates  developed  in  the  next 
section. 

Simple  Plume  Approximation 

In  a  pure  ideal  plume,  the  gravitational  force  acts  on  density 
differences  between  the  source  fluid  and  its  environment  and. 
for  a  buoyant  source  fluid,  an  upwardly  expanding  plume 
results.  In  a  pure  jet,  fluid  enters  an  ambient  fluid  of  equal 
density,  and  the  flow  is  governed  entirely  by  the  initial  momen¬ 
tum.  Many  flows,  called  buoyant  jets,  fall  somewhere  between 
the  two,  with  both  initial  momentum  and  either  positive  or 
negative  buoyancy  relative  to  the  ambient  fluid.  Although  jets 
and  plumes  behave  differently  near  the  source,  all  buoyant  jets 
will  emulate  ideal  plumes  at  distances  far  from  the  source 
[Turner,  1973).  In  the  following  formulation  we  use  simple 
plume  theory  to  examine  quantitatively  the  hypothesis  that  our 
vent  field  plume  can  be  considered  an  ideal  plume  exiting  into  a 
uniform,  quiescent  environment. 

The  behavior  of  these  plumes  can  be  approximated  for  small 
temperature  differences  by  simple  plume  theory  [Fischer  el  at., 
1979]  in  which  the  plume  is  characterized  by  a  given  initial 
buoyancy  flux  from  a  point  source: 

So  =  (og/fl/(C»  (I) 

where  II  is  the  heat  flux  added  by  the  source,  a  is  the  density  of 
the  fluid,  a  is  the  coefficient  of  thermal  expansion  ( I  48  x  10'* 
°C"’  at  ambient  temperature),  g  is  gravitational  acceleration, 
and  Cp  is  the  specific  heat  (4200  J  kg  '°C''  at  ambient  tempera¬ 
ture)  In  a  simple  plume,  with  no  initial  momentum  flux,  dimen¬ 
sional  analysis  results  in  equations  relating  the  centerline  axial 
velocity.  IF.  to  initial  buoyancy  flux  So  and  height  Z  above  a 
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Fig  6.  Dive  1 386  vertical  velocity,  temperature.  conductivity,  pressure,  and  density  profiles  through  the 
vent  field.  The  descent  at  the  beginning  of  the  dive  can  be  seen  at  I6.S  hours. 


point  source: 

W7—  Ci(Bo/Z)'r3  (2) 

and  centerline  temperature  T,  to  initial  buoyancy  flux  and 
height: 

r=  Gtffl.’Z'V'Vg)  (31 

The  dimensionless  constants  Ci  and  Cj  have  been  obtained 
empirically:  4.7  for  G  by  Rouse  el  at  ,  [1952]  and  9.1  for  Ca 
[ Chen  and  Rodi.  1980].  Real  plumes  do  not  exit  from  point 
sources,  but  Morion  el  at.  [1956]  have  demonstrated  that  all 
plumes  exiting  from  finite  sources  can  be  modeled  as  flows 
generated  from  virtual  sources  a  distance  Z„  below  the  real 
source.  To  obtain  other  plume  parameters,  the  equations  relat¬ 
ing  initial  volume  flux.  Q„,  and  specific  momentum  (lux  M,  to 
buoyancy  flux  and  height  are 

<?.  =  CsB^'^Zo13  (4) 

Mo  =  Gflo^',Z.4,  (5) 

In  addition,  the  plume  radius  can  be  related  to  Q  and  W: 

/?  =  [(>/<*  HO  J’'1  (6) 

where  Cj  and  Ca.  obtained  using  an  assumption  of  Gaussian 
profiles  for  temperature  and  velocity  [Fischer  el  at..  1979],  are 
0. 1 5  and  0.3 1 5,  respectively.  (The  value  of  Ca  is  slight  ly  different 
from  that  given  by  Fischer  et  al.  because  we  have  taken  into 
account  the  contribution  of  both  mean  and  turbulent  fluxes.) 
The  initial  values  of  plume  parameters  such  as  buoyancy  flux 


(Bo),  volume  discharge  (Qo).  specific  momentum  discharge 
(Mo),  radius  (Ro).  and  velocity!  H'oJcan  be  calculated  given  an 
initial  temperature  (To),  and  either  the  centerline  temperature 
or  velocity  at  a  known  height.  An  average  centerline  tempera¬ 
ture  anomaly  of  0.27  ±  0  09’ C  at  24  m  above  the  bottom  and  an 
average  centerline  velocity  of  0. 15  ±  0.03  m/ sat  24  m  have  been 


Fig  7  Simultaneous  temperature  profiles  as  a  function  of  height  from 
dive  1386  Profiles  were  taken  al  22  5.  25  4,  30  1 .  34  2.  35  6.  38  2.  and 
40  5  m  above  the  Alvin  sampling  basket  Bars  on  right  tndicat*  aO  20°  C 
temperature  differential  The  abruptness  of  the  temperature  peaks  is 
due  m  the  sharpness  of  the  hot/cold  boundary  characteristic  of  plumes, 
coupled  with  the  length  of  sampling  interval  compared  with  the  width  of 
the  plume. 
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estimates  is  negligible  because  at  24  in  the  (calculated)  initial 
vertical  velocity  of  63  cm  /s  has  been  reduced  to  a  (measured) 
value  15  cm/s,  which  is  still  a  factor  rd  4  times  ihc magiiiludc  of 
the  ambient  cross  current  I’lumc  dellcctions  at  this  height  are 
estimated  to  be  of  the  order  of  10  m,  which  is  consistent  with  the 
observed  horirontal  distribution  of  anomalies. 

The  maximum  height  of  the  plume  in  a  stratified  environ¬ 
ment  can  be  calculated  from  the  heal  flow  and  density  stratifica¬ 
tion.  Specifically,  from  Campbell  el  al,  [1985], 

Zm.,  =  S^BJn)",N-3^  (8) 


Fig  8.  (Top)  Volume  versus  temperature  for  NaCI  solution  lower 
points  (0-40”  C)  from  Fofonoff  and  V/ illard  [  1 984 ].  upper  points  from 
C7ien(l98l]  Solid  curve  depicts  function  used  for  modeling  (Bottom) 
Enthalpy  versus  temperature  for  HiO  solution  from  Keentn  [1969]. 

-j-ng  the  maximum  values  from  the  data  shown  in 
Figures  4-6.  A  vent  exit  temperature  of  338°C  was  measured 
with  the  submarine's  high  temperature  probe  at  one  of  the 
orifices.  Using  the  above  equations  and  the  average  centerline 
temperature,  the  following  values  for  the  relevant  plume 
parameters  are  obtained: 

Zo  =  0  33  m 
=  0.63  m/s 
«o  =  0  033  m 
B,  =  7.9  x  I0‘4  m4/s3 
Q„  =  2.2  x  I0'3  m3/s 
M.  =  6.1  x  10  4  m4/sJ 

These  values  are  consistent  (within  10%)  with  those  calcu¬ 
lated  from  the  average  centerline  velocity.  These  initial  values 
can  now  be  used  to  check  the  validity  of  certain  assumptions 
made  when  the  simple  equations  presented  above  were  applied 
to  this  vent.  These  assumptions  relate  to  the  effects  of  source 
geometry  initial  momentum,  viscosity,  ambient  stratification, 
and  shear  flow  on  plume  behavior. 

The  simple  plume  relationships  presented  above  may  be  app¬ 
lied  to  observed  plumes  provided  that  the  comparisons  with 
measurements  are  made  beyond  the  distance  lo  =  Qol  Mo'1*. 
within  which  the  source  geometry  is  important,  and 
U  =  within  which  initial  momentum  influences 

the  plume  [Fischer  ei  al  .  1979).  Using  the  above  values  for  the 
hydrothermal  vents  studied  here,  it  is  found  that  l<j  —  0  09  m 
and  lm  =  0.14  m.  The  data  in  this  experiment  were  collected 
well  beyond  these  distances  The  effects  of  viscosity  can  be 
neglected  provided  that  the  Reynolds  number 

Re  =  {irR*W*)',7/v  (7) 

(where  v  is  the  kinematic  viscosity)  is  greater  than  about  4000 
[Labus  and  Symons,  1972).  The  Reynolds  number  at  the  exit  ot 
the  hydrothermal  vent  is  about  4  x  10*  (v  =  lx  1 0‘*  m2/s). 

The  effect  of  the  3-4  cm/s  lateral  current  on  the  heat  flow 


where  Zm„  is  the  maximum  height  off  the  bottom  and  N  is  the 
buoyancy  frequency 

N  =  [{-gdp)Hpdz)]'  2  (9) 

where  p  is  a  reference  density  of  seawater 

Using  data  obtained  on  dive  descent  and  ascent,  an  average 
density  gradient  of  4  8  x  I0‘s  kg/ m4  for  the  bottom  300  m  was 
recorded  (averaging  in  50  m  intervals).  With  a  buoyancy  flux  of 
7  9  x  10  4  m4/s3and  /V  =  6.8  x  I0  4H7,  the  plume  should  reach 
maximum  height  of  150  m.  at  which  point  the  plume  density 
reaches  ambient  densiiv  and  it  begins  to  spread  laterally.  Since 
the  plume  does  not  begin  to  spread  until  150  m.  simple  plume 
theory  (which  ignores  density  stratification  effects)  can  be  used 
at  heights  below  this  and  is  valid  for  the  range  of  measurements 
in  this  experiment  (up  to  50  m). 

Nonlinear  Plume  Model 

The  discussion  above  justifies  the  neglect  of  plume  momen¬ 
tum,  ambient  stratification,  and  horirontal  currents  in  analvr- 
ing  the  plume  data  A  potential  problem  with  using  the  simple 
plume  model  is  that  the  environment  of  hydrothermal  plumes 
comprises  seawater,  high  pressure  (2500  x  l04Pa),  and  large 
temperature  gradients  ( 2-350°  C)  The  specific  volume  (  4")  and 
the  enthalpy  (F)  are  nonlinear  functions  of  temperature  (see 
Figure  8)  and,  in  the  temperature  and  pressure  ranges  at  a 
hydrothermal  vent,  linearapproximalionsto  Fand  Fcannot  be 
used  ( i.c. .  the  coefficient  of  thermal  expansion  a  and  the  specific 
heat  Cp  of  seawater  are  not  constant)  As  a  consequence,  the 
plume  buoyancy  flux  is  not  conserved  but  varies  with  height 
Koisovmos  [  1975]  has  shown  that  corrections  should  be  made 
in  plume  behavior  models  for  even  moderate  differences  in 
temperature.  Instead  of  a  posteriori  use  of  nonlinear 
temperature-density  or  temperature-enthalpy  relationships  (as 
was  done  by  Converse  el  al  [1984])  to  correct  a  model  based  on 
buoyancy  flux  conservation,  a  plume  model  based  on  enthalpy 
conservation  that  takes  into  account  the  nonlinear  dependence 
of  enthalpy  and  volume  with  temperature  must  be  used  This 
more  general  approach  to  the  buoyant  plume  problem  uses  the 
Morion  el  al.  [1956]  entrainment  hypothesis  coupled  with  an 
assumption  of  similarity  of  horirontal  profiles  of  velocity  and 
temperature.  This  method,  which  uses  a  differential  formula¬ 
tion.  is  described  in  the  appendix.  This  approach  also  allows 
plume  and  environmental  parameters  such  as  ambient  stratifi¬ 
cation  or  cross  llow  to  enter  the  problem  as  initial  conditions, 
boundary  conditions  and  coefficients  in  a  system  of  differential 
equations,  although  such  calculations  are  not  presented  in  this 
paper. 

1  he  numerical  solution  o(  Ihc  differential  model  yields  the 
same  results  as  the  simple  plume  model  when  the  coefficients  of 
thermal  expansion  and  specific  heat  are  held  constant  When 
they  are  allowed  to  vary  as  expected,  a  comparison  of  simple 
plume  theory  with  nonlinear  formulation  indicates  that  the 
nonlinear  variation  of  density  and  enthalpv  with  temperature 
has  some  eflecl  on  the  plume  temperature  and  velocity  profiles 


Lithe  et  al  Hvuroihermai.  Vent  Pi  ume  F'iow 


Temperature  *C 

Fig.  9.  Height  versus  temperature  for  initial  heat  flow  values  of  (from 
left  to  right)  I.  2,  3.  and  4  MW  obtained  using  Ihe  nonlinear  model. 
Peak  temperature  values  from  dives  1384  (point  A).  1385  (point  B)and 
1386  (point  C)  taken  at  24  m  above  the  submarine. 


Temperature  °C 

Fig  II  Height  versus  velocity  for  inilial  heat  flow  values  of  (from  led 
to  right)  1 .  2.  3.  and  4  MW  obtained  using  Ihe  nonlinear  model  Peak 
temperatures  from  dive  1386  taken  at  40  5.  38  2.  35.6.  34  2.  30  I.  25.4 
and  22.5  m  above  the  submarine. 


when  applied  to  a  vent  field.  For  a  vertical  velocity  of  14.4  cm/s 
at  24  m.  the  sir""1'  plume  model  predicts  a  heat  flow  of  2.02 
MW.  while  the  nonlinear  model  predicts  1 .97  M  W.  a  difference 
of  3%.  For  a  temperature  anomaly  at  24  m  of  0.34I°C,  the 
simple  plume  model  yields  3.40  M  W,  whereas  the  nonlinear  one 
produces  2.76  MW,  a  difference  of  20%.  While  these  differences 
are  modest  compared  with  possible  errors  in  the  measurements, 
the  nonlinear  model  provides,  at  reasonable  computational 
cost,  a  theoretically  consistent  prediction  of  plume  temperature 
and  velocity.  Accordingly,  the  nonlinear  model  will  be  used  in 
the  next  section  to  estimate  heat  flows. 

Heat  Flow  Estimates 

Heat  flow  values  using  the  nonlinear  model  can  be  obtained 
by  Fitting  either  T  versus  Z  or  W  versus  Z  to  ihe  peak  tempera¬ 
ture  and  velocity  data.  Heat  flow  estimates  resulting  from  using 
T and  W'are  independent  of  each  other,  and  their  comparison 
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Fig.  10  Height  versus  velocity  for  initial  heat  flow  values  of  (from  left 
to  right)  I.  2,  3,  and  4  MW  obtained  using  the  nonlinear  model  Peak 
velocity  values  from  dives  1384  (point  A).  1385  (point  B).  and  1386 
(point  C)  taken  al  24  m  above  ihe  submarine. 


can  be  used  as  a  measure  of  confidence  in  the  final  heat  flow 
estimate  Maximum  peak  values  of  temperature  anomaly  and 
vertical  velocity  from  each  of  three  dives  were  used  in  the  model: 
it  is  assumed  that  the  lesser  peaks  on  Ihe  F(or  HO  versus  time 
records  were  due  to  passing  through  the  plume  hut  not  through 
its  center.  To  obtain  the  heat  flow,  the  maximum  (emperatures 
(Figure  9)  and  velocities  (Figure  10)  measured  wilh  the  flowme¬ 
ter  from  each  of  three  dives  are  plotted  with  the  calculated 
temperature  versus  height  and  velocity  versus  height  curves  for 
various  inilial  heat  flows  generated  by  the  nonlinear  enthalpy 
and  thermal  expansion  model  I  he  seven-thermistorchain  used 
on  dive  1386  yields  additional  lempcralure  maxima  at  different 
heights  ( Figure  1 1 ).  A  summary  of  the  data  and  resultant  heat 
flow  is  given  in  Table  I. 

The  Final  estimate  of  the  total  heat  flow  must  take  into 
account  the  need  of  the  plume  models  for  an  input  of  centerline 
temperature  or  velocity.  Since  rcal-lime  temperature  was  not 
available  during  the  dive,  we  cannot  he  sure  that  the  peak  values 
were  cenlerline  values.  It  is  likely  that  some  of  the  peak  temper¬ 
ature  values  were  recorded  off-centerline  and  thus  are  lower 
than  the  true  cenlerline  temperatures.  T  his  would  imply  that  the 
maximum  temperature  recorded  during  the  cxneriment  should 
be  chosen  and,  consequently,  the  maximum  heat  flow  obtained 
wilh  the  nonlinear  plume  model.  However,  the  models  also 
assume  that  the  centerline  value  is  a  mean  value,  averaged  over 
a  long  time  to  remove  the  effects  of  turbulent  fluctuations. 
Unfortunately,  again  because  electronics  problems  in  Alvin 
eliminated  real-time  readouts,  we  were  not  able  to  position  the 
sensor  in  one  spot  of  the  plume  for  the  necessary  amount  of 
time.  Therefore  the  temporal  temperature  variability  of  the 
plume  due  to  turbulent  fluctuations  is  not  known.  If  the  abso¬ 
lute  peak  temperature  measured  during  the  experiment  is 
chosen  as  a  representative  mean  centerline  temperature,  too 
high  a  heat  flow  estimate  may  result.  However,  because  the 
temperature  and  velocity  values  are  so  sensitive  to  lateral  posi¬ 
tion  in  ihe  plume,  it  is  more  likely  that  the  heal  flow  values  are 
biased  on  the  low  side.  The  average  value  of  all  Ihe  heat  flow 
estimates  is  2.8  MW  and  the  maximum  value  is  4  6  MW  Given 
the  uncertainties  mentioned  above,  a  reasonable  heat  flow  esti¬ 
mate  for  the  vent  field  is  3.7  ±  0  8  MW. 

.  Heal  flow  can  also  be  calculated  directly  by  mulliplying  point 
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TAHLE  I  Nonlinear  Model  Ileal  How  (Maximum  Value  Hrom  Each  Dive) 


Dive 

Height*  m 

1  emperature 
Anomaly  °C 

Heat 

1  lev  MW 

Velocity  cm/* 

Ileal 

J  low  MW 

I.1H4 

24  1) 

f) 

2  2* 

IN  K 

4  21 

1 3K  5 

24  (1 

O  f*X 

o  ‘>r 

12  1 

1  2| 

1386 

24  0 

0  Wl 

2  Ht 

14  4 

2.01 

40  5 

0  201 

4  4 

38  2 

0  230 

4  6 

35  6 

0  23) 

4  1 

34  2 

0  210 

3  1 

30  1 

0  179 

19 

25  4 

0  2)2 

1  8 

22  5 

0  453 

3  7 

•  Above  submarine 

t  Value  derived  from  data  obtained  by  flowmeter  instrumentation. 


velocity  and  temperature  measurements  This  is  a  relatively 
inaccurate  procedure,  in  part  because  temperature  and  velocity 
are  correlated  in  turbulent  plumes,  and  errors  due  to  sparse 
sampling  of  turbulent  fluctuations  are  multiplied.  To  illustrate, 
temperature  and  velocity  are  known  at  24  m.  and  we  estimate 
the  plume  radius  to  be  about  2.5  m  at  this  height  (using  simple 
plume  relations.  f?  =  0.  IZ).  Multiplying  temperature  and  veloc¬ 
ity  at  each  point  and  taking  the  average  of  this  value,  we  can 
obtain  an  estimate  of  heat  flow  (using  C0  =  4200  J/  kg'’  °C 
density  1000  kg/ m3)  ranging  from  0.04  to  0.59  MW  for  the  three 
dives.  These  estimates  are  up  to  a  factor  of  50  times  lower  than 
of  any  of  the  temperature  or  velocity  heat  flow  estimates,  and 
they  illuminate  some  of  the  difficulties  in  using  the  technique  of 
averaging  flow  through  a  surface  for  turbulent  plumes  without 
sufficient  spatial  resolution 

Conclusions 

We  obtained  the  first  simultaneous  velocity,  temperature, 
and  conductivity  data  to  be  taken  in  the  convective  flow  eman¬ 
ating  from  a  hydrothermal  vent  field.  The  single  temperature, 
velocity,  and  conductivity  profiles  at  24  m  and  the  seven  simul¬ 
taneous  temperature  profiles  at  various  heights  indicate  that  the 
flow  field  forms  a  single  plume. 

The  flow  characteristics  and  surrounding  environment  are 
such  that  simple  plume  formulas  can  be  applied  to  predict  such 
plume  behaviors  as  maximum  height  and  total  heat  flow.  The 
nonlinear  effects  of  large  temperature  variations  have  been 
considered  here  as  far  as  they  relate  to  using  plume  theory  to 
predict  heat  flow  The  good  agreement  between  heat  flow  esti¬ 
mates  from  vertical  velocity  and  from  temperature  anomalies 
supports  the  heat  flow  estimate  of  3.7  ±  0.8  MW  for  this  vent 
field. 

Appendix:  Calculations  for  Plumfs  Wiih  Variable 
Nonlinear  Parameters 

The  primary  variables  of  the  plume  model  are  the  total  fluxes 
of  volume,  specific  momentum,  and  enthalpy,  integrated  over 
the  plume  cross  section: 

Volume  flux 

Q  =  jo2rrrw  dr  (10) 

Specific  momentum  flux 

M  =  /o2rrrw 2  dr  (II) 

Enthalpy  flux 

£  —  jl2rrrw  Ah  dr  (12) 

where  r  is  radial  distance  from  the  plume  centerline,  m  is  vertical 
velocity,  and  Ah  is  specific  enthalpy. 

The  time-averaged  tracer  concentration  (temperature,  den¬ 
sity,  or  velocity)  across  a  turbulent  plume  exhibits  essentially^ 


Gaussian  distributions  [Fischer  ei  al .  1979]  These  meanveloc- 
ity  and  concentration  profiles  are  self-similar  with  height,  ena¬ 
bling  tracer  distributions  to  be  expressed  only  in  terms  of 
maximum  (centerline)  values  and  some  measure  of  width  (the 


point  off-axis  where  the  vertical  velocity  ir  falls  to  e~'  =  0  37  of 
the  centerline  velocitv  W  for  Gaussian  distributions). 

iv  =  W  exp  (-r2/ R2)  (13) 

A/r  =  A exp  [  r?/(  (14) 

where  A  =  I  2[Fi.uherel  al. .  1979]  With  these  assumptions,  the 
integrated  flux  parameters  become 

(7=rrR2l  V  (15) 

M  —  K„rr  R2  H’*/2  (lb) 

£  =  A.ttA2/?2  «'A/w/(  I  +  A2)  (17) 

where  the  constant  coefficients  K„  =  0  89  and  A,  =  I  24  have 


been  included  to  account  for  the  additional  transport  of 
momentum  and  enthalpy  by  turbulent  fluctuations  rather  than 
by  mean  motions  [ Fischer  ei  al  .  1979). 

Following  Fischer  ei  al  [1979],  differential  conservation 


equations  for  Q.  M.  and  £  are  given  by 

dQ/dz  =  2rroRW 

(18) 

dM/dz  =  rrgtfA2/?2 

(19) 

dEJdz  =  0 

(20) 

Equation  1 18)  states  that  increase  in  specific  volume  flux  (7  is 
the  result  of  entrainment,  which  is  in  turn  proportional  to  the 
local  centerline  plume  velocity  where  the  constant  of  propor¬ 
tionality  <7  is  the  entrainment  coefficient.  The  value  of  cr  for 
buoyant  icts  lies  between  that  of  jets  (p |  =  0  0535)  and  plumes 
(p„  =  0  1)833).  Priestly  and  Ball  (1955 ]  found  that,  to  a  good 
approximation 

(7  —  o  Ki  (P|«t  Ppium«)(/7i/  Rip)  (21) 

where  Rip  is  the  plume  Richardson  number,  with  a  value  of 
0  557  [Fischer  el  al..  1979],  Ri  is  the  local  Richardson  number 
given  by 

Ri=  |[4(2>r)’2A2/(l  +  A2>Xg/M/»*]|”  (22) 

Equation  (19)  stales  that  the  vertical  buoyancy  force  acting  per 
unit  height  of  the  jet  equals  the  rate  of  change  of  vertical 
momentum  flux  The  quantity  6  is  the  specific  density  anomaly 
(p.  -  PmmdiPm  at  the  plume  centerline.  For  the  nonlinear  plume 
this  must  be  calculated  as  a  function  of  the  centerline  tempera¬ 
ture  (or  enthalpy). 

Equation  (20)  expresses  conservation  of  enthalpy  flux  E 
across  any  horizontal  plane  in  the  plume  The  maximum 
enthalpy  anomaly  A7;m.,  is  in  general  a  nonlinear  function  of 
the  plume  centerline  temperature  difference  ATm*..  Only  for 
small  values  of  A  Tm„  and  thus  constant  C„  is 

£  =  pCprrA2H  A  rm../(  I  +  A2)  (23) 


Little  £ t  al  Hyoroi  mermai  Vrsr  Pt  hme  Fi  ow 


The  system  of  equations  (( 1 8)-< 20))  can  be  solved  relatively 
easily  by  using  an  explicit,  first-order,  numerical  formulation 
Specific  volume  values  for  temperatures  of  NaCl  solutions 
greater  than  40° C  were  obtained  from  Chen  (1981}  and 
enthalpy  values  for  pure  water  from  A eenen  ei  al.  (1969). 
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ABSTRACT 

A  twelve  day  time  series  of  temperature  data  was  collected  from  a 
point  three  centimeters  above  a  diffuse  flow  area  at  a  Guaymas  Basin 
hydrothermal  site.  Using  concurrent  tidal  gauge  data  from  the  town  of 
Guaymas  it  is  shown  that  the  effects  of  tidal  currents  can  be  strong 
enough  to  dominate  the  time  variability  of  a  temperature  signal  at  a 
fixed  point  in  hydrothermal  flow  and  are  a  plausible  explanation  for  the 
variations  seen  in  the  Guaymas  Basin  temperature  data. 
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INTRODUCTION 


The  high  heat  flow,  high  permeability  and  low  sedimentation  of  the 
world's  Mid-Ocean  Ridge  system  support  hydrothermally  active  regions  of 
the  seafloor  in  localized  sites  encompassing  areas  as  large  as  10  in  . 

The  process  of  hydrothermal  circulation  transfers  heat  and  material  from 
the  lithosphere  into  the  oceans,  influencing  oceanic  composition  as  well 
as  creating  massive  ore  deposits  on  the  seafloor.  The  significance  and 
global  impact  of  this  process  are  determined  by  the  heat  and  mass  flux 
from  individual  sites,  the  number  of  hydrothermal  sites,  and  the 
variation  of  these  fluxes  over  time. 

Recent  efforts  have  focussed  on  characterizing  temporal  changes  in 
hydrothermal  activity.  Hydrothermal  vent  fields  evolve  continually,  in 
flow  rates,  size,  geometry  and  biological  composition  (Hekinian  et  al , 
1983;  Johnson  and  Tunnicliffe,  1985),  the  most  direct  evidence  being  both 
the  existence  of  numerous  extinct  fields  and  apparently  young  fields 
uncolonized  by  biota.  Possible  mechanisms  for  reducing  flow  include  the 
blocking  of  flow  channels  by  hydrothermal  precipitates,  exhaustion  of  the 
underlying  heat  source,  reduction  of  porosity  by  mineralization,  and 
obstruction  of  flow  paths  due  to  tectonic  activity.  On  the  other  hand, 
increase  in  flow  may  result  from  tectonic  activity  and  magma  injection. 
Periodic  oscillation  in  flux  may  be  produced  by  tidal  forcing  either 
indirectly  through  tectonic  activity  or  directly  through  pressure  effects 
on  the  flow.  Cataclysmic  changes  are  possible  based  on  models  of  heat 
transfer  near  the  magma-water  interface  of  the  convective  system  which 
predict  sudden  changes  in  flow  due  to  the  behavior  of  seawater  as  it 
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reaches  the  two -phase  region  of  pressure  and  temperature  (Cann  et  al. , 
1985)  . 

Quantifying  even  the  instantaneous  heat  and  mass  flux  from 
individual  sites  is  difficult  due  to  the  complexity  of  flow  systems 
typically  encountered  in  hydrothermal  fields.  In  a  single  field  the  ^low 
can  range  in  character  from  warm,  10°C  water  slowly  seeping  through  large 
masses  of  tube  worms  and  bacterial  mats,  to  hot,  350°C  water  jetting  from 
narrow  sulfide  chimneys  only  0.05  m  in  diameter.  The  most 
straightforward  flux  measurements  to  make  are  those  of  well-defined 
plumes  emanating  from  single  chimneys.  These  measurements  are  achieved 
by  obtaining  temperature,  flow  rate  and  exit  area  at  a  vent  orifice  using 
conventional  thermistors  and  mechanical  flow  meters  (Converse  et  al . , 
1984;  Macdonald  et  al.,  1980).  Only  slightly  more  complex  is  measuring 
temperature  and  velocity  above  a  multiple  outlet  sulfide  structure  and 
using  plume  theory  to  derive  heat  flux  from  the  entire  edifice  (Little  et 
al.,  1987). 

Much  more  difficult  to  quantify  is  the  diffuse  flow  emanating  from 
cracks,  fissures  and  biological  masses  at  velocities  of  less  than  about 
0.1  m/s  because  the  velocity  and  temperature  signals  are  low  and  unevenly 
distributed  over  a  large  area.  Diffuse  flow  is  visually  identified  by 
the  presence  of  shimmering  water,  upward  moving  particles,  hydrothermal 
mineral  deposits  and  an  abundance  of  biota.  Measuring  instantaneous, 
point  temperate 'e  is  reasonably  simple  although  typical  non-uniform 
distributions  requ5 re  many  measurements.  Velocity  is  quite  difficult  to 
obtain  because  the  flow  is  heterogeneous  over  scales  as  small  as  a  few 
centimeters  and  the  flow  velocity  is  low  and  on  the  same  order  of 
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magnitude  as  ocean  currents  (0.05-0.10  m/s)  which  can  dominate 
hydrothermal  fluid  related  flow. 

Measurements  of  temporal  changes  in  hydrothermal  activities  require 
some  means  of  obtaining  time -series  data  within  the  difficult  constraints 
of  deployment,  recovery,  power  consumption,  data  storage,  and  instrument 
reliability.  Little  (in  prep)  has  explored  the  use  of  passive  acoustic 
monitoring  for  this  purpose.  However,  until  this  method  is  proven, 
temperature  may  be  the  parameter  of  choice,  especially  for  low 
temperature  diffuse  flow. 

When  making  temperature  measurements  in  a  hydrothermal  environment, 
care  must  be  taken  to  separate  out  the  effects  of  changes  in  hydrothermal 
flow  and  temperature  from  changes  in  ambient  currents.  In  this  article 
we  examine  a  time -series  of  temperature  data  collected  near  a  diffuse 
flow  area  at  a  Guaymas  Basin  hydrothermal  site.  We  will  show  that  the 
effects  of  tidal  currents  can  be  strong  enough  to  dominate  the  time 
variability  of  a  temperature  signal  at  a  fixed  point  in  hydrothermal 
flow,  and  are  a  plausible  explanation  for  the  variations  seen  in  the 
Guaymas  Basin  temperature  data.  This  study  points  out  the  necessity  to 
measure  lateral  currents  whenever  temperature  measurements  made  in 
hydrothermal  vent  fields  are  used  to  examine  time  variability  of  venting. 

MEASUREMENTS  AND  RESULTS 

The  experiment  site  was  a  hydrothermally  active  area  located  in  the 
Southern  Trough  of  Guaymas  Basin  in  the  Gulf  of  California  at  27°00'N, 
111°24'W  (figure  Al)  where  vents  flow  through  a  thick  layer  of  pelagic 
sediment  before  breaking  out  onto  the  seafloor.  In  an  effort  to 
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characterize  temporal  variability  of  diffuse  hydrothermal  flow,  a 
temperature  probe  was  placed  0.03  m  above  the  sediment  surface  in  an 
approximately  5x5  m  zone  characterized  by  bacterial  mat,  whitish  mineral 
deposition  and  patches  of  shimmering  water.  A  few  meters  from  the 
thermistor  site  was  a  particularly  dense  bacterial  mat  in  and  above  which 
temperatures  were  greater  than  30°C  (the  upper  limit  of  the  temperature 
probe)  and  perhaps  as  high  as  230  °C  as  evidenced  by  charring  of  the 
recording  instrument's  wooden  base.  In  the  vicinity  of  the  study  area 
were  broad  areas  of  hydrothermal  discharge  covered  with  living  bacterial 
mat  whose  temperatures  were  measured  from  the  submarine  at  less  than  11°C 
above  ambient  at  a  depth  of  0.08  m  below  the  sediment  surface  and  78°C 
above  ambient  at  a  depth  of  0.50  m  in  the  sediments.  Numerous  270°C 
vents  were  associated  with  mounds  and  chimneys  that  rose  up  to  25  m 
above  this  sediment  pond  area  (Lonsdale  and  Becker,  1985). 

The  thermistor  was  deployed  with  a  submersible  and  left  in  place  for 
12.5  days  from  July  17,  1985  to  July  28,  1985,  during  which  time 
temperatures  were  recorded  every  5  minutes  for  a  total  of  3790  points. 

The  ambient  temperature  in  this  area  as  measured  during  descent  of  the 
instrument  package  was  2.85°C.  The  maximum  recorded  temperature  for  this 
deployment  was  20°C ,  which  occurred  during  deployment  probably  as  the 
probe  touched  or  came  very  close  to  the  sediment  surface.  We  will  use 
this  value  for  sediment  surface  temperature  in  calculations  that  follow. 

Tidal  high  and  low  water  surface  elevations  were  obtained  from  the 
town  of  Guaymas ,  Mexico,  approximately  100  km  away,  for  the  same  time 
period.  A  cubic  spline  is  subsequently  used  to  interpolate  this 
temporally  unevenly  spaced  data.  Tidal  currents  measured  at  10  m  above 
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bottom  in  the  troughs  of  Guaymas  Basin  have  reversing  northwest -southeast 
flow  with  velocities  of  up  to  0.12  m/s  (Lonsdale  and  Becker,  1985). 
Current  speed  was  not  recorded  at  the  vent  site  and  will  be  derived  from 
tidal  heights  recorded  at  the  town  of  Guaymas. 

The  most  prominent  feature  of  the  raw  temperature  data  (figure  A2) 
is  the  strong  semi-diurnal  oscillation,  which  exhibits  a  temperature 
range  from  a  minimum  of  3.05°C  to  a  maximum  of  4.87°C.  The  concurrent 
tide  data  (also  figure  A2)  show  that  the  dominant  tide  in  this  region  is 
diurnal.  Closer  examination  reveals  that  the  temperature  reaches  a  peak 
when  the  tide  is  at  a  minimum  or  maximum,  and  reaches  a  low  when  the  tide 
is  at  an  inflection-point. 


DISCUSSION 

From  the  observations  follows  the  hypothesis  that  temperature  is 
influenced  by  tidal  current  speed,  a  process  which  can  have  a  90°  phase 
lag  to  tidal  height  and  which  peaks  twice  as  often  as  tidal  elevation. 

To  test  this  hypothesis  we  will  compare  the  frequency  spectra  of 
temperature  and  tide  to  determine  whether  correlations  at  frequencies 
higher  than  diurnal  exist.  Following  this,  we  will  estimate  current 
speed  from  tidal  height,  use  this  in  a  simple  model  of  boundary  layer 
flow  and  compare  the  results  with  the  measured  temperature  data. 

To  compare  temperature  and  tidal  frequencies  under  the  supposition 
that  tidal  periods  are  twice  that  of  temperature,  we  have  presented  the 
frequency  spectra  with  the  tidal  period  axis  divided  by  two  (figure  A3). 
The  tidal  spectrum  was  calculated  after  cubic  spline  interpolation  (see 
figure  A2) .  The  temperature  spectrum  was  produced  by  first  dividing  the 
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data  set  into  fourteen,  256  point  segments  by  sub-sampling  the  full  set 
every  fourteen  points.  The  new  segments,  with  points  separated  by  70 
min. ,  are  used  to  produce  an  averaged  spectrum  with  a  minimum  resolvable 
period  of  2.33  hours.  The  power  levels  are  not  normalized  so  comparison 
only  can  be  made  of  the  peak  period  locations,  not  amplitudes.  There  are 
striking  similarities  between  the  two  spectra,  especially  at  the  12  and  6 
hour  periods  but  also  at  the  shorter  periods  such  as  4  and  3  hours, 
indicating  that  tidal  current  speed,  which  peaks  twice  as  often  as  tidal 
height,  may  be  influencing  temperature. 

Now  we  turn  to  a  boundary  layer  flow  model  to  examine  the  effects  of 
current  speed  on  temperature  at  a  fixed  point  above  a  heated  surface. 

When  a  current  passes  over  a  heated  surface,  a  thermal  boundary  layer  is 
produced  whose  thickness  increases  with  distance  downstream  from  the 
leading  edge  of  the  plate.  There  are  two  ways  that  a  vector  current  can 
influence  temperature  at  a  point  inside  a  thermal  boundary  layer.  First, 
if  the  heated  surface  is  asymmetrical  and  finite  and  the  sensor  is 
located  off  center  then  rotating  the  direction  of  the  current  vector 
while  maintaining  a  constant  current  magnitude  will  change  the  effective 
distance  between  the  sensor  and  leading  edge  of  the  surface.  As  boundary 
layer  thickness  is  a  function  of  this  distance,  a  current  rotation  will 
result  in  a  change  of  boundary  layer  thickness  and  hence  the  recorded 
temperature.  Such  temperature  variations  would  have  frequencies 
correlated  with  the  tides.  However,  it  would  be  pure  chance  if  the 
phases  of  variation  coincided  -  there  would  be  no  reason  to  expect  a 
temperature  low  at  a  tidal  inflection  point  since  the  timing  of  the 
variation  would  depend  only  on  sensor  position  within  the  hot  plate. 
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The  second  mode  of  current  influence  on  temperature  depends  on  the 
fact  that  thermal  boundary  layer  thickness  is  a  function  of  current 
magnitude,  and  temperature  at  a  fixed  point  above  the  surface  depends  on 
boundary  layer  thickness.  We  suggest  that  this  mode  is  operating  and 
temperature  variations  seen  are  due  to  the  fact  that  the  thermistor  is  in 
an  area  where  a  thickening  and  thinning  of  the  boundary  layer  is  caused 
by  changes  in  current  magnitude. 

To  examine  how  temperature  relates  to  flow  we  first  assume  a  general 
form  for  the  dependence  of  the  boundary  layer  thickness,  5,  on  flow 
velocity,  U: 

S  -  CU'H  (1) 

This  functional  relationship  is  consistent  with  a  thermal  boundary  layer 
development  in  a  flow  with  a  uniform  velocity  and  a  constant  diffusivity 
(Fischer  et  al.,  1979).  The  constant  C  will  be  determined  by  calibration 
with  the  temperature  and  velocity  measurements,  and  will  be  limited  in 
resolution  to  an  order  of  magnitude  because  surface  temperature  and 
bottom  current  speed  are  not  well  constrained.  These  approximations  do 
not  affect  comparisons  of  temporal  variability  so  much  as  comparisons  of 
absolute  temperatures. 

To  relate  temperature  to  the  position  of  a  measurement  point  within 
the  boundary  layer  we  assume  a  parabolic  for  the  temperature  profile: 

T(v)  °  Ts-(Ts-Ta) [2y  -  yjh  (2) 

S  6 

Where  Ts  is  the  surface  temperature  (-20°C) ,  Ta  is  the  ambient 
temperature  («*2.85°C),  and  y  is  distance  from  surface  ( -0.03m).  This 
profile  reflects  the  general  characteristics  of  the  temperature 


30 


variations  within  a  thermal  boundary  layer  (Schlicting,  1968).  A  simpler 
linear  functionality  yielded  less  satisfactory  results. 

The  flow  speed  over  the  heated  surface,  U,  was  assumed  to  be 
proportional  to  the  first  derivative  of  tidal  height  as  measured  at 
Guaymas .  To  obtain  a  functional  form  for  U,  the  interpolated  tide  data 
was  harmonically  decomposed  using  a  least  squares  sine  fit  such  that: 

h  =  KqS  A(  i)  sin[  9  ( i)  t  +  <&(i)]  (3) 

i 

where  A(i),  0(i)  and  4>{i)  are  the  amplitude,  frequency,  and  phase  of  each 
component  (see  table  1) .  The  flow  speed  U  is  then  modelled  as  being 
linearly  dependent  on  the  absolute  value  of  the  time  derivative  of 
equation  3  for  h: 

U  -  Kq  abs (  S  A(i)tf(i)cos[«(i)t  +  *(i)]  )  +  K1  (4) 

i 

The  constants  K  and  are  determined  to  match  estimated  values  of 
o  1 

the  maximum  and  minimum  current  speeds,  U  and  U  .  at  the  height  of 

r  max  min  ° 

the  sensor.  U  is  estimated  from  the  maximum  free  stream  velocity 
max  J 

measured  at  0.12  m/s  at  10  m  above  bottom  (Lonsdale  and  Becker,  1985). 
Using  a  logarithmic  velocity  profile  between  two  points  above  the  bottom: 

Ul  -  lnfzl) - ln(zo)  (5) 

U2  ln(z2) - ln(zo) 

where  zo  -roughness  height/25  (Townsend,  1980).  For  a  roughness  height 
of  0.05  m  we  estimate  zo  to  be  0.002m.  Using  zl-0.03m,  z2-10  m,  U2— 0.12 
m/s  in  equation  3  results  in  an  order  of  magnitude  estimate  for  U1  -  Umax 

-0.05  m/s.  This  value  of  U  may  be  used  directly  in  conjunction  with 

'  max  J  J 
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TABLE  1 


Period  Amplitude  Phase 

(hours)  (meters)  (degrees) 


4. 

.00 

0. 

.010 

166 

4. 

,14 

0. 

,002 

230 

6. 

,10 

0. 

,010 

206 

6. 

,21 

0. 

.002 

344 

6. 

.27 

0. 

,009 

273 

8. 

,18 

0. 

,013 

120 

11. 

.97 

0, 

.939 

230 

12. 

.00 

1. 

.032 

57 

12. 

.42 

0, 

.317 

13 

12. 

.66 

0. 

.164 

289 

23. 

.93 

1. 

.451 

3 

24. 

.07 

1. 

.184 

188 

25. 

.82 

0 

.204 

102 

26, 

.87 

0 

.160 

60 

327. 

.90 

0. 

.101 

40 

661 

.30 

0 

.140 

289 
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the  observed  temperature  minimum  Tm^n  —  3.05°C  and  equation  1  and  2  to 

_  3  3 y2 

calculate  C=7 . 5x10  m  '  /s  .  Using  the  same  equations  with  this  value 

of  C  and  T  -  T  -  4.87°C  yields  an  estimate  of  U  .  =  0.03m/s. 

Correspondingly,  the  values  of  Kq  -  9.27  and  -  0.03  m/s  may  be 

obtained.  In  addition,  with  the  above  values  for  C  and  U,  6  .  -  0.033  m 

mm 

and  6  -  0.043  m. 

max 

The  time  series  of  velocities  obtained  from  equation  4  may  now  be 
substituted  into  equation  1  and  2  to  obtain  an  estimate  of  temperature  at 
the  sensor  height  (figure  A4a  and  A4b) .  A  detailed  examination  shows 
extremely  good  correlation  between  predicted  and  actual  temperatures 
between  0  and  150  hours.  Following  this,  a  notable  difference  arises 
where  from  150-290  hours  maximum  temperature  peaks  decrease  and  the 
variations  get  slightly  out  of  phase.  This  is  probably  due  to  the  fact 
that  tides  can  have  a  rotational  component  which  is  not  considered  here, 
ie  they  are  not  exactly  reversing  and  the  minimum  current  increases 
during  these  periods.  This  causes  disagreement  between  the  estimated 
T(h) ,  which  is  forced  to  a  constant  minimum  set  by  scaling  to  v  .  . 


CONCLUSIONS 

The  important  measurements  of  temperature  variation  over  space  and 
time  for  the  purpose  of  estimating  heat  flux  and  flux  variability  from 
hydrothermal  vents  should  be  pursued,  however,  this  analysis  shows  that 
temperature  variability  measured  near  a  source  of  diffuse  hydrothermal 
flow  can  be  produced  by  tidal  currents.  A  time  series  of  temperature 
from  a  hydrothermal  area  in  the  Guaymas  Basin  can  be  explained  in  a  large 
part,  if  not  entirely,  by  invoking  boundary  layer  theory  and  tidal  cross 
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Figure  Al 

Map  showing  location  of  hydrothermal  vent  site  (small  box)  in  the 
Southern  Trough  of  Guaymas  Basin  in  the  Gulf  of  California.  Units  on 
contours  are  in  100  m.  Town  of  Guaymas,  marked  on  inset,  is  about 
100  km  away.  35 
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Figure  A3 

Frequency  spectra  from  raw  temperature  data,  solid  line,  and  cubic 
spline  fit  to  tides,  dotted  line.  The  frequency  axis  for  the  tidal 
spectrum  has  been  divided  by  two  to  show  good  comparison  of  estimated 
tidal  current  spectrum  with  temperature  spectrum.  The  y-axis  units  are 
useful  only  for  comparing  relative  power  levels  within  a  single 
spectrum. 
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Figure  Ai 

a)  Raw  temperature  (solid  line)  and  temperature  estimated  from 
boundary  layer  flow  model  and  tidal  currents  (dotted  line), 

shows  remarkable  similarity  in  both  relative  amplitude  and  temporal 
variation  for  the  first  half  of  the  thermistor  deployment. 

b)  Raw  temperature  (solid  line)  and  temperature  estimated  from 
boundary  layer  flow  model  and  tidal  currents  (dotted  line) 

from  second  half  of  data  record  shows  differences  in  relative  amplitude 
and  less  agreement  in  phase.  These  differences  are  probably  due  to  the 
use  of  the  first  order  approximation  that  currents  are  the  first 
derivative  of  height. 


SOUND  PRODUCTION  BY  HYDROTHERMAL  VENTS 

ABSTRACT 

Theoretical  examination  of  hot,  turbulent,  buoyant  jets  exiting  from 
hydrothermal  chimneys  revealed  acoustic  source  mechanisms  capable  of 
producing  sound  at  levels  higher  than  ambient  ocean  noise.  Pressure 
levels  and  frequency  generated  by  hydrothermal  jets  are  dependent  on 
chimney  dimensions,  fluid  velocity  and  temperature  and  therefore  can  be 
used  to  monitor  changes  in  these  parameters  over  time. 

A  laboratory  study  of  low  Mach  number  jet  noise  and  amplification  by 
flow  inhomogene 1 ties  confirmed  theoretical  predictions  for  homogeneous 
jet  noise  power  and  frequency.  The  increase  In  power  due  to  convected 
flow  Inhomogeneltles ,  however,  was  lower  In  the  near  field  than  expected. 

Indirect  evidence  of  hydrothermal  sound  fields  (Reidesel  et  al.  , 
1982;  Bibee  and  Jacobson,  1986)  showing  anomalous  high  power  and  low 
frequency  noise  associated  with  vents  is  due  to  processes  other  than  jet 
noise . 

Direct  measurements  of  hydrothermal  vent  sound  fields  show  that 
frequency  and  power  levels  are  comparable  to  ambient  ocean  noise  and  are 
consistent  with  jet  quadrupole  near-field  sound  which  has  been  amplified 
by  the  dipole  behavior  of  convected  flow  inhomogeneities.  This  near- 
field  amplification  Is  not  as  great  as  that  predicted  for  the  far-fleld 
but  Is  consistent  with  theoretical  considerations  for  near-fleld  dipole 
and  quadrupole  behavior. 
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INTRODUCTION 


The  global  transport  of  heat  and  chemicals  from  the  earth's  interior 
through  the  seafloor  and  into  the  biosphere  depends  on  the  movement  of 
water  through  oceanic  lithosphere.  At  the  crests  of  the  world's  mid¬ 
ocean  ridge  (MOR)  system,  cold  seawater  seeps  into  the  permeable, 
unsedimented  crust  and  is  heated  by  the  shallow  magma  chamber  and 
intruded  lavas  of  newly  formed  oceanic  lithosphere  in  the  process  called 
hydrothermal  circulation.  The  hot  water,  rich  in  minerals,  rises 
buoyantly  in  narrow  upflow  zones  from  crustal  depths  of  2-4  km 
(Macdonald,  1982)  and  exits  at  hydrothermal  vents  located  on  the  ridge 
axis  (figure  Bl).  Upon  mixing  with  the  cold,  ambient  ocean  bottom  water, 
metals  and  sulfides  precipitate,  depositing  zinc,  copper  and  Iron 
sulfides  and  forming  the  massive  sulfide  deposits  characteristic  of 
hydrothermal  vent  sites.  It  is  becoming  increasingly  apparent  that 
hydrothermal  systems  permeate  a  large  portion  of  the  Mid-Ocean  Ridge 
system  and  may  significantly  contribute  to  the  chemical  composition  of 
the  oceans  and  atmosphere. 

To  date,  active  hydrothermal  vents  have  been  found  (Hoagland  and 
Broadus ,  1987)  in  the  Pacific  Ocean  on  medium  to  fast  spreading  ridges: 
at  the  Galapagos  Rift  area  near  86°W  (Ballard  et  al . ,  1982),  along  the 
East  Pacific  Rise  (EPR)  at  20°S,  11°N  (McConachy  et  al.,  1986),  13°N,  and 
21°N  (Ballard  et  al.,  1981;  Ballard  et  al.,  1984;  Francheteau  and 
Ballard,  1983;  Hekinian  et  al . ,  1983b;  RISE,  1980),  In  the  Guaymas  Basin 
(Peter  and  Scott,  1985),  on  the  Explorer  Ridge  (Scott  et  al . ,  1984),  the 
Juan  de  Fuca  Ridge  (CASM,  1985;  Hammond  et  al.,  1984;  Crane  et  al . , 
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1985;  Normark,  1986),  the  Gorda  Ridge  (USGS,  1985),  the  Lau  Basin 
(Hawxins,  1986;  Craig,  1987),  and  in  the  Sea  of  Okhotsk  (Ocean  Science 
News,  1986).  Numerous  extinct  vent  sulfide  chimneys  have  also  been  -ound 
in  these  areas.  On  the  slow  spreading  Mid-Atlantic  Ridge  (MAR)  active 
vents  have  been  located  near  the  TAG  (Trans -Atlantic  Geotraverse)  area  at 
26°N  and  at  23°N  (Rona,  1984;  Rona,  1985). 

o 

Vent  fields,  ranging  in  area  from  10  to  10,000  m  ,  form  in  the 
narrow  active  volcanic  regions  of  ridge  crests.  The  fields  are 
characterized  by  sulfide  structures  up  to  15  m  high,  surficial  sulfide 
deposits  and  sedimentation,  and  are  usually  found  hosting  a  suite  of 
organisms  whose  primary  food  source  is  hydrogen  sulfide  from  vent  fluid. 
The  flow  regime  is  often  quite  complex  and  depends  on  the  subsurface 
geological  structure  as  well  as  sediment  cover.  The  highest  velocity  and 
hottest  flow  emanates  from  sulfide  chimneys,  with  typical  velocities  of 
1-3  m/s,  temperatures  of  250-350  °C,  and  exit  areas  of  10-30  cm^ 

(Converse  et  al.,  1984;  Little  et  al.,  1987).  This  type  of  flow 
originates  from  near  the  top  of  the  magma  chamber,  with  hot  fluid  moving 
upward  through  the  crust  in  a  relatively  unobstructed  path  (Cann,  1982). 
The  rest  of  the  vent  field  encompasses  areas  of  diffuse  flow  emanating 
from  cracks  and  fissures.  In  this  case,  the  high  temperature  fluid  has 
premixed  with  ambient  seawater  that  has  seeped  several  tens  of  meters 
into  the  subsurface.  The  fluid  is  therefore  cooler,  10-30°C,  and  has 
precipitated  most  of  its  sulfides  before  exiting  the  sea  floor. 

The  full  impact  of  hydrothermal  fluxes  of  heat  and  chemicals  on 
biogeochemical  cycles  can  only  be  assessed  by  determining  both  the  total 
spatial  extent  and  temporal  variability  of  hydrothermal  vents,  neither  of 


41 


which  are  presently  well  known.  Vent  sites  are  typically  found  by  towing 
cameras,  nephelometers ,  CTDs  (conductivity,  temperature,  depth 
recorders) ,  or  taking  hydrocasts  to  locate  the  plume  of  hydrothermal 
fluid  found  above  active  vents.  Measurements  of  instantaneous 
morphology,  chemical  composition,  mineralogy,  temperature,  flow  rates  and 
biota,  have  generally  been  done  with  a  submersible,  a  method  which 
results  in  short  time  samples  and  sparse  lateral  coverage  along  the  ridge 
axis  (Backer  et  al . ,  1985;  Hekinian  et  al.,  1983a;  Karsten  et  al.,  1984; 
Tivey  and  Delaney,  1985;  Hessler  et  al.,  1985).  In  addition,  both  the 
difficulty  in  location  vents  and  the  harsh  environment  they  comprise  have 
resulted  in  an  unusual  dearth  of  information  on  vent  activity.  The  high 
temperatures  and  particle  precipitation  at  vents  make  conventional  direct 
sensing  instruments  such  as  thermistors  and  flow  meters  Ineffective  for 
long-term  measurements  (none  have  operated  for  more  than  a  few  days; 
Converse  et  al.,  1984).  Also,  these  instruments  need  to  be  precisely 
placed  by  a  submarine  with  a  manipulator  arm  inside  an  active  sulfide 
chimney,  and  the  results  are  strongly  dependent  upon  the  position  of  the 
sensor  in  the  chimney. 

There  is  at  present  a  strong  scientific  need  for  methods  with  which 
to  measure  the  long-term  (days  to  years)  variations  in  hydrothermal  vent 
activity.  Calculations  of  vent  lifetime  based  on  such  Indirect  estimates 
as  vent  fauna  age  and  size  of  heat  source  in  the  ocean  crust  range  from 
one  to  104  years  (Converse  et  al.,  1984;  Hekinian  et  al . ,  1983b; 
Macdonald,  1982).  Hydrothermal  vent  fields  evolve  continually,  in  flow 
rates,  size,  geometry  and  biological  composition  (Hekinian  et  al.,  1983a; 
Johnson  and  Tunnicliffe,  1985),  the  most  direct  evidence  being  the 
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existence  of  numerous  extinct  fields.  Possible  mechanisms  for  reducing 
flow  include  the  blocking  of  flow  channels  by  hydrothermal 
mineralization,  exhaustion  of  underlying  heat  source  (Cann  and  Strens, 
1982),  reduction  of  porosity,  and  closing  of  flow  paths  by  tectonic 
activity.  On  the  other  hand,  flow  may  increase  as  a  result  of  tectonic 
activity  or  magma  injection.  Periodic  oscillation  in  flow  may  be 
produced  by  tidal  forcing  either  indirectly  through  tectonic  activity  or 
directly  through  pressure  effects  on  flow.  Cataclysmic  changes  are 
possible  based  on  models  of  heat  transfer  at  the  base  of  the  convective 
system  which  predict  sudden  changes  in  flow  due  to  the  behavior  of 
seawater  as  it  reaches  the  high  temperature  two-phase  region  of  pressure 
and  temperature  (Cann  et  al.,  1986).  Measurements  of  flow  variations  are 
needed  to  constrain  these  and  other  models  of  hydrothermal  processes . 

Long-term  monitoring  of  fluid  fluxes  within  the  severe  environmental 
conditions  at  vents  requires  the  development  of  a  remote  sensing 
technique  that  can  be  located  away  from  hydrothermal  fluid.  This  article 
will  examine  the  applicability  of  a  method  to  acoustically  monitor  fluid 
flow  by  recording  the  sound  field  generated  by  hydrothermal  vents.  The 
sound  field  of  a  hot,  turbulent  jet  is  known  to  be  related  to  fluid 
efflux  velocity,  fluid  temperature,  and  orifice  diameter  (Lighthill, 
1952).  The  monitoring  of  jet  sound  power  and  frequency  can  be  used  to 
detect  changes  in  fluxes  at  a  particular  site.  In  addition,  a  detailed 
understanding  of  sound  source  mechanisms  and  vent  geometry  can 
theoretically  allow  vent  sound  fields  to  be  used  to  determine  the 
absolute  velocity,  diameter,  and  temperature  of  fluid  discharge.  Long¬ 
term  monitoring  of  this  sound  field  will  give  information  on  the 
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variation  of  these  parameters  over  time  and  thus  provide  data  on  the 
variations  of  vent  fluid  fluxes  over  the  life  cycle  of  a  particular  vent 
field. 

The  use  of  passive  remote  acoustic  monitoring  with  available 
technology  alleviates  many  of  the  degradation  problems  associated  with 
instrumentation  placed  directly  in  hydrothermal  flow.  In  this  paper  we 
will  examine  evidence  for  hydrothermal  sound  generation  and  assess  the 
technique  of  determining  the  long-term  behavior  of  heat,  chemical  and 
mass  fluxes  from  the  active  hydrothermal  vents  found  along  Mid-Ocean 
Ridges  by  acoustically  monitoring  the  sound  produced  by  hot,  turbulent, 
hydrothermal  vent  jets.  Included  in  this  presentation  are  a  theoretical 
examination  of  sound  source  mechanisms,  a  laboratory  study  of  jet  noise, 
and  a  summary  of  results  of  noise  measurements  made  near  hydrothermal 
vent  sites  (Riedesel  et  al.,  1982;  Bibee  and  Jacobson,  1986;  Little,  in 
prep. ) 

THEORETICAL  SOURCE  MECHANISMS 

Sound  production  by  moving  fluids  occurs  through  the  conversion  of 
kinetic  energy  into  acoustic  energy  (Dowling  and  Ffowcs  Williams,  1983). 
The  kinetic  energy  of  a  moving  fluid  is  attributed  to  hydrodynamic 
motions  and  inertial  effects  and  depends  on  fluid  density  and  velocity, 
whereas  acoustic  energy  is  stored  in  the  elastic  forces  of  the  fluid  and 
depends  on  its  compressibility.  Pressure  variations  occur  in  all 
unsteady  flows;  sometimes  they  are  confined  to  the  flow  and  balance 
fluctuations  of  local  momentum,  and  sometimes  they  propagate  away  from 
the  flow  as  sound  (Ffowcs  Williams,  1969).  Propagating  and  non- 
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propagating  pressure  fields  (often  designated  sound  and  pseudo-sound, 

respectively)  are  sensed  in  the  same  way  by  hydrophones,  but  they  behave 

differently.  The  non-propagating  pressure  field  defines  the  near-field 

and  does  not  radiate;  it  is  a  hydrodynamic  pressure  field  in  which 

compressibility  effects  are  negligible.  The  propagating  pressure  field 

makes  up  the  far-field,  or  radiated  pressure  field,  which  behaves  as  a 

classical  acoustic  wave  whose  speed  depends  on  the  compressibility  of  the 

fluid.  The  non-propagating  pressure  level  is  proportional  to  the 

2 

dynamic  head  of  the  flow,  which  is  pV  ,  where  p  is  fluid  density  and  V  is 
fluid  velocity. 

The  basic  physical  processes  which  transfer  inertial  to  elastic 
energy  and  generate  a  radiating  pressure  field  in  a  fluid  primarily  do  so 
through  mechanisms  (figure  B2)  described  as  monopoles,  dipoles,  and 
quadrupoles  (or  more  generally,  multipoles).  The  monopole  radiator 
arises  through  volume  or  mass  fluctuations,  such  as  expanding  and 
collapsing  bubbles.  Also,  if  mass  flux  from  a  hydrothermal  vent  is  not 
constant,  but  varies  in  a  regular  way,  then  conditions  exist  for  monopole 
sound  production.  The  dipole  radiator  is  generated  by  external  force 
fluctuations  which  result  in  variations  of  momentum,  such  as  vibrations 
of  an  object  in  a  fluid.  Vibrations  of  the  chimney  induced  by  the  flow 
of  turbulent  hydrothermal  fluid  would  produce  dipole  sound  radiation. 

The  quadrupole  is  produced  by  fluctuations  of  the  rate  of  momentum  across 
a  fixed  surface,  such  as  turbulent  shear  stress  in  the  mixing  region  of  a 
jet.  In  the  mixing  region  of  turbulent  jets  where  shear  stress  is  high, 
kinetic  energy  is  converted  to  sound  through  fluctuations  in  momentum 
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flux.  Quadrupole  sources  also  can  be  considered  formed  from  the  opposing 
dipoles  of  vortices  found  on  the  edges  of  the  jet  (Powell,  1964). 

Generally,  in  dividing  multipole  pressure  fields  into  the  two 
regions  of  differing  behavior  (radiating  and  non- radiating) ,  the  near¬ 
field  is  chosen  as  less  than  one  wavelength  from  the  source  and  the  far- 
field  as  beyond  one  wavelength,  although  the  change  is  not  abrupt.  A 
single  monopole  radiates  a  pressure  wave  equally  in  all  directions.  The 
amplitude  decays  as  1/r,  where  r  is  radial  distance  from  the  source  and 
it  has  no  near-field  component.  A  dipole,  made  of  two  equal  but  opposite 
polarity  monopoles,  has  a  two-dimensional  radiation  pattern  and  is  a  less 
efficient  radiator  due  to  cancellation  effects  of  the  two  monopoles.  The 

near- field  of  a  dipole,  however,  has  a  higher  amplitude  than  that  of  a 

o 

monopole  of  the  same  source  strength,  and  decays  as  1/r  until  r-A  (A  « 

wavelength).  In  the  far-field,  the  amplitude  of  the  dipole  pressure  wave 

decays  as  1/r.  A  quadrupole  has  a  three-dimensional  radiation  pattern 

and,  being  formed  from  two  equal  but  opposite  dipoles,  is  the  least 

efficient  radiator.  However,  the  near-field  of  a  quadrupole  has  the 

3 

highest  amplitude  of  these  multipoles  and  it  decays  as  1/r  .  Its  far- 
field  decays  as  1/r  but  is  the  weakest  of  the  three  poles.  Monopole 
radiation,  if  present,  will  dominate  the  far-field  acoustic  signal, 
followed  by  dipole  radiation.  Quadrupole  radiation  will  be  significant 
only  at  high  fluid  exit  velocities  if  the  other  sources  are  not  present, 
or  in  the  near- field. 

At  hydrothermal  vents  it  is  possible  that  all  three  sources  are 
present  and  that  they  are  site  specific.  Monopole  radiation  can  be 
produced  by  cavitation,  boiling,  pulsating  exit  flows,  or  resonance  in 
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subsurface  cavities.  Dipole  radiation  will  arise  from  vibrations  of  the 
chimney,  interaction  of  the  turbulent  flow  with  a  rigio  surface,  or 
convection  of  flow  inhomogeneities.  Quadrupole  radiation  will  emanate 
from  the  shear  stress  produced  in  the  turbulent  mixing  region  of  the  free 
hydrothermal  jet.  Estimates  of  the  contribution  of  these  sources  to  a 
hydrothermal  acoustic  pressure  field  follow. 

Sound  production  by  hydrothermal  vents  is  theoretically  a  function 
of  fluid  exit  velocity  (V) ,  temperature  (T) ,  density  (p) ,  sound  speed 
(c) ,  orifice  diameter  (D) ,  and  sulfide  chimney  dimension  (L) .  The 
dependence  of  sound  on  these  parameters  varies  according  to  the 
particular  source  mechanism  involved.  These  parameters  are  difficult  to 
measure,  and  vary  from  site  to  site,  but  estimates  have  been  obtained  in 
several  experiments  (Converse  et  al.,  1984,  Little  et  al . ,  1987;  Tivey 
and  Delaney,  1986).  The  following  nominal  values  are  used  in  this 
discussion  of  acoustic  source  mechanisms: 


Velocity : 

V  -  2  m/s 

Temperature : 

T  -  350  °C 

Density: 

p  -  1000  kg/m 

Diameter: 

D  -  0.05  m 

Length : 

L  -  3  m 

Sound  Speed: 

c  -  1500  m/s 

Monopole 

i)  Pulsating  Exit  Flow 
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The  sound  generated  from  this  type  of  pulsating  flow  from  a  pipe  can  be 
treated  as  if  it  was  produced  by  the  motions  of  a  baffled  piston  (Ross, 
1976).  The  pressure  field  from  such  a  source  is  approximated  as: 

r(.r;  -  3^DT.f(o) 

(1) 

where  is  the  rms  pressure  fluctuations,  f(0)  -  1  for  A  >  0.5m.  The 
pressure  produced  by  the  mean  flow  of  a  jet  is  approximately  pV  where  V 
is  the  mean  exit  velocity.  If  we  assume  that  is  10%  of  the  mean  flow, 
then 

P£  -  0.1 pV2  -  40C  Pa  (2) 

and 

P(r)  -  12  Pa*m 

r  (3) 

ii)  Subsurface  Cavity  Resonance 

Hydrothermal  vents  are  generally  found  on  fresh  basalt  flows  which 
often  contain  numerous  drained  lava  lakes  and  collapse  pits.  Recent 
efforts  to  drill  into  zero-age  crust  at  a  hydrothermal  site  were  hampered 
by  the  presence  of  subsurface  cavities  (ODP,  1986).  These  cavities  are 
capable  of  resonating  and  producing  considerable  sound  if  excited  at  the 
proper  frequency.  Fluid  moving  in  and  out  of  a  cavity  provides  mass 
while  the  flexibility  of  the  cavity  walls  provides  the  spring  action  to 
sustain  tesonarce.  Flow  excitation  of  the  cavity  occurs  through 
interaction  of  turbulent  flow  with  the  walls  as  the  fluid  enter  the 
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volume.  If  the  dominant  frequency  of  turbulence  approaches  the  resonance 
frequency  of  the  cavity,  the  cavity  will  begin  to  pulsate,  which,  in 
turn,  will  strengthen  production  of  that  very  frequency  of  turbulence 
(Laufer  and  Yen,  1983).  Such  flow  excited  resonance  will  radiate  a 
strong  tonal  component  (Ross,  1976).  The  frequency  of  such  acoustic 
radiation,  being  dependent  on  cavity  size,  rigidity,  and  geometry,  is 
difficult  to  predict  since  very  little  information  on  subsurface  cavity 
structure  in  hydrothermal  systems  exits.  The  magnitude  of  radiated  power 
depends  on  how  close  the  bandwidth  of  the  excitation  frequency  is  to  the 
resonant  frequency  of  the  cavity.  This  type  of  sound  would  be  extremely 
site  specific  in  hydrothermal  systems. 

iii)  Cavitation 

Cavitation  is  the  formation  of  a  macroscopic  bubble  at  a  liquid- 

liquid  or  liquid-solid  interface  caused  by  a  local  drop  in  pressure 

(Ross,  1976;  Urick,  1975).  Such  local  pressure  drops  occur  when  a  moving 

fluid  is  forced  to  accelerate  around  a  bend  or  past  a  fixed  object.  The 

magnitude  of  this  pressure  drop  is  a  function  of  the  density  and  velocity 

2 

of  the  fluid  and  is  approximately  equal  to  pV  /2 ,  which  for  hydrothermal 
vents  is  about  2000  Pa.  Cavitation  inception  depends  on  the  existence  of 
submicroscopic  voids,  called  cavitation  nuclei,  in  the  liquid.  When 
pressure  outside  a  nucleus  drops  below  the  surface  tension  of  the  bubble, 
rapid  vaporization  can  occur  and  the  bubble  will  rupture,  expand  and  then 
collapse,  producing  significant  monopole  sound.  The  rates  of  growth  and 
collapse,  which  determine  acoustic  frequency  and  power,  are  affected  by 
the  interaction  of  the  pressure  and  velocity  fields  at  the  moving 
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boundary,  surface  tension,  evaporation,  dissolved  gas  content,  heat 
conduction,  viscosity  and  compressibility.  Dissolved  gas  in  the  fluid 
has  the  effect  of  cushioning  the  collapsing  bubble,  causing  multiple 
rebounds,  and  reducing  sound  production. 

The  p<»ak  frequency  of  cavitation  is  inversely  related  to  bubble 
radius  a  which  is  likely  to  be  small  given  the  depth  and  large  static 
pressure  at  vents: 

f  -  1_  /P 

2a  J  p  (4) 

where 


P  -  p(a)-p(«) , 


(5) 


p(®)  is  ambient  pressure,  p(a)  is  the  pressure  in  the  liquid  just  outside 
the  liquid- gas  surface,  and  p  is  ambient  fluid  density.  Low  frequency 
broad  band  noise  dominates  when  the  gas  content  of  the  bubble  is  high. 
Sharp,  high  amplitude  peaks  are  characteristic  of  cavitation  due  to  water 
vapor  only.  The  inception  of  cavitation  has  been  experimentally 
determined  to  depend  on  a  cavitation  parameter: 


K  -  (P(«) -Pv'Pg)/(0 . 5pV  ; 


(6) 


where  P  is  vapor  pressure  in  the  bubble,  and  P  is  gas  pressure  in  the 
^  & 

bubble.  If  K  is  greater  than  4,  cavitation  will  not  occur  (Ross,  1976), 
Hydrothermal  vents  have  been  found  at  depths  of  1500-3700  m  which 
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7  7 

correspond  to  ambient  pressures  of  1.5x10  -3.7x10  Pa.  The  vapor 
pressure  of  seawater  at  350°C  is  about  l.SxlO'7  Pa.  The  dissolved  gas 
content  will  effectively  increase  the  pressure  inside  the  bubble.  It  is 
therefore  possible  that  cavitation  occurs  at  shallow  vent  sites  where 
P(®)-Pv  <  8000  Pa.  However,  it  is  highly  unlikely  for  cavitation  to  be  a 
source  of  sound  at  greater  depths . 

iv)  Boiling 

Boiling  differs  from  cavitation  in  that  the  enlargement  of  bubbles 
occurs  due  to  an  increase  in  temperature  and  hence  vapor  pressure,  rather 
than  a  decrease  in  local  outside  pressure.  Sound  production  will  depend 
on  the  temperature  of  the  hydrothermal  fluid  and  the  depth  of  the  vent. 
Most  systems  will  not  exhibit  appreciable  sound  due  to  boiling.  Some  of 
the  shallower  vent  sites  could,  however,  experience  two-phase  separation 
and  generate  some  sound  through  growth  and  collapse  of  bubbles.  This 
process,  being  monopole  in  nature,  could  be  of  primary  importance  in 
hydrothermal  sound  generation,  if  present.  Recent  studies  north  of 
Iceland  (Jon  Olafsson,  personal  communication)  have  visually  documented 
boiling  at  discharge  sites  of  hydrothermal  vents  at  a  depth  of  100  m. 
These  vents  will  have  a  distinctive  acoustic  signature. 

Dipole 

i)  Chimney  Vibrations 

It  is  possible  for  turbulent  flow  through  sulfide  chimneys  to  induce 
vibrations  in  their  structures.  The  frequency  of  oscillation  depends  on 
the  stiffness  and  length  of  the  chimney.  The  amplitude  depends,  as  in 
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cavity  resonance,  on  mean  flow  velocity  and  on  how  close  the  forcing 
frequency  is  to  the  natural  frequency  of  the  chimney.  A  way  to  estimate 
pressure  from  a  vibrating  chimney  is  to  use  the  general  description  of 
dipole  resonance  coupled  with  a  frequency  estimate  from  theories  of 
structural  vibrations  for  circular  beams  (Ross,  1976).  The  dipole  source 
term  Do  due  to  a  fluctuating  force  F  at  frequency  f  is: 

D  -  F 

2*if  (7) 

Maximum  pressure  amplitude  for  dipole  radiation  is  given  by 
P  -  -(2fff)2Dn  -  2*ifF 

4?rrc  4jrrc  (8) 

The  frequency  of  resonance  in  a  circular  beam  is  given  by 
f  -  .I2g-l);rdc,b 

(9> 

where  c^  is  sound  velocity  in  the  beam,  is  beam  density,  and  d  is  beam 

3 

diameter  (-1  m) .  Using  a  value  of  cb-3000m/s,  -  4000  kg/m  and  m-1 , 

equation  9  results  in  a  frequency  of  40  Hz.  Force  F  is  estimated  by 
assuming  that  10%  of  the  available  pressure  from  fluid  flow  (4000  Pa  from 
above)  acts  upon  the  inside  wall  of  the  chimney  of  length  L  and  diameter 
0.05  m,  which  results  in  a  force  of  approximately  240N.  Using  this  in 
equation  8  produces  a  fluctuating  chimney  pressure  estimate  of 

P  -  3 . 2  Pa»m 

r  (10) 
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Lighthill 's  (1952)  formulation  of  the  aeroacoustic  problem  showed 

that  the  mixing  region  of  a  jet  could  be  equated  to  a  volume  of 

quadrupoles  with  strength  proportional  to  the  stress  tensor  in  the  moving 

fluid.  Such  a  mathematical  formulation  allows  radiated  jet  acoustic 

8 

power  levels  to  be  estimated  and  shown  to  be  proport:. onal  to  V  for  low 
Mach  number  flow  (Lighthill,  1952;  1954). 

The  frequency  of  jet  noise  depends  on  the  size  of  the  turbulent 
eddies  in  the  jet,  which  scale  with  jet  diameter  (Lighthill,  1963). 

Since  the  jet  expands  laterally  and  decreases  speed  with  increasing 
distance  from  the  orifice,  acoustic  power  and  frequency  from  a  given 
section  of  jet  will  both  decrease.  The  resultant  far-field  spectrum  will 
have  a  peak  frequency  near  0.8V/D  (Lush,  1971)  with  a  power  fall-off  of 
6-9  db  per  octave  (frequency  doubling)  below  the  peak  and  2-3  db  per 
octave  above  the  peak.  In  the  near  field,  the  frequency  spectrum 
produced  will  be  similar  but  depend  on  relative  proximity  to  either  the 
mixing  region  near  the  orifice,  which  produces  high  frequencies,  or  the 
fully  developed  turbulent  region  further  downstream  which  produces  lower 
frequencies . 

The  power  and  frequency  distribution  are  both  directional  in  jets, 
with  higher  frequencies  being  stronger  at  90°  to  the  jet  axis  and  lower 
frequencies  dominating  at  low  angles  (<45°)  to  the  jet  axis  (Lush,  1971). 

Sound  production  by  turbulent  jets  is  a  complicated  process  and 
extensive  theoretical  and  laboratory  work  has  examined  the  exact 
directionality,  frequency,  and  power  of  sound  from  low  to  high  Mach 
number  (M)  jets  (Lighthill,  1954;  Ffowcs  Williams,  1977;  Lush,  1971; 


Laufer  and  Yen,  1983;  Goldstein,  1984;  Powell,  1964;  Dowling  et  al . , 

1978;  Mankbadi,  1985;  Cohen  and  Wygnanski,  1987).  There  are  however, 
very  few  experimental  results  on  very  low  Mach  number  jets  (Jorgensen, 

.  3 

1961).  Hydrothermal  vents,  with  M-10  ,  fall  into  this  category. 

Therefore,  we  can  use  theoretical  studies  of  jet  radiation  efficiency  to 
predict  intensity  of  jet  noise  from  vents.  The  maximum  pressure  field  of 
low  Mach  number  jets  in  the  far-field  is  proportional  to  (Lush,  1971): 


P  « 

c  r 


at  a  peak  frequency  of 
f  «  V/D 

The  absolute  magnitude  of  the  near  and  far  pressure  fields  of  a 
turbulent  jet  generating  quadrupole  radiation  can  be  approximated 
(Ross,  1976;  Ffowcs  Williams,  personal  communication): 


p  -  10 ~2pV2 

-2  2  3 

p  -  10  (D/r) J 


-  40  Pa 


for  r  <  D  near-field 


-  5xl0’^/r^  Pa*m^  for  r  <  X  near-field 


p  -  10  (D/X)^D/r  -  4x10  ^/r  Pa*m  for  r  >  X  far-field 

where  X-cD/V  in  equation  15  returns  equation  12. 


(11) 


(12) 


by 


(13) 

(14) 

(15) 


Other  Effects 

i)  Pipe  Resonance 

Flow  through  a  hydrothermal  chimney  can  set  up  resonant  internal 
pressure  waves,  much  like  an  organ  pipe,  whose  frequency  depends  on  sound 
speed  in  the  fluid  and  on  the  length  of  the  pipe: 


f  -  nc/L 


(16) 
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where  n  is  an  integer.  Power  level  depends  on  how  similar  the  frequency 
of  the  driving  pressure  oscillations  are  to  the  resonant  frequencies  of 
the  pipe,  and  is  difficult  to  predict  in  hydrothermal  systems  since 
nothing  is  yet  known  about  variations  in  fluid  flow. 

ii)  Flow  Inhomogeneities 

If  a  turbulent  flow  is  of  non-uniform  density,  the  most  efficient 
source  is  not  the  velocity  quadrupole  term  but  a  dipole  order  term  whose 
radiated  power  scales  as  (Morfey,  1973).  This  is  due  to  the  fact  that 
accelerated  density  differences  found  in  heterogeneous  turbulent  flow  are 
effectively  force  fluctuations  resulting  in  variations  of  momentum,  which 
radiate  as  more  efficient  dipoles.  The  radiation  efficiency  of  this  type 
of  flow  is  represented  by  a  ratio  of  dipole  pressure  field  to  quadrupole 
pressure  field  (Morfey,  1973): 

Pd/Pq  -  c pA/V  (17) 

where  A  is  the  specific  volume  difference  between  jet  and  ambient  fluid. 
The  amplification  is  inversely  proportional  to  Mach  number  and  as  such, 
may  be  a  very  strong  effect  at  hydrothermal  vents.  For  hydrothermal 
fluid  at  350°C,  the  specific  volume  difference  is  approximately  0.5x10  ^ 
kg/m^ ,  and: 

Pd/Pq  -  4xl02  (18) 
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Using  the  pressure  calculated  above  for  quadrupole  radiation,  equation 
15,  yields  a  far- field  inhomogeneous  flow  pressure  level  of: 

-3 

P  -  1.5x10  /r  Pa*m  (19) 

The  complexities  and  variability  of  hydrothermal  vent  systems  allow 
for  a  wide  range  of  possible  sound  generating  mechanisms.  Sources  and 
nominal  sound  frequency  and  power  levels  are  tabulated  in  table  1  using 
values  for  vent  parameters  given  above.  Levels  indicate  that  some  vents 
may  produce  significantly  more  sound  than  others  based  on  geometry  and 
structure  alone.  In  addition,  the  range  of  velocities  seen  at  vents  will 
give  rise  to  an  additional  several  orders  of  magnitude  range  in  power 
levels  as  most  mechanisms  are  heavily  dependent  on  velocity.  This 
theoretical  analysis  suggests  that  monitoring  changes  in  vent  fields 
through  sound  is  feasible,  although  specific  sources  and  power  levels  may 
vary  considerably  depending  on  vent  site.  On  the  other  hand,  determining 
exact  values  for  velocity,  diameter,  length  and  temperature  will  depend 
so  heavily  on  chimney  dimensions  and  internal  structure  that  a  complete 
morphologic  description  of  the  vent  field  would  be  required  to  separate 
out  the  different  source  mechanisms. 

LABORATORY  EXPERIMENT 

In  order  to  determine  whether  theoretical  estimates  of  jet  noise 
production  through  quadrupole  radiation  and  amplification  by  flow 
inhorao gene i ties  were  valid  at  the  low  Mach  numbers  found  at  vents,  we 
conducted  a  simple  laboratory  investigation.  Our  measurements  were  made 
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in  the  near  field  since  no  experimental  investigations  exist  in  the 
literature  on  jet  noise  and  the  effects  of  flow  inhomogeneities  on  near¬ 
field  pressure  levels  in  flows  with  velocity  and  spatial  scales 
commensurate  with  hydrothermal  vents .  The  aim  of  this  study  was  to 
substantiate  theoretical  estimates  of  jet  sound  pressure  levels  and  use 
the  results  to  help  design  a  field  experiment  to  directly  measure 
hydrothermal  vent  sound  fields  (Little,  in  prep.). 

In  the  laboratory  we  concentrated  on  pure  jet  noise  production,  a 
process  dependent  on  quadrupole  radiation  for  homogeneous  fluids  and  on 
dipole  radiation  for  non-homogeneous  fluids .  Comparisons  were  made  of 
several  velocities  a,.’  nozzle  dimensions  for  cold  jets  (homogeneous)  to 
determine  effects  of  these  parameters  on  power  and  frequency.  Sound 
amplification  due  to  inhomogeneous  flow  was  investigated  using  both  hot 
jets  into  cold  ambient  fluid,  and  fresh  water  jets  Into  a  saline 
solution. 

Methods 

The  experimental  tank  holding  ambient  water  was  1.2  m  x  1.5  m  x  2.4 

m  and  constructed  of  wood  with  one  glass  face.  The  jet  orifice  was 

mounted  on  the  bottom  and  fed  through  a  pipe  running  down  the  side  of  the 

tank  (figure  B3).  Flow  through  the  jet,  driven  by  tap  pressure,  was 

-  5  3 

monitored  to  within  3x10  m  /s  by  a  flow  meter  in  series  with  the  pipe. 
To  reduce  the  effects  of  external  vibrations,  the  hydrophone  was 
suspended  in  the  tank  on  a  long  wire  cable.  The  hydrophone  itself  was  a 
deep  sea  hydrophone  (the  same  one  used  in  the  field  experiment  (Little  in 
prep.)),  15  cm  x  5cm,  and  was  calibrated  to  within  6  db  by  comparison  to 
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a  known  receiver.  Data  was  recorded  digitally  at  500  Hz  and  2000  Hz 
after  passing  through  an  appropriate  anti-aliasing  filter. 

Four  configurations  were  examined:  a  cold  (10°C)  jet  with  velocity 

of  3  m/s  passing  through  an  orifice  of  2.7  cm  diameter  into  cold  ambient 

water  with  sound  measured  at  a  range  of  49  cm  (to  center  of  hydrophone) ; 

a  cold  jet  with  a  velocity  of  7.5  m/s  passing  through  a  1.05  cm  diameter 

orifice  monitored  at  27  cm;  a  hot  (55°C,  p-985kg/m^)  jet  at  3  m/s  and 

orifice  diameter  of  2.7  cm  monitored  at  49  cm,  and  a  fresh  water  jet  with 

velocity  of  3  m/s  passing  through  an  orifice  of  2.7  cm  diameter  into  cold 

3 

salty  water  with  salinity  55  ppt  (density:  1041  kg/m  ),  monitored  at  a 
distance  of  49  cm. 

Results 

The  laboratory  environment  was  noisier  than  optimal  for  a  rigorous 
test  of  theory  (figure  B4)  but  the  results  achieved  here  agree  with 
theoretical  predictions  of  jet  noise  within  the  order  of  magnitude 
resolution  of  this  experiment  (summarized  in  table  2) .  Frequency  was  as 
predicted  by  equation  12  for  jet  noise  in  all  configurations.  Power 
measured  was  within  an  order  of  magnitude  of  that  predicted  by  the  near 
field  quadrupole  equation  14  for  the  two  cold  jet  arrangements.  The 
higher  velocity  jet  produced  somewhat  less  power  than  anticipated,  which 
may  be  due  to  the  dominant  wavelength  approaching  the  dimensions  of  the 
tank.  The  amplification  due  to  convected  flow  inhoraogeneities  was  less 
than  half  that  predicted  by  theory  (equation  17)  for  the  hot  jet,  and  ten 
times  less  than  predicted  for  the  saline  solution.  This  is  attributed  to 
the  fact  that  our  measurements  were  made  in  the  very  near  field  where 
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amplification  effects  of  flow  inhomogeneities,  causing  dipole  radiation, 
on  the  pressure  field  may  be  less  dominant  than  in  the  far-field. 

FIELD  OBSERVATIONS 

The  potential  difficulty  in  detecting  hydrothermal  vent  acoustic 
signals  lies  not  so  much  in  the  sensitivity  of  hydrophones  and  recording 
instruments,  as  in  the  intensity  and  variability  of  ambient  ocean  noise. 
Typical  ambient  ocean  noise  has  been  extensively  studied  (Urick,  1986; 
Burdick,  1984;  Hodgkiss  and  Anderson,  1980;  Wentz,  1962;  Akal  et  al . , 
1986)  but  absolute  levels  depend  on  the  exact  location  and  time  of 
sampling.  Noise  varies  considerably  depending  on  depth  of  water, 
weather,  shipping,  season,  local  topography  and  is  punctuated  by 
transient  sources  such  as  local  ships  and  biological  sources,  making 
estimates  for  specific  sites  difficult. 

The  lowest  frequency  band,  0.01-5  Hz  (figure  B5) ,  exhibits  high 
power  levels  and  is  dominated  by  microseisms:  low  frequency  pressure 
disturbances  caused  by  non-linear  interactions  of  ocean  surface  waves 
(Webb,  1984),  and  teleseismic  events. 

Power  in  the  band  from  5-100  Hz,  produced  by  distant  shipping,  is 
dependent  on  sound  that  has  travelled  tens  to  hundreds  of  kilometers  and 
is  strongly  influenced  by  wave  guide  propagation  effects  such  as  sound 
channelling  due  to  a  velocity  minimum  in  the  water  column.  Sound  which 
travels  great  distances  through  the  ocean  experiences  a  broadening  of  its 
spectral  peak. 

The  mid-band,  100-10,000  Hz,  is  a  function  of  local  sea  state  and 
wind  related  noise  caused  by  spray,  breaking  waves  and  falling  water 
droplets.  Sound  due  to  weather  within  a  few  kilometers  of  the 
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mt  :asureuient  si  te  will  dominate  pressure  levels  on  the  seafloor  In  this 
frequency  band. 

Finally,  the  high  frequency  band  above  10,000  Hz  is  dominated  by 
thermal  noise  caused  by  agitation  of  water  molecules. 

Transient  sources  also  contribute  to  noise  at  a  given  location. 

These  include  local  ships  and  submarines,  which  are  characterized  by  high 
energy  narrow  band  peaks,  often  including  harmonics  of  a  fundamental 
frequency,  anywhere  from  5  to  200  Hz.  Biological  sources  can  also 
contribute  to  the  noise  field,  with  whales  and  dolphins  capable  of 
producing  high  amplitude,  short  duration  sounds  at  frequencies  from  18- 
100  Hz  (Watkins,  1981). 

Published  Field  Data 
East  Pacific  Rise 

The  first  evidence  of  vent  sound  generation  came  from  data  collected 
near  hydrothermal  vents  at  21°N  on  the  EPR  where  an  array  of  ocean  bottom 
hydrophones  was  set  out  to  study  microearthquake  activity  at  an  active 
vent  site  (Riedesel  et  al . ,  1982).  In  comparing  two  hydrophones,  located 
300  m  and  2  km  from  the  vents,  it  was  noted  that  the  automatic  gain 
control  on  the  instrument  nearer  the  vents  set  itself  to  a  gain  level  16- 
64  times  less  sensitive  than  that  of  the  instrument  far  away.  This 
indicated  that  the  intensity  of  ambient  sound  near  the  vent  was 
consistently  louder  than  the  sound  further  away,  suggesting  that  the 
major  source  of  ambient  noise  in  this  area  was  the  hydrothermal  vents 
(Riedesel  et  al.,  1982).  In  addition,  the  month  long  record  of  gain 
level  shows  daily  variations  of  noise  level  that  are  correlated  between 
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the  two  hydrophones  (figure  B6).  An  inspection  of  a  time  series  record 
(figure  B7)  reveals  a  major  difference  in  amplitude  at  very  low 
frequencies  (.16  Hz).  These  frequencies  are  usually  associated  with 
surface  generated  microseisms  (Webb,  1984)  but  it  is  not  clear  why  the 
amplitude  should  be  so  different  at  a  receiver  separation  of  only  1.5  km. 
These  data  are  suggestive  of  a  low  frequency  vent  sound  source  but  the 
evidence  is  inconclusive. 

Juan  de  Fuca  Ridge 

Recently,  anomalous  high  ambient  acoustic  noise  was  observed  in  the 
caldera  of  Axial  Seamount  at  46°N  on  the  Juan  de  Fuca  Ridge  (Bibee  and 
Jacobson,  1986).  Noise  levels  from  2-30  Hz  varied  up  to  25  db  over  6  km 
as  measured  on  four  separate  receivers  (figure  B8)  and  it  was  suggested 
that  the  noise  source  was  at  or  near  the  sea  floor.  With  the  assumption 
of  spherical  spreading,  the  source  was  placed  at  400m  from  one  of  the 
instruments,  in  the  northern  part  of  the  caldera.  This  placed  it  within 
200m  of  a  low  temperature  vent  field  mapped  on  a  previous  expedition. 
Interestingly,  10  db  variations  in  sound  power  level  were  observed  over 
the  course  of  the  eight  day  record.  An  instrument  placed  several  hundred 
meters  from  the  known  high  temperature  vent  field  Ashes,  located  in  the 
southwestern  part  of  the  caldera,  recorded  no  evidence  of  anomalous 
acoustic  noise. 

ARGORISE  Towed  Hydrophone 

In  December,  1985,  we  attached  a  deep-sea  hydrophone  to  the  ARGO 
televiewer  system  during  the  survey  of  the  East  Pacific  Rise  (EPR) 


61 


between  latitudes  10°  and  12°N  (ARGORISE) .  Noise  was  recorded  in  real 
time  on  the  ship  via  ARGO's  conducting  cable.  Digital  samples  were  taken 
at  1000  Hz  every  hour  and  when  in  the  vicinity  of  hydrothermal  areas. 

The  vehicle  was  passively  towed  by  the  surface  ship  and  flow  noise  was 
minimized  by  reducing  the  ARGO's  horizontal  velocity  to  less  than  0.2 
kts.  Unfortunately,  a  heave  compensator  was  not  able  to  remove  all  the 
vertical  motions  induced  by  the  surface  swell;  thus,  only  about  10%  of 
the  recorded  data  was  uncontaminated  by  flow  noise.  Few  vigorous  black 
smokers  were  encountered  during  the  200+  hours  of  slaloming  down  the 
ridge  crest,  but  five  hydro thermally  active  areas,  as  defined  by  biota, 
were  crossed.  Comparison  of  sound  recordings  near  (within  several 
meters)  and  far  (further  than  1000  m)  from  these  areas  reveals 
consistently  elevated  sound  levels  near  hydrothermal  sites  (figure  B9)  at 
frequencies  generally  between  15  and  30  Hz.  Unfortunately,  absolute 
power  levels  were  unobtainable  but  noise  levels  near  vents  were  up  to  a 
factor  of  10  times  those  away  from  vents.  Contamination  from  vehicle 
flow  noise  prevents  more  robust  further  analysis  of  the  spectrum. 

Ashes  Vent  Field.  Juan  de  Fuca  Ridge.  Acoustic  Monitor  Experiment 

In  September,  1987,  high  quality  recordings  of  noise  within  a  few 
tens  of  meters  of  an  active  vent  were  made  using  two  hydrophones  emplaced 
by  the  submersible  Alvin  in  Ashes  Vent  Field,  Axial  Seamount,  on  the  Juan 
de  Fuca  Ridge  (45°55'N,  130°02'W)  (Little,  in  prep).  The  hydrophones 
were  suspended  38m  apart  on  a  vertica?  :able  with  the  lower  one  placed  2- 
3m  from  the  vent  orifice.  No  continuous  vent  signals  were  deemed 
significant  based  on  a  criterion  of  90%  probability  of  detection  and  5% 
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probability  of  false  alarm.  However,  a  small  signal  near  40  Hz,  with  a 

-4  2 

power  level  of  1x10  Pa  /Hz  was  noticed  on  two  records  taken  near 

Inferno  vent  (figure  BIO) .  This  vent  had  an  exit  diameter  of  4  cm  and 

velocity  of  2  ±  0.5  m/s  and,  using  equation  12  above,  the  expected  peak 

frequency  for  jet  noise  is  40+10  Hz,  which  is  consistent  with  the 

measurements.  Predicted  power  level  for  the  near- field  of  such  a  jet 

-  7  2 

is,  from  equation  14  with  r-2.3m,  1.7x10  Pa  /Hz.  The  power  level 

3 

suggests  the  occurrence  of  jet  noise  amplification  of  about  10  due  to 
convected  density  inhomogeneities.  Theoretically,  amplification  due  to 

_  3 

convected  flow  inhomogeneities  for  a  density  difference  of  0.5x10 

3 

kg/m  (corresponding  to  the  measured  exit  temperature  of  329°C)  can  be  a 

factor  of  up  to  1.6x10^  above  homogeneous  quadrupole  radiation,  or  a 

-2  2 

power  level  of  3x10  Pa  /Hz.  Measured  power  is  less  than  predicted  for 
near  field  radiation  with  such  amplification  and  we  attribute  the 
difference  to  the  fact  that  measurements  were  made  in  the  near -field 
where,  as  seen  in  the  laboratory,  amplification  effects  are  reduced. 

Summary  of  Field  Observations 

The  high  amplitude,  low  frequency  noise  observed  at  21°N,  if  it  is 
of  hydrothermal  origin,  would  have  to  be  produced  by  pulsations  of  flow 
or  large  cavity  resonance  (or  an  unidentified  mechanism)  since  sound 
source  processes  associated  with  turbulent  flows  would  generate  much 
higher  frequencies  than  observed. 

The  ARGORISE  data  is  consistent  with  black  smoker  inhomogeneous  jet 
sound  production  in  chat  the  frequencies  were  within  the  expected  range 
and  signals  were  above  ambient  only  within  a  few  meters  of  vents.  The 
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fact  that  no  anomalous  signals  were  detected  by  Bibee  and  Jacobson  (1986) 
near  the  Ashes  vent  field  is  consistent  with  Little's  Ashes  experiment  in 
which  the  hydrophone  had  to  be  within  several  meters  of  the  black  smoker 
to  detect  a  signal.  Even  at  that  distance  the  probability  of  detecting  a 
signal  was  lower  than  90%.  The  anomalous  sound  recorded  in  the  northern 
part  of  the  caldera  is,  in  light  of  these  results,  not  due  to  a  high 
temperature,  high  velocity  black  smoker.  It  may  be  more  related  to 
subsurface  hydrothermal/tectonic  interactions  where  the  cracking  of  rocks 
would  be  able  to  produce  higher  power  at  lower  frequencies  than  the 
simple  movement  of  water. 


CONCLUSIONS 

The  need  to  determine  long-term  variations  in  hydrothermal  fluid 
flow  prompted  this  study  of  sound  generation  at  hydrothermal  vents. 
Theoretical  examination  of  hot,  turbulent,  buoyant  jets  exiting  from 
hydrothermal  chimneys  illuminated  several  acoustic  source  mechanisms 

capable  of  producing  sound  at  levels  higher  than  ambient  ocean  noise. 
Among  these  were  fluctuating  flow,  pipe  resonance,  cavitation,  chimney 
vibrations  and  jet  noise  with  amplification  due  to  convected  density 
inhomogeneities.  Pressure  levels  and  frequency  generated  by  these 
mechanisms  are  dependent  on  chimney  dimensions,  fluid  velocity  and 
temperature  and  therefore  can  be  used  to  monitor  changes  in  these 
parameters  over  time. 

A  laboratory  study  of  low  Mach  number  jet  noise  and  amplification  by 
flow  inhomogeneities  confirmed  theoretical  predictions  for  homogeneous 
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jet  noise  power  and  frequency.  The  increase  in  power  due  to  corvected 
flow  inhomogeneities,  however,  was  lower  in  the  near  field  than  expected. 

Indirect  evidence  of  hydrothermal  sound  fields  (Reidesel  et  al., 
1984;  Bibee  and  Jacobson,  1986)  showing  anomalous  high  power  and  low 
frequency  noise  associated  with  vents  (21°N,  northern  Axial  caldera),  is 
due  to  processes  other  than  jet  noise.  Possible  hydrothermal ly  induced 
processes  include  low  frequency  fluid  source  pumping,  cavity  resonance, 
and  flow/rock  interactions  in  the  form  of  fracturing. 

Direct  measurements  of  hydrothermal  vent  sound  fields  (ARGORISE, 
Ashes  vent  field)  show  that  frequency  and  power  levels  are  consistent 
with  jet  quadrupoie  near-fiela  sound  which  has  been  amplified  by  the 
dipole  behavior  of  convected  flow  inhomogeneities.  This  near- field 
amplification  is  not  as  great  as  predicted  for  the  far- field  but  is 
consistent  with  theoretical  considerations  for  near-field  dipole  and 
quadrupoie  behavior. 

Hydrothermal  vents  generate  sound  but  source  mechanisms  at  vents 
will  be  site  specific,  depending  on  vent  geometry  and  velocity.  This 
coupled  with  great  variability  in  ambient  ocean  noise  means  that  some 
vents  will  generate  sound  which  is  more  easily  detectable  than  others. 
The  nominal  vent,  however,  produces  sound  at  a  level  close  to  ambient 
ocean  noise  and  as  such,  needs  a  complex  hydrophone  array  for 
quantitative  flow  monitoring. 
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Source  Mechanism 

Monopole 

Table  1 
Frequency 

Hz 

Pressure  Level 
Pa 

Power  at  3m 
Pa  /Hz 

Pulsation 

<  3000 

12.0/r 

16 

Cavity  Res. 

<  1000 

? 

<0(1)* 

Pipe  Res. 

n*250 

? 

<0(1) 

Cavitation 

>  1000 

? 

<0(1) 

Boiling 

? 

? 

? 

Dipole 

Chimney  Vibr. 

40 

3.2/r 

>1 

Flow  Inhoraog. 

far  field 

40 

2x10" 3/r 

>5x10' 7 

near  field 

40 

2/r3 

>6x10' 3 

Quadrupole 

far  field 

40 

4x10" 6/r 

>2x10' 12 

near  field 

40 

-3  3 

5x10  Vr 

>4x10' 8 

*  0()-  Of  order  of  magnitude  () 

Table  2 

Theoretical 

Measured 

Jet  Tank  Vel.  Dia. 

m/s  cm 

Freq.  Power  Amp.  Freq. 

Hz  Pa  /Hz  Hz 

Power  Amp. 

Pa  /Hz 

Cold  Cold  3.0  2.80 

88  2-8xl0'4 

70-90 

5 . OxlO'4 

Cold  Cold  7.5  1.05 

550  5-50xl0‘4 

450-650 

l.OxlO'6 

Hot  Cold  3.0  2.80 

88  2-8xl0'4  x 

55  70-90 

1 . 2x10" 3  2.5 

Cold  Salt  3.0  2.80 

88  2-8xl0'4  x 

380  70-90 

1.0x10' 3  2.0 
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Figure  B1 

Schematic  or  hydrothermal  circulation  in  oceanic  lithosphere  above 
magma  chamber  showing  isotherms  cai  d  by  convection. 


MONOPOLE 

DIPOLE 

QUADRUPOLE 


Figure  B2 

Important  jet  sound  generating  parameters,  including  chimney  length  L 
diameter  D,  exit  velocity  V.,  density  p. ,  and  ambient  density  p  ,  and 
monopole,  dipole,  and  quadripole  sourceJ mechanisms . 
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Figure  B3 

Laboratory  configuration  with  bottom  mounted  jet  and  cable  mounted 
hydrophone . 
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Power  PaA2/Hz  Power  PaA2/Hz 


Lab:  Cold  Jet  -  Ambient 


xlO3 


xlO  3  Lab:  Hot  Jet  -  Ambient 


Figure  BAa 

Mean  and  95%  confidence  limits  for  difference  In  sound  power  between 
ambient  and  a  laboratory  Jet  of  cold  water  at  3  m/s  through  a  2.7  cm 
diameter  nozzle.  Expected  peak  frequency  Is  about  85  He.  Variation 
near  60  and  120  Hz  Is  electronic  noise. 


Figure  BAb 

Mean  and  95%  confidence  limits  for  difference  In  sound  power  between 


ambient  and  a  laboratory  jet  of  hot  (130°  F)  water  Into  a  cold  tank 


3  m/s  through  a  2.7  cm  nozzle.  Expected  peak  frequency  ls~about~85 
Variation  near  60  and  120  Hz  Is  electronic  noise.  Hot  jet  sound  Is 
•  lavaUd  Above  cold  Jet  sound. 
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HYDROPHONE  PREAMP  GAIN  VS.  TIME 


Figure  B6 

Plot  of  gain  In  hydrophone  preamps  vs.  days  on  station  for  two 
hydrophones  near  the  black  smoker  vent  field  at  21°N  on  the  East 
Pacific  Rise.  The  upper  curve  Is  from  the  hydrophone  2km  from  the 
vents  and  the  lower  from  the  hydrophone  300m  from  the  vents .  The 
signals  are  correlated  and  indicate  louder  sound  near  the  vents  than 
far  away. 
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Figure  B7 

Time  series  of  data  from  two  hydrophones  near  21°N  on  the  EPR.  Larger 
amplitude  low  frequency  record  is  from  the  hydrophone  300m  from  the 
vents  and  the  other  from  the  hydrophone  2km  away. 
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DB  re  1  MICROPASCAL/SQRT(HZ) 


HYDROPHONE  NOISE  LEVELS  NEAR  85/07/16  2200 


FREQUENCY  (HZ) 


Figure  88 a 

Noise  levels  in  Axial  Volcano  from  Blbee  and  Jacobson's  1985 
experiment.  OBS  19  vas  in  the  northern  part  near  a  low 
temperature  vent  field.  OBS  3  was  near  Ashes  vent  field. 

Figure  B8b 

Site  map  of  above  experiment  with  Ashes  vent  field  in  south¬ 
western  part  of  caldera  and  CASH  low  temperature  field  in  northern 
part.  74 


Power  (PaA2/Hz) 


Inferno  1,1 1,  Hell  4,8,9 


Figure  BIO 

Low  noise  period  from  two  hydrophone  deployment  on  Ashes  vent 
field,  Axial  volcano  (Little,  in  prep).  Upper  solid  curve  is 
from  upper  hydrophone  on  first  deployment,  upper  dotted  is  from 
upper  hydrophone  on  second  deployment.  Lower  solid  curve  is 
from  the  lower  hydrophone  on  the  first  deployment,  lower  dotted 
curve  is  from  lower  hydrophone  on  the  second  deployment.  Notice 
in  particular  the  peak  on  the  lower,  first  deployment  at  39  Hz, 
due  to  jet  noise  at  Inferno  vent. 
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SOUND  ASSOCIATED  WITH  HYDROTHERMAL  VENTS: 


JUAN  DE  FUCA  RIDGE  FIELD  EXPERIMENT 


ABSTRACT 


High  quality  acoustic  noise  measurements  were  obtained  by  two 
hydrophones  located  3  m  and  40  m  from  an  active  hydrothermal  vent  on 
Axial  Seamount,  Juan  de  Fuca  Ridge,  in  an  effort  to  determine  the 
feasibility  of  monitoring  hydrothermal  vent  activity  through  flow  noise 
generation.  Most  of  the  noise  field  could  be  attributed  to  ambient  ocean 
noise  sources  of  microseisms,  distant  shipping  and  weather,  punctuated  by 
local  ships  and  biological  sources.  Water/rock  interface  wavos  were 
detected  which  showed  high  pressure  amplitudes  near  the  seafloor  and, 
decaying  with  vertical  distance,  produced  low  pressures  at  40  m  above  the 
bottom.  The  rapidity  of  pressure  fall  off  with  height  implies  a  local 
source . 

Detection  of  vent  signals  was  hampered  by  unexpected  spatial  non- 

stationarity  due  to  shadowing  effects  of  the  caldera  wall.  No  continuous 

vent  signals  were  deemed  significant  based  on  a  criterion  of  90% 

probability  of  detection  and  5%  probability  of  false  alarm.  However,  a 

-4  2 

small  signal  near  40  Hz,  with  a  power  level  of  1x10  Pa  /Hz  was  noticed 
on  two  records  taken  near  the  Inferno  black  smoker.  The  frequency  of 
this  signal  is  consistent  with  predictions  and  the  power  level  suggests 
the  occurrence  of  jet  noise  amplification  due  to  convected  density 
inhomogeneities. 
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INTRODUCTION 


At  mid-ocean  ridges  the  circulation  of  seawater  brings  heat  and 
chemicals  up  from  the  depths  of  oceanic  lithosphere,  releases  part  of 
this  load  onto  the  seafloor  and  injects  the  rest  into  the  water  column. 
The  full  impact  of  hydrothermal  circulation  on  ocean  composition,  global 
heat  flux  and  the  generation  of  economic  ore  deposits  has  not  been 
assessed  because  neither  the  full  spatial  nor  temporal  distributions  of 
venting  are  known.  Thirty  or  so  vent  sites  of  varying  size  have  been 
identified  both  in  the  Pacific  and  Atlantic  oceans  since  their  first 
discovery  in  1977  (Corliss  et  al.,  1979;  Hoagland  and  Broadus ,  1987). 
Although  instantaneous  measurements  of  hydrothermal  characteristics  such 
as  morphology,  temperature,  salinity,  chemistry,  heat  flux,  and 
biological  composition  have  been  taken  at  many  of  them  (Converse  et  al., 
1982;  Hekinian  et  al . ,  1983;  Hessler  and  Smithey,  1983;  Ballard  et  al., 
1984;  Hammond  et  al.,  1984;  Norraark  et  al . ,  1986;  Rona  et  al.,  1986; 
McConachy  et  al . ,  1986;  Tivey  and  Delaney,  1986;  Little  et  al.,  1987; 
Craig  et  al . ,  1987),  long-term  measurements  of  any  type  greater  than  a 
few  days  are  extremely  rare  (Johnson  and  Tunnicliffe,  1985;  Little  et 
al.,  submitted).  This  deficiency  is  due  to  the  inaccessibility  and 
severe  environmental  conditions  found  at  vents  eg.  high  temperatures, 
high  pressures  and  reactive  chemicals.  An  understanding  of  the  local 
scale  fluid  flow  characteristics  and  changes  in  these  over  time  is 
required  to  develop  models  of  hydrothermal  processes  ranging  from 
subsurface  water- rock  interactions  (Cann  and  Strens,  1986)  to  biological 
dispersion  and  growth  (Hessler  et  al.,  1985;  Van  Dover,  1986),  The 
current  scientific  need  to  remotely  monitor  long-term  changes  in  flow 
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velocity  has  justified  a  feasibility  test  of  a  passive  acoustic  monitor 
at  a  high  temperature  vent  field.  This  method  utilizes  the  sound 
generated  by  moving  fluid  at  hydrothermal  chimneys,  the  frequency  and 
amplitude  of  which  depend  upon  the  fluid  velocity  and  density,  orifice 
diameter  and  chimney  structure  (Lighthill,  1952;  Dowling  et  al . ,  1977; 
Ffowcs  Williams,  1969;  Morfey,  1973). 

In  theory  an  exact  description  of  vent  sound  could  be  used  to 
determine  flow  velocity,  orifice  diameter  and  possibly  fluid  density 
(Little  and  Stolzenbach,  submitted).  In  practice,  however,  no 
unequivocal  proof  exists  that  vents  generate  sound  at  levels  that  are 
detectable  in  the  deep  ocean,  although  some  data  is  available  that 
suggests  this  is  so  (see  review  in  Little  and  Stolzenbach) .  No 
measurements  of  vent  sound  have  been  made  of  adequate  quality  to 
determine  source  mechanisms  and  permit  estimation  of  vent  parameters.  In 
this  paper  an  experiment  is  described  that  was  designed  to  determine  the 
feasibility  of  detecting  hydro thermally  generated  sound  in  the  ocean.  In 
September,  1987,  high  quality  recordings  of  noise  within  a  few  tens  of 
meters  of  an  active  vent  were  made  using  two  hydrophones  emplaced  by  the 
submersible  Alvin  in  Ashes  Vent  Field,  Axial  Seamount,  on  the  Juan  de 
Fuca  Ridge  (45°55'N,  130°02'W).  The  results  of  this  experiment  are 
presented  in  this  paper. 


REVIEW  OF  AMBIENT  OCEAN  NOISE 

The  potential  difficulty  in  detecting  hydrothermal  vent  acoustic 
signals  lies  not  so  much  in  the  sensitivity  of  hydrophones  and  recording 
instruments,  as  in  the  intensity  and  variability  of  ambient  ocean  noise. 
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Typical,  ubiquitous,  deep-water  ambient  ocean  noise  spectra  can  be 
separated  into  four  frequency  bands  caused  by  four  different  source 
mechanisms,  figure  Cl,  (Urick,  1975;  Burdick,  1984;  Wentz,  1962). 

The  lowest  frequency  band,  0.01-5  Hz,  exhibits  high  power  levels  and 
is  dominated  by  microseisms:  low  frequency  pressure  disturbances  caused 
by  non-linear  interactions  of  ocean  surface  waves  (Webb,  1984),  and 
teleseismic  events.  Temporal  variations  are  as  rapid  as  a  few  minutes 
for  interface  waves  travelling  near  the  sea  floor,  and  as  long  as  a  few 
hours  for  variations  in  sea  state  and  swell. 

Power  in  the  band  from  5-100  Hz,  produced  by  distant  shipping,  is 
dependent  on  sound  that  has  travelled  tens  to  hundreds  of  kilometers  and 
is  strongly  influenced  by  wave  guide  propagation  effects  such  as  sound 
channelling  due  to  a  velocity  minimum.  Travelling  great  distances 
through  the  ocean  broadens  the  spectral  peak  and  eliminates  waves 
travelling  outside  a  few  degrees  of  horizontal.  Time  variations  in  this 
band  are  slow,  on  the  order  of  hours  to  days  and  depend  on  changes  in 
ship  traffic  and  large  scale  temperature  and  salinity  structure  (in  the 
sound  channel) . 

The  mid-band,  100-10,000  Hz,  is  a  function  of  local  sea  state  and 
wind  related  noise  caused  by  spray,  breaking  waves  and  falling  water 
droplefs.  Sound  due  to  weather  within  a  few  kilometers  of  the 
measurement  site  will  dominate  pressure  levels  on  the  seafloor  in  this 
frequency  band.  Changes  with  time  in  power  level  are  on  the  order  of 
hours  to  days  and  are  dependent  on  weather  patterns . 

Finally,  the  high  frequency  band  above  10,000  Hz  is  dominated  by  noise 
caused  by  thermal  agitation  of  water  molecules. 
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Transient  sources  also  contribute  to  noise  at  a  given  location. 

These  include  local  ships  and  submarines,  which  are  characterized  by  high 
energy  narrow  band  peaks,  often  including  harmonics  of  a  fundamental 
frequency,  anywhere  from  5  to  200  Hz.  Biological  sources  can  also 
contribute  to  the  noise  field,  with  whales  and  dolphins  capable  of 
producing  high  amplitude,  short  duration  sounds  at  frequencies  from  18- 
100  Hz  (Watkins,  1981). 

Noise  in  the  ocean  is  typically  non- stationary  both  in  time  and 
space  (Burdick,  1984;  Hodgkiss  and  Anderson,  1980).  The  noise  field 
varies  considerably  over  time  due  to  the  rich  variety  of  sources  (Akal  et 
al . ,  1986)  and  measurements  made  hours  apart  can  show  striking 
dissimilarity  (as  will  be  seen  below).  Sound  recorded  at  near  bottom 
hydrophones  is  subject  to  additional  propagation  effects  of  local 
topography,  seafloor  heterogeneity,  and  the  interface  of  the  ocean- 
seafloor  boundary.  These  effects  can  produce  severe  spatial  variability 
in  the  sound  field  over  short  distances  through  diffraction,  scattering 
and  exponential  decay  of  interface  waves  with  distance  from  a  boundary. 

In  the  face  of  these  difficulties ,  and  to  make  the  study  of  noise  in  the 
ocean  tractable,  the  necessary  assumption  of  stetionarity  is  often 
validated  by  making  comparisons  over  closely  spaced  times  and  distances. 


SITE  DESCRIPTION 

Ashes  Vent  Field  (Hannington  and  Scott,  1985;  Hammond  et  al,  1986) 
is  located  75  m  from  the  south-western  wall  of  Axial  Seamount  on  the  Juan 
de  Fuca  Ridge  at  45°55'N,  130°02'W  (CASM,  1983;  CASM,  1985;  Ernbley  et 
al .  ,  1988),  figure  C2 .  The  site  is  noted  for  its  smooth  floor  and 
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absence  of  cracks  and  fissures,  which  makes  it  an  ideal  spot  in  which  to 
work  with  the  submersible  and  emplace  a  vertical  instrumented  cable.  The 
caldera  floor,  approximately  4  km  from  east  to  west  and  10  km  north  to 
south,  is  1540  m  below  the  sea  surface  which  nominally  places  it  in  the 
deep  sound  channel  for  this  latitude  (Burdick,  1984).  In  the  Ashes  Vent 
Field,  there  are  two  main  black  smokers  with  accompanying  sulfide  and 
anhydrite  chimneys,  "Hell"  and  "Inferno",  separated  by  35  m  (figure  C3) . 
In  addition  to  these  there  are  several  lower  temperature  and  velocity 
white  smokers  and  many  patches  of  diffuse  flow  within  the  60x60  m  area  of 
hydrothermal  activity.  At  the  time  of  acoustic  sampling,  the  chimney  of 
Inferno  was  3  m  high,  topped  by  one  site  of  black  smoke  efflux  sampled  at 
a  temperature  of  326°C  with  visually  determined  velocity  of  1-2  m/s  and 
diameter  of  4  cm.  There  were,  in  addition,  several  other  sites  of  clear 
fluid  discharge,  one  located  20  cm  from  the  black  smoker,  with  a 
temperature  of  126°C,  velocity  of  less  than  1  m/s,  and  diameter  of  1  cm. 
The  others  were  located  near  the  bottom  of  the  edifice.  The  chimney  of 
Hell  was  2  m  high  and  hosted  a  single  black  smoker  at  its  top  with  an 
orifice  diameter  of  5  cm  and  a  visually  estimated  velocity  of  0.5-1. 5 
m/s . 

INSTRUMENTATION  AND  DEPLOYMENT 

The  hydrothermal  acoustic  monitoring  instrument  consisted  of  two 
hydrophones  suspended  on  a  cable  beneath  a  float  (figure  C4)  and  attached 
to  a  microprocessor  controlled  digitally  recording  system  (Mellinger  et 
al .  ,  1986).  Designed  for  detecting  an  unknown  vent  acoustic  signature, 
the  system  had  a  16  bit  analog  to  digital  converter  and  a  programmable  1- 
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10-100-1000  gain  amplifier  resulting  in  a  dynamic  range  of  156db.  One 
tunable  8-pole  Butterworth  filter  for  each  channel  provided  anti-aliasing 
for  sampling  rates  of  300  and  2400  Hz.  Calibration  to  absolute  sound 
power  levels  was  obtained  by  comparison  to  a  known  receiver  and  is 
accurate  to  within  6  db  for  the  bandwidth  15-1200  Hz  and  15db  for  the 
bandwidth  1-15  Hz.  Calibration  of  relative  response  between  the  two 
channels  revealed  less  than  a  4  db  difference  for  all  frequencies  of 
interest  (1- 1200Hz)  in  this  experiment. 

The  deployment  scheme  was  designed  to  accommodate  three  major 
constraints:  minimizing  the  _cs  of  the  submersible,  recording  without 
ships  or  submarines  in  the  vicinity,  and  limited  memory  storage 
capability.  To  this  end,  the  procedure  was  to  drop  the  instrument  from 
the  ship  onto  the  vent  site  using  an  acoustic  navigation  net,  dive  on  the 
mooring  with  a  submersible,  pick  it  up,  and  move  it  to  within  5  m  of  an 
active  chimney.  The  instrument  then  turned  on  that  evening  at  a  preset 
time  and  recorded  through  the  night  c.fter  the  submarine  and  ship  had  left 
the  area.  The  following  day,  the  instrument  was  acoustically  released 
and  recovered  by  the  surface  ship.  The  data  was  then  transferred  to  a 
portable  computer  and  the  instrument  readied  for  further  deployments. 


DATA  DESCRIPTION 

Two  fully  successful  hydrothermal  acoustic  monitor  deployments  were 
obtained  on  Atlantis  II/Alvin  voyage  118,  leg  21.  For  the  first 
deployment,  Alvin  dive  #1917,  the  lower  hydrophone  was  placed  2  m 
horizontally  and  1  m  vertically  from  the  north-east  side  of  black  smoker 
Inferno.  The  upper  hydrophone  was  39.5  m  above  the  lower  one.  The 
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instrument  turned  on  at  23:00  on  9/23/87  local  time  and  recorded  two 
consecutive  sets  of  data  every  hour  until  03:00  on  9/24/87.  The  first 
set  was  9216  points  long  taken  at  2400  Hz  with  anti-alias  filter  at  800 
Hz.  The  second  set  was  8192  points  long  recorded  at  300  Hz  with  anti¬ 
alias  filter  at  100  Hz.  A  total  of  14  sets  were  obtained,  7  at  each  of 
the  two  frequencies. 

On  the  second  deployment  the  lower  hydrophone  was  positioned  2  m 
horizontally  and  2  m  vertically  above  black  smoker  Hell,  35  m  from 
Inferno,  on  dive  #1923.  The  upper  hydrophone  was  38.9  m  above  the  lower. 
The  instrument  began  recording  at  18:00  on  9/29/87  local  time  and  sampled 
once  an  hour  until  06:00  on  9/30/87.  Each  sample  was  8192  points  long 
taken  at  300  Hz  with  anti-alias  filter  at  100  Hz,  for  a  total  yield  of  13 
sets . 

GENERAL  DATA  REDUCTION 

The  output  of  a  hydrophone  placed  in  the  ocean  results  from  a 
combination  of  system  noise  and  pressure  fluctuations  in  the  ocean.  Due 
to  the  fact  that  all  the  sources  of  fluctuations  are  not  completely 
known,  only  a  statistical  description  is  permitted,  based  on  observations 
over  an  extended  time  period.  In  producing  our  statistical  description 
of  the  hydrophone  output  and  in  analyzing  the  sound  we  will  use  several 
assumptions.  First,  we  will  assume  that  the  noise  is  temporally 
stationary  over  our  individual  sample  periods.  This  enables 
interpretation  of  fourier  transforms  of  the  time  series  as  representative 
of  the  distribution  of  power  over  frequency.  To  help  validate  this 
assumption,  examinations  of  the  time  series  will  be  used  to  eliminate 
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impulsive  signals  and  gross  differences  with  time.  The  second  assumption 
is  that  the  signal  from  the  vent  is  constant  over  the  duration  of  the 
experiment.  Visual  observations  of  the  smokers  support  this  assumption 
as  the  flow  velocity  appeared  to  remain  constant  from  dive  to  dive. 

Power  spectra  are  calculated  for  this  analysis  using  the  Welsh 
method  (Oppenheim,  1975).  The  8192  point  records  were  divided  into  16, 
512  point  sets,  each  set  multiplied  by  a  Hanning  window,  and  used  with  an 
FFT  algorithm  to  compute  the  power  spectral  densities  which  are 
subsequently  averaged  together.  This  gives  a  power  estimation  accuracy 
of  1/VN  P  -  1/4  P  where  N  is  the  number  of  sets  and  P  is  the  power  level 
at  a  given  frequency.  The  mathematical  spectral  estimates  are  thus 
accurate  to  +  0.125xP.  The  calibration  between  the  two  hydrophones  is 
better  than  a  factor  of  1.3xP,  since  the  measured  differences  (for  which 
corrections  have  been  made  in  the  plotted  spectra)  were  about  4  db  or  a 
factor  of  1.5  for  the  frequencies  in  the  bandwidth  of  this  system.  The 
frequency  resolution  of  1/T,  where  T  is  the  time  length  of  a  set,  is  0.6 
Hz  for  the  300  Hz  samples,  and  4.7  Hz  for  the  2400  Hz  samples. 

In  addition  to  comparisons  of  simple  power  spectra,  the  coherence 
(coherence  squared)  and  phase  between  the  two  channels  will  be  examined. 
The  coherence  level  reveals  the  amount  of  signal  common  to  both 
receivers  while  the  phase  differences  at  a  given  frequency  provides 
information  on  the  angle  of  incidence  of  the  incoming  wave.  The 
direction  of  wave  travel  can  be  obtained  by  looking  at  linear  trends  in 
phase  as  a  function  of  frequency.  Phase  is  calculated  here  such  that  if 
a  broadband  signal  impinges  on  the  array  at  an  angle  8  from  above,  the 
phase  will  be  a  linearly  increasing  trend  in  frequency  whose  slope  is 
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dependent  on  6  and  receiver  separation.  If  the  waves  arrive  from  below, 
phase  will  decrease  with  increasing  frequency. 


ANALYSIS  OF  NOISE  FIELD 

Two  power  spectra  encompassing  bandwidths  of  1-800  Hz  from  the  first 
deployment  (Inferno)  and  1-100  Hz  for  the  second  (Hell),  presented  in 
figure  C5  and  C6 ,  represent  two  of  the  lower  noise  periods  of  the 
experiment.  The  most  obvious  feature  in  these  spectra  is  the  difference 
in  power  levels  between  the  two  hydrophones  in  the  bandwidth  10-200  Hz, 
with  the  upper  hydrophone  receiving  more  power  than  the  lower.  Further, 
a  comparison  with  ambient  ocean  noise  curves  (figure  Cl,  from  Urick) 
reveals  a  marked  similarity  in  power  levels  and  spectral  shape.  We 
hypothesize  that  most  of  the  noise  recorded  in  the  caldera  is  due  to  the 
ambient  ocean  noise  sources  described  above  and  that  the  major  difference 
in  signal  level  is  due  to  a  shadowing  and  reflecting  effect  by  the  nearby 
caldera  wall.  The  following  analysis  examines  this  hypothesis  in  detail. 

High  frequency:  100-800  Hz 

Sound  in  this  frequency  band  is  dominated  by  sea  state  and  wind 
force  and  will  experience  changes  on  a  time  scale  of  several  hours  to 
days.  An  examination  of  the  time  history  for  the  Inferno  deployment, 
figure  C7 ,  reveals  power  variations  of  about  a  factor  of  three  for  the 
bandwidth  200-800  Hz  over  the  course  of  the  deployment.  The  ship's  log 
records  a  constant  sea  state  of  3-5  feet  and  wind  increasing  from  11  to 
21  knots  over  the  deployment.  This  is  consistent  with  the  general  trend 
of  increasing  power  seen  in  figure  C7 . 
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Weather  related  sound  propagates  downward  from  the  sea  surface  with 
sea  conditions  directly  overhead  exerting  the  most  influence  on  sound  at 
the  seafloor  below.  Since  sound  waves  are  travelling  to  the  receiver 
array  from  above,  one  would  expect  little  or  no  caldera  wall  effects  on 
this  bandwidth  if  the  sources  are  indeed  weather  related.  Figure  C8 
supports  this  as  there  is  very  little  difference  in  power  between  the  two 
receivers  above  about  250Hz.  Although  the  coherence  is  generally  weak 
(figure  C9),  it  is  to  be  expected  for  wavelengths  short  as  compared  to 
receiver  separation.  A  look  at  the  cross-correlation  reveals  the  highest 
correlation  for  a  lag  corresponding  to  an  end- on  wave  approach  from  above 
(figure  CIO),  suggesting  sources  at  the  sea  surface. 

Medium  Frequency  Band:  5-100  Hz 

Sound  from  distant  shipping,  travelling  horizontally  large  distances 
(hundreds  to  thousands  of  kilometers),  often  in  the  sound  channel,  will 
depend  on  both  traffic  density  and  speed,  and  propagation  path 
characteristics.  Time  variations  are  expected  to  be  on  the  order  of 
hours  but  may  be  punctuated  by  local  sources  such  as  ships  overhead 
(figure  Cll) .  It  is  because  these  waves  from  long  distant  sources  travel 
very  close  to  horizontal  that  the  caldera  wall  can  influence  the  sound  at 
the  vent  field.  All  the  sound  arriving  from  the  west  must  diffract 
around  the  edge  of  the  wall  to  reach  the  lower  receiver  as  it  is  in  the 
geometric  shadow  of  eastward  travelling  acoustic  waves  (figure  C3).  In 
addition,  westward  travelling  waves  can  reflect  off  the  wall  and 
increase  the  power  levels  at  the  upper  receiver,  while  the  lower 
receiver,  being  so  close  to  the  bottom,  will  intercept  fewer  reflected 
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waves.  Several  effects  are  seen  in  the  power  spectra  which  support  this 
hypothesis.  First,  the  upper  hydrophone  has  up  to  four  times  as  much 
power  in  the  affected  bandwidth  as  the  lower  (figure  C12a,b,c).  Second, 
the  effect  should  drop  off  at  low  frequencies  as  the  wavelength  of  sound 
approaches  the  dimensions  of  the  wall.  Wavelengths  near  75  m  and  longer, 
20  Hz  and  below,  will  be  less  influenced  by  the  wall  and  hence  reach  both 
receivers  more  equally,  as  can  be  seen  in  figure  C12a,b.  Power 
differences  should  also  be  reduced  at  higher  frequencies,  above  200  Hz, 
because  the  predominant  sound  at  these  frequencies  is  not  travelling 
horizontally,  as  was  seen  above  (figure  C12c) .  Third,  the  second 
deployment  should  show  slightly  greater  differences  since  Hell  vent  is  15 
m  closer  to  the  caldera  wall  than  Inferno;  this  is  seen  in  figure  C12a, 
and  C12b.  Finally,  the  phase  difference  of  the  waves  at  the  two 
receivers  should  be  roughly  zero  and  constant  as  the  sound  hits  the  array 
broadside  from  westward  travelling  waves.  Figure  C13  presents  coherence 
and  phase  for  a  record  taken  during  Hell  deployment.  The  generally  low 

coherence  is  expected  if  one  receiver  samples  part  of  a  sound  field  not 

sampled  by  the  other.  However,  the  coherence  is  markedly  above  random, 
and  the  phase  significant  as  can  be  seen  by  comparing  them  to  coherence 

and  phase  calculated  between  two  totally  uncorrelated  samples  taken  hours 

apart  (figure  C14) .  The  phase  between  20-100  Hz  on  Hell  8  is  quite 
constant  and  near  zero  (compare  to  uncorrelated  phases  in  figure  C14) , 
implying  that  the  sound  impinges  at  right  angles  to  the  vertical  array 
and  is  thus  travelling  horizontally. 

Low  Frequency  Band:  0.25-10  Hz 
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There  appear  to  be  two  major  types  of  sound  in  this  band,  one, 
microseisms  caused  by  local  sea  surface  waves,  and  the  other  of  unknown 
origin  which  produces  evanescent  interface  waves  travelling  along  the 
seafloor  boundary  whose  amplitudes  decay  exponentially  with  distance  from 
the  boundary  (Dowling  and  Ffowcs  Williams,  1983).  The  power  falls  off 
rapidly  with  height,  implying  short  wavelengths  and  therefore  a  local 
source  generating  these  interface  waves.  The  power  levels  in  this  band 
for  the  interface  waves  is  quite  variable  and  intermittent.  On  the  other 
hand,  mictoseism  power  changes  more  slowly  as  it  is  tied  to  changing  sea 
state.  The  upper  hydrophone,  being  less  influenced  by  interface  waves, 
should  have  slower  and  lower  amplitude  variability  than  the  lower 
hydrophone.  Figures  C15-C17  depict  this  effect  and  show  the  striking 
difference  between  microseism  and  interface  waveforms.  A  record  from 
Inferno,  figure  C15,  and  one  from  Hell,  figure  C16,  show  nearly  identical 
low  frequency  waves,  in  phase,  as  expected  from  long-wavelength 
microseisms.  The  second  record  from  Hell  (figure  C17),  taken  four  hours 
earlier,  shows  high  amplitude  interface  waves  on  the  lower  hydrophone  and 
typical  microseisraic  amplitude  waves  on  the  upper,  exemplifying  the 
effect  that  interface  wave  amplitude  dies  away  exponentially  from  the 
seafloor.  The  power,  coherence  and  phase  confirm  the  conclusions  of  the 
visual  examination  of  the  time  series.  In  figures  C18a  and  C18b  can  be 
seen  the  high  coherence  and  constant  phase  for  microseisms.  In  figures 
C19a  and  C19b  are  shown  the  low  coherence  and  variable  phase  for  the 
interface  waves . 

In  summary,  the  general  power  spectral  level,  time  variability, 
power  differences  and  wave  phase  relationships  suggest  that  the  bulk  of 
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recorded  sound  is  attributable  to  the  identifiable  ocean  noise  sources  of 
microseisms,  distant  shipping,  and  ocean  surface  weather.  In  addition, 
several  records  present  striking  evidence  of  evanescent  interface  waves 
generated  by  a  local  source. 

Isolated  Events 

Several  individual  noise  events  were  recorded  during  the  course  of 
this  experiment.  One  record  taken  during  Inferno  deployment  captured  a 
monochromatic  19  Hz  sound  signal  which  we  attribute  to  a  whale  (figure 
C20) .  The  signal  has  high  power  (figure  C21) ,  coherence  (figure  C22a) 
and  the  phase  shows  that  the  signal  is  coming  in  end- on  from  above 
(figure  C22b) .  The  separation,  d,  of  the  receivers  is  such  that  for  this 
frequency,  d-hA,  so  sound  arriving  from  above  will  have  a  phase  delay 
between  the  two  receivers  of  180°,  as  is  seen  in  figure  C22b.  Curiously, 
the  amplitude  is  higher  on  the  lower  hydrophone  than  on  the  upper.  This 
may  be  due  to  the  reflected  wave  from  the  hard  bottom  causing 
constructive  interference  on  the  lower  hydrophone  and  destructive 
interference  on  the  upper  hydrophone  h  wavelength  away. 

Similar  power  and  phase  phenomena  are  seen  on  a  record  taken  when 
the  Atlantis  II  is  the  sound  source.  Figure  C23  shows  the  power  spectra 
of  noise  with  the  ship  overhead  and  slightly  eastward  of  the  vent  site. 
There  is  very  high  coherence  in  the  signals  and  a  linear  phase  with 
slope  describing  waves  coming  in  end-on  (figure  C24) .  At  20  Hz  the  180° 
phase  shift  is  evident,  similar  to  that  of  the  whale. 

Hydrothermal  Noise  Detection 
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Neither  Hell  nor  Inferno  vents  produce  enough  sound  to  isolate 
unambiguously  their  entire  spectrum  levels  from  ambient  ocean  noise.  In 
fact,  as  seen  above,  overall  noise  levels  are  higher  far  from  the  vents 
than  they  are  close  to  them.  This  unexpected  result  makes  direct 
comparisons  of  simultaneous  but  spatially  separated  noise  spectra 
difficult.  In  order  to  identify  regions  of  the  spectra  which  may  have 
vent  signal,  we  must  identify  a  background  noise  level  for  this  site.  It 
had  been  hoped  that  the  far  receiver  would  provide  this  information,  but 
geometrical  influences,  ie  the  caldera  wall,  proved  to  be  too  significant 
and  complex.  As  an  alternative,  we  have  chosen  to  represent  background 
ambient  ocean  noise  with  the  lowest  recorded  noise  spectra.  Anything 
recorded  significantly  above  this,  as  defined  below,  will  be  considered  a 
signal  which  we  will  attempt  to  attribute  to  known  sources. 

Since  all  analyses  are  within  a  statistical  framework,  it  is  useful 
to  assign  significance  by  examining  the  probability  that  a  giver,  signal 
is  real  and  not  a  chance  fluctuation.  The  effective  signal  to  noise 
level  will  be  described  here  using  a  value  called  the  detection  index 
(Burdick,  1984;  Bangs  and  Schultheiss,  1973;  Owsley  and  Swope,  1981). 

This  is  an  array  output  signal  to  noise  ratio  which  depends  on  the  number 
of  receivers,  the  sample  length,  and  the  signal  field  and  noise  field. 

It  describes  both  the  probability  of  detection  and  the  probability  of 
false  alarm  (mistakenly  detecting  a  signal  when  its  not  present)  for  a 
given  signal  and  noise  input  to  the  hydrophones  (see  appendix).  A  high 
detection  index  provides  simultaneously  high  confidence  in  correctly 
detecting  a  signal  when  it  is  present  and  in  not  falsely  detecting  one 
when  it  is  absent.  We  will  calculate  this  detection  index  by  assigning 
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the  lowest  power  level  record  to  be  noise,  comparing  the  other  records  to 
these  values  and  designating  as  signal  anything  above  the  noise.  The 
magnitude  of  the  resulting  detection  indices  will  represent  our 
confidence  in  a  given  signal.  We  will  deem  as  significant  any  detection 
index  over  10  -  this  value  represents  a  probability  of  detection  of  90% 
and  probability  of  false  alarm  of  5%. 

The  average  of  the  lowest  noise  spectra,  one  set  from  Hell  and  the 
other  from  Inferno,  assigned  to  represent  ambient  ocean  noise  (figure 
C25),  are  remarkably  similar  considering  they  were  recorded  several  days 
apart.  The  lower  hydrophone  near  Hell  will  be  used  for  comparison  to 
Inferno  and  vice  versa,  based  on  the  logic  that  the  vents  are  not  exactly 
identical  and  any  differences  in  sound  level  due  to  venting  should  be 
apparent  given  the  35  m  separation  between  vents. 

When  ws  compare  the  quietest  record  at  Hell  with  background  Inferno 

noise  we  find  that  the  detection  indices  at  all  frequencies  are  less  than 

ten,  as  shown  in  figure  C26.  Also,  when  the  quietest  record  from  Inferno 

is  compared  to  background  at  Hell,  the  indices  are  all  less  than  ten 

(figure  C26.  We  must  conclude  that  the  sound  produced  by  the  vents  are 

below  the  90%  confidence  detection  limit  of  this  array  in  this  noise 

field,  a  detection  limit  that  is  equivalent  to  a  SNR  of  about  10  in  a  5 

-4  2 

Hz  band,  or  about  10  Pa  /Hz. 

We  will  note  hei.ep  however,  that  Inferno  appears  to  be  generating  a 
small  sound  signal  at  about  40  Hz,  as  seen  in  figure  C25,  of 
approximately  1x10'  Pa^/Hz.  This  signal  is  evident  on  the  two  lowest 
noise  300  Hz  bandwidth  records.  The  other  five  contain  high  noise  on 
both  channels  at  this  frequency  as  well  as  others,  probably  due  to  an 
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undocumented  ship  or  submarine  in  the  area,  and  no  conclusions  can  be 


drawn  from  these.  The  frequency  and  power  level  of  Inferno  vent  are 
consistent  with  theoretical  estimates  of  noise  from  a  jet  (Little  and 
Stolzenbach,  in  prep.).  A  velocity  of  2  ra/s  from  the  0.04  n  orifice 
would  generate  sound  with  a  peak  frequency  of  40  Hz.  Sound  pressure 
level  at  2.3  m  in  the  near  field  region  of  a  cold  turbulent  jet  is 
predicted  to  be  2x10  ^Pa  at  peak  frequency,  (a  power  level  of  about  5x10" 

8  2  o 

Pa  /Hz).  A  hot,  350  C  jet  exiting  into  cold  seawater  can  produce  up  to 

a  factor  of  10^  higher  power  levels,  through  the  amplification  of  sound 

by  convecting  flow  inhomogeneities  (Morfey,  1973).  The  data  from  the 

3 

Inferno  deployment  show  a  factor  of  about  10  elevation  in  power  levels 
over  expected  cold,  jet  sound,  which  we  attribute  to  the  high  temperature 
and  consequently  low  density  of  the  exiting  hydrothermal  fluid. 

Comparison  with  other  Data 

Noise  measurements  made  in  Axial  Caldera  two  years  earlier  (Bibee 
and  Jacobson,  1986)  are  comparable  to  the  data  from  this  experiment.  The 
1985  instrument  placed  closest  to  the  vents  detected  sound  about  a  factor 
of  10  in  power  lower  than  our  lowest  noise  record.  Anomalously  high 
noise  from  the  northern  part  of  the  caldera  measured  on  the  1985 
experiment  remains  unexplained  but,  in  light  of  our  experiment,  it  is  not 
from  a  typical  black  smoker  vent  field  like  Ashes. 


CONCLUSIONS 

High  quality  data  collected  on  this  experiment  were  used  to 
characterize  the  noise  field  associated  with  Ashes  Vent  Field  on  Axial 
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Seamount.  Narrow  band  temporal  variability  was  dominated  by  local  ship 
traffic,  submarines  and  whales.  Broadband  variability  is  attributed  to 
changes  in  distant  shipping  or  propagation  paths  thereof,  and  to  changes 
in  local  weather. 

The  very  low  frequency  signals  (0.5-2Hz)  appear  to  be  from  two 
distinct  sources,  one  continuous  microseisms  and  the  other  intermittent 
local  events.  Microseisms  which  are  equal  in  amplitude  on  both  receivers 
appear  in  almost  all  records.  Sporadically  appearing  interface  waves 
show  high  pressure  amplitudes  near  the  seafloor  and,  decaying  rapidly 
with  vertical  distance,  produce  very  low  signal  levels  at  40  m  above  the 
bottom. 

No  continuous  vent  signals  were  deemed  significant  based  on  a 

criterion  of  90%  probability  of  detection  and  5%  probability  of  false 

alarm.  However,  a  small  signal  near  40  Hz,  with  a  power  level  of  1x10 
2 

Pa  /Hz  was  noticed  on  two  records  from  the  Inferno  deployment.  This 
frequency  of  this  signal  is  consistent  with  predictions  and  the  power 
level  suggests  the  occurrence  of  jet  noise  amplification  due  to  convected 
density  inhcmogeneities . 

Detection  of  vent  signals  was  hampered  by  the  unexpected  spatial 
non-stationarity  where  receivers  38  m  apart  had  about  a  factor  4  ambient 
noise  level  difference  due  to  effects  of  local  topography.  This  made 
absolute  power  level  comparisons  impossible  and  the  resulting  low  signal 
coherences  reduced  confidence  in  phase  and  directionality  results. 

Future  experiments  will  need  to  include  more  receivers  to  improve 
the  effective  SNR.  Several  hydrophones  should  be  located  close  enough  to 
the  vent  (-5  m)  to  allow  beamforming  through  coherence  and  phase 
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calculations.  In  addition,  to  improve  SNR,  a  more  active  vent  site 
should  be  chosen  since  acoustic  power  output  increases  rapidly  with 
increasing  fluid  exit  velocity. 

To  generalize  from  one  vent  site,  it  appears  that  the  method  of 
monitoring  hydrothermal  vent  fluxes  through  passive  acoustics  will 
require  more  than  a  simple  surface  ship  deployment  of  a  two -element 
hydrophone  array.  It  may  be  most  useful  on  high  output  vent  sites  rather 
than  on  low  to  moderate  hydrothermal  areas . 
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APPENDIX 


The  signal  detection  ability  of  a  hydrophone  array  in  a  noise  field 
is  calculated  by  determining  the  raw  input  signal  level  required  to 
achieve  a  given  probability  of  detection  while  maintaining  an  acceptable 
probability  of  false  alarm.  The  probability  of  detecting  a  signal  within 
a  gaussian  noise  distribution  depends  on  the  amplitude  of  the  signal, 
amplitude  and  coherence  length  of  noise,  number  of  observations,  and 
number  of  receivers.  A  matrix  description,  P^,  of  the  signal  waveform 
for  a  receiver  array  with  a  point  source  located  at  the  hydrothermal 
chimney  orifice,  can  be  generated  from  a  description  of  the  waveform  at 
each  receiver  (from  Owsley  and  Swope,  1981),  V^: 

V^.  —  [  a^e  ^  t  ,  a2e  ^  ^  }  da) 

Here  f  is  frequency,  j  is  J(- 1) ,  rl  and  r2  are  receiver -source  separation 
distances  (figure  C27) ,  c  is  sound  speed  and  a  is  an  amplitude  factor 
accounting  for  spherical  spreading  loss: 

(2a) 


and, 

P  -  V  V  + 
f  f  f 


(3a) 


The  correlation  length  of  the  ambient  noise  field  in  the  ocean  will 
determine  how  difficult  noise  removal  will  be  given  a  fixed  receiver 
separation.  We  will  assume  that  the  noise  field  is  isotropic  both  for 
simplicity  and  since  the  results  are  not  significantly  different  from 
either  surface  generated  or  azimuthally  distributed  noise  fields.  For 
isotropic  noise,  the  correlation  function  q  is  (Burdick,  1984): 


q  —  sin(2?rd/A)/(2?rd/A) 


(4a) 
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where  d  is  receiver  separation  and  A  is  wavelength.  The  array  noise 
correlation  function  for  two  receivers  will  then  be: 


Q£  •  [;  ?] +  •  [i  ?] 

where  e  is  zero  for  perfectly  bandlimited,  perfectly  predictable,  noise 
but  for  real,  ocean  systems  with  a  random  uncorrelated  noise  component,  « 
must  be  non-zero.  For  isotropic  noise,  £  is  approximated  by  0.1  (J. 
Krolik,  personal  communication) . 

The  effective  signal  to  noise  ratio  output  (SNRq)  of  a  system  in  a 
coherent  noise  field  is  (from  Bangs  and  Schultheiss,  1973): 

SNRq  -  (Sf/Nf ) 2  trace (  (?fQf " X) 2 ]  )  (6a) 

where  N  is  number  of  sample  sets  and  S^/N^  is  the  SNR  as  a  function  of 
frequency  input  to  each  receiver.  Equation  6a  is  used  to  calculate  the 
output  SNR^  given  an  input  S^/N^.  For  the  100Hz  bandwidth,  S^/N^.  is 
calculated  using  an  average  of  the  quietest  records  from  the  lower 
hydrophone  of  the  opposite  deployment  as  the  value  of  ambient  noise  N^. 
This  must  be  done  because  the  upper  hydrophone  of  both  deployments  had 
such  high  noise  levels.  The  assumption  inherent  in  this  procedure  is 
that  the  two  vents  are  not  producing  exactly  the  same  signal.  The  SNRq 
is  then  calculated  in  5.3Hz  frequency  bands.  The  SNRq  equals  Burdick's 
(1984)  detectability  index  d  .  We  have  chosen  as  significant  a  dQ 
greater  than  10,  which  is  equivalent  to  a  probability  of  detection  of  90% 
and  a  probability  of  false  alarm  of  5%  (Burdick,  1984). 
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Figure  C2 

Location  of  Ashes  Vent  Field  In  Axial  Seamount,  on  the  Juan  de  Fuca 
Ridge.  The  field  Is  just  east  of  the  50-100  m  high  scarp  defining 
the  caldera  wall  (Embley  et  al.,  1988).  An  earlier  seismic 
experiment  had  a  hydrophone  near  CASH  vent  field  In  the  northern 
part  of  the  caldera. 


Plan  View 


- — West  Distance  m  East — ► 


Side  View 


150  “MO  130  120  110  100  90  80  70  60  50  40  30  20  10  0 

♦—West  Distance  m  East- 


Figure  C3 

a)  Plan  view  of  Ashes  Vent  Field  showing  black  smokers  Hell  and 
Inferno,  and  white  smokers  Mushroom,  Hillock,  and  Virgin  Mound, 
with  caldera  wall  60-80  m  to  the  west  of  deployment  sites. 

b)  Side  view  of  deployment  showing  relationship  between  wall,  vents 
and  hydrophone  array.  Eastward  travelling  waves  are  prevented  from 
reaching  lower  hydrophones  by  the  shadowing  effect  of  this  wall. 
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Figure  C4 

Detail  of  instrumented  cable  and  deployment  scheme.  Lower  hydrophone 
was  2-3  m  from  vent  orifice,  upper  hydrophone  was  38-40  m  above  lower. 


101 


Power  (PaA2/H.-. 


Inferno  10  Low  Noise 


Figure  C5 

Low  noise  period  power  spectral  density  from  Inferno  deployment  (record 
10)  of  bandwidth  l-800Hz.  Solid  line  is  for  lower  hydrophone  (close  to 
the  vent)  and  dotted  line  is  for  upper  hydrophone.  Differences  in  the 
bandwidth  30-250  Hz  are  attributed  to  effects  of  the  caldera  wall  on 
ambient  ocean  noise.  The  mathematical  spectrum  estimation  for  all 
spectra  shown  in  this  paper  is  accurate  to  +0.13  x  power  level.  In 
addition,  the  spectra  shown  Include  the  calibration  between  upper  and 
lower  hydrophones,  a  difference  which  was  less  than  a  factor  of  1.6  for 
the  frequencies  presented  here. 
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Figure  C6 

Low  noise  period  power  spectral  density  from  Hell  deployment  (record  4) 
of  bandwidth  l-100Hz.  Solid  line  is  for  lower  hydrophone  (close  to 
vent)  and  dotted  line  is  for  upper  hydrophone.  This  details 
differences  in  sound  level  caused  by  the  wall  in  the  bandwidth  0.5- 
100  Hz. 
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Figure  C7 

Hourly  time  history  of  spectra  from  Inferno  deployment  (lower 
hydrophone)  of  bandwidth  100-800HZ  showing  Increasing  tnergy  toward 
end  of  deployment  which  corresponds  to  increasing  wind  velocity. 
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Figure  C8 

Low  noise  period  power  spectral  density  from  Inferno  deployment  (record 
10)  showing  difference  in  upper  (dotted  line)  and  lower  (solid  line) 
hydrophone  power  levels  out  to  250  Hz,  followed  by  similar  levels. 

This  due  to  the  caldera  wall  affecting  horizontally  travelling  20-250 
Hz  waves  and  not  influencing  250-800  Hz  vertically  travelling  waves. 
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Figure  CIO 

Cross  correlation  between  upper  and  lower  hydrophone  for  high  pass 
filtered  (above  200  Hz)  Inferno  low  noise  period  (record  10).  Greatest 
magnitude  of  correlation  occurs  at  a  lag  of  64  points,  which 
corresponds  to  the  receiver  separation,  implying  that  the  waves  are 
impinging  on  the  array  from  end-on  (from  above). 
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Figure  Cll 

a)  Hourly  tlae  history  of  spectra  from  Inferno  deployment  for  bandwidth 
5-100  Hr,  shoving  a  general  Increase  during  the  middle  of  the 
deployment,  probably  In  response  to  an  Increase  In  ship  activity, 

b)  Hourly  time  history  of  spectra  from  Hell  deployment  for  bandwidth 
5-100  Hr,  shoving  relatively  constant  levels  except  for  the  first  and 
last  hour  during  which  time  the  Atlantis  II  was  in  close  proximity 

to  the  vent  elta.  inQ 
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Figure  C12 

Ratio*  of  upper  and  lower  hydrophone  power  levels  for  all  records  from 
a)  Hell  deployment,  b)  Inferno  deployment:  1-10  Hz  and,  c)  Inferno 
deployment:  1:800  Hz.  Notice  that  the  quotient  decreases  with  both 
decreasing  and  Increasing  frequency,  centered  around  typical  distant 
shipping  frequencies  of  30-200  Hz.  This  difference  Is  due  to  effect* 
of  the  caluera  wall. 
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Figure  C13 

a)  Power  spectral  density  for  low  noise  period  of  Kell  deployment 
(record  8).  Solid  line  is  lower  hydrophone,  dotted  line  Is  upper 
hydrophone . 

b)  Coherence  between  upper  and  lower  hydrophones  for  low  noise  period 
of  Hell  deployment  (record  8).  Coherence,  although  low.  Is  significant 
(compare  to  figure  14) . 

c)  Phase  for  same  period  showing  relatively  constant  and  zero  phase 
produced  by  waves  Impinging  on  the  array  broadside  because  they  were 
travelling  horizontally. 
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Hell  uncomelated 


Figure  C14 

«)  Coherence  calculated  between  records  taken  fron  upper  and  lover 
hydrophones  at  two  different  tines  (totally  uncorrelated).  Mote  the 
lev  level  of  thla  randoa  coherence. 

b)  Phase  for  sane  records.  Note  the  coopletely  random  phase 
(especially  as  compared  to  figure  13b). 
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Figure  C15 

Time  series  pressure  level  from  Inferno  deployment  (record  1),  upper 
trace  Is  from  upper  hydrophone  and  vice  versa.  Note  similarity  in 
phase  and  pressure  level  for  low  frequency  oscillations  on  both 
hydrophones  due  to  microselsms. 
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Figure  C16 

Time  series  pressure  level  from  Hell  deployment  (record  9),  upper 
trace  Is  from  upper  hydrophone  and  vice  versa.  Note  similarity  In 
phase  and  pressure  level  for  low  frequency  oscillations  on  both 
hydrophones  due  to  atcroselsas. 
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Figure  Cl 7 

Time  series  pressure  level  from  Hell  deployment  (record  5),  upper 
trace  is  from  upper  hydrophone  and  vice  versa.  Note  large  amplitude 
erratic  oscillations  on  lower  hydrophone  not  seen  on  upper  hydrophone. 
Upper  hydrophone  shows  typical  microseisra  oscillations  while  lower 
records  passage  of  Interface  waves  which  poorly  penetrate  water  column. 
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Hell  9  0.5-10  Hz 


Flgur*  Cl 8 

a)  Upper  part  ahowa  power  levels  froa  Hell  (record  9)  In  0.5-10  Hz 
bandwidth.  Note  similarity  In  power  between  upper  hydrophone  (dotted 
line)  and  lower  hydrophone  (solid  line)  due  to  nlcroselsa  energy. 

b)  Lower  part  shows  high  coherence  between  the  two  hydrophones  from 
0.5-3  Ha  characteristic  of  nlcroselsns. 

c)  This  shows  the  constant  and  zero  phase  between  upper  and  lower 

hydrophone .  ik 


10" 


Hell  5  0  3- 10  Hi 


Hell  3 


Figure  C19 

a)  Upper  pare  shows  high  power  level  on  lower  hydrophone  (solid  line) 
and  lower  power  on  upper  hydrophone  (dotted  line)  for  Hell  deployment 
(record  5). 

a)  Low  coherence  between  the  two  hydrophones  is  seen  since  pressure 
levels  at  lower  hydrophone  are  due  to  Interface  waves  not  seen  on 
upper  hydrophone . 

b)  This  shows  non-constant  and  non-zero  phase  as  is  expected  between 

the  upper  and  lower  hydrophones  since  they  are  receiving  different 
signals.  116 
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Figure  C20 

Upper  trace  is  pressure  level  from  upper  hydrophone,  lower  trace  Is 
pressure  level  at  lower  hydrophone.  The  beginning  of  this  period 
from  Inferno  (record  3)  a  monochromatic,  19  Hz  signal  is  detected, 
which  we  attribute  to  a  whale. 


Power  (PaA2/Hz) 


Figure  C21 

Power  spectral  density  for  Inferno  (record  3),  solid  line  is  from 
lower  hydrophone,  dotted  line  is  from  upper  hydrophone,  showing  peak  at 
19  Hz  attributed  to  a  whale.  The  amplitude  is  higher  on  lower 
hydrophone  due  to  reflections  off  the  hard  bottom  causing  constructive 
interference  on  the  lower  phone  and  destructive  on  the  upper. 
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Figure  C23 

Power  spectral  density  from  Hell  deployment  (record  12)  when  the 
Atlantis  II  is  overhead  and  slightly  to  the  north  east  of  the  vent 
site.  Peaks  and  harmonics  are  typical  of  near  ship  noise. 
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Figure  C24 

a)  Coherence  between  two  receivers  for  Hell  (12)  when  ship  Is  near. 
Notice  extremely  high  coherence  at  ship's  frequencies. 

b)  Transfer  phase  between  the  two  receivers  for  this  record  at 
high  power  and  coherence  freque'  ties.  Linear  trend  shown  by  dotted 
line  would  be  produced  by  waves  arriving  from  nearly  vertical.  Phase 
from  ship  noise  falls  mainly  on  this  line,  or  +  180°  from  It  due  to 
reflections . 
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Figure  C25 

Low  noise  records  used  for  calculation  of  signal  detection  index, 
average  of  Inferno  records  1  and  11,  (upper  solid  line  from  upper 
hydrophone,  lower  solid  line  from  lower  hydrophone)  and  Hell  records 
4,  8,  and  9  (upper  dotted  line  from  upper  hydrophone,  lower  dotted 
line  from  lower  hydrophone) .  Note  small  peak  in  lower  Inferno  power 
spectral  density  near  38  Hz;  this  is  attributed  to  jet  noise  from 
Inferno  vent. 
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Figure  C26 

«)  Signal  detection  index  for  low  noise  Hell  deployment  (record  11) 
note  that  nowhere  le  the  SNR  greater  than  10,  the  value  which 
correspond*  to  a  90%  confidence  in  *lgnal  detection  with  a  5*  false 
alarm  rat*. 

b)  Detection  index  for  low  noise  Inferno  deployment  (record  1) ,  not* 
that  nowhere  is  the  SNR  greater  than  10.  the  value  which  corresponds 
to  a  90%  confidence  in  signal  detection  with  a  3%  false  alarm  rat*. 


OCEAN  NOISE  MEASUREMENTS  NEAR  THE  TAG  HYDROTHERMAL  AREA, 
MID-ATLANTIC  RIDGE,  LATITUDE  26°N 

ABSTRACT 

In  light  of  the  potential  utility  of  a  remote,  passive  acoustic 
technique  for  monitoring  variations  in  hydrothermal  vent  fluxes,  we  have 
examined  ambient  noise  in  the  frequency  band  1-30  Hz  at  a  range  of  0.75- 
14  km  from  the  TAG  (Trans -Atlantic  Geotraverse)  hydrothermal  area  on  the 
Mid-Atlantic  Ridge  near  26°N,  the  site  of  an  extremely  active  high 
temperature  hydrothermal  vent  field  (Rona,  1986).  The  ambient  noise 
field  exhibits  great  temporal  and  spatial  variations  attributed  in  part 
to  typical  ocean  noise  sources  such  as  distant  shipping  and  microseisms. 
Power  spectral  levels  as  measured  at  each  of  six  ocean  bottom  hydrophones 
(OBH)  were  used  to  estimate  the  location  of  point  sources  of  sound  in  the 
area,  if  any. 

The  hydrothermal  vent  did  not  produce  enough  sound  to  be  located  as  a 
point  source  using  data  from  the  OBH  array.  The  only  consistently 
identifiable  point  source  found  with  the  data  set  was  generating  sound  in 
a  0.8-3. 5  Hz  bandwidth  and  located  outside  the  median  valley.  It  appears 
to  be  harmonic  tremor  associated  with  the  tip  of  a  ridge  on  the  western 
side  of  the  spreading  axis  and  may  be  volcanic  in  origin. 
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INTRODUCTION 


High  temperature  hydrothermal  fluid  travels  through  fractured 
oceanic  crust  in  narrow  upflow  zones  to  exit  at  the  seafloor  in  flows 
ranging  from  high  velocity  buoyant  jets  emerging  from  tall  chimneys  to 
diffuse  seepage  from  local  cracks  and  fissures.  This  diversity  of  flow 
makes  detailed  characterization  of  vent  fields  over  space  and  time 
extremely  difficult.  Without  this  information,  precision  modelling  of 
hydrothermal  systems  from  chemical,  geological,  fluid  dynamical  and 
biological  standpoints  is  nearly  impossible.  The  lack  of  long-term  flow 
information  motivated  a  study  of  the  use  of  passive  acoustics  for 
monitoring  flow  characteristics  (Little  and  Stolzenbach,  in  prep).  The 
sound  field  produced  by  a  turbulent  jet  is  theoretically  heavily 
dependent  on  fluid  velocity  and  exit  geometry  (Lighthill,  1952).  Such  a 
sound  field  could  be  used  to  monitor  changes  in  these  jet  parameters  over 
time.  A  technique  to  monitor  and  understand  changes  in  sound  output  as 
they  relate  to  fluid  flow  would  give  researchers  a  much  needed  tool  for 
studying  long-term  temporal  variations  in  hydrothermal  activity.  The 
feasibility  of  this  method  will  depend  on  an  ability  to  detect  and 
identify  hydrotherraally  generated  sound  in  a  background  of  ambient  ocean 
noise.  Recently  published  papers  have  reported  anomalous  sound  fields, 
in  the  bandwidth  1-30  Hz,  associated  with  hydrothermal  centers,  one  at 
21°N  on  the  East  Pacific  Rise  (EPR)  (Reidesel  et  al.,  1982)  and  one  on 
the  Juan  de  Fuca  Ridge  (Bibee  and  Jacobson,  1986).  The  mechanism 
responsible  for  sound  generation  in  these  two  cases  has  not  been 
identified.  In  light  of  these  results  and  the  potential  utility  of  a 
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passive  acoustic  monitoring  technique  for  hydrothermal  vents,  we  have 
examined  ambient  noise  collected  near  the  TAG  (Trans -Atlantic 
Geotraverse)  hydrothermal  area  on  the  Mid- Atlantic  Ridge  near  26°N,  the 
site  of  an  extremely  active  high  temperature  hydrothermal  vent  field 
(Rona,  1986).  A  microearthquake  study  of  the  TAG  area  (Kong  et  al, 

1986) ,  conducted  prior  to  the  discovery  of  the  high  temperature  vent 
field,  conveniently  provided  a  data  set  with  which  to  examine  the  ambient 
sound  field  in  the  frequency  band  1-30  Hz  at  a  range  of  0.75  to  14  km 
from  the  vent  site.  The  purpose  of  examining  this  data  set  was  first,  to 
determine  whether  the  hydrothermal  vent  could  be  identified  as  a  source 
of  sound,  and  second,  whether  any  other  point  sources  of  sound  could  be 
located  within  the  study  area. 

AMBIENT  OCEAN  NOISE 

The  level  of  ambient  ocean  noise  in  a  given  frequency  band 
determines  the  degree  of  difficulty  associated  with  identifying  and 
detecting  an  unknown  sound  source.  Noise  in  the  ocean  is  typically  non¬ 
stationary  both  in  time  and  space  (Burdick,  1984;  Hodgkiss  and  Anderson, 
1980).  The  noise  field  varies  considerably  over  time  due  to  the  rich 
variety  of  sources  (Akal  et  al.,  1986)  and  measurements  made  hours  apart 
can  show  striking  dissimilarity.  Spatially,  the  noise  field  power  levels 
can  be  orders  of  magnitude  different  depending  on  topography,  water  sound 
velocity  structure,  and  proximity  of  receiver  to  the  seafloor  or  sea 
surface  (Urick,  1986;  Burdick,  1984;  Wentz,  1962).  On  average,  typical, 
ubiquitous,  deep-water  ambient  ocean  noise  spectra  in  the  frequency  range 
covered  by  this  experiment  are  caused  by  two  different  source  mechanisms. 
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The  lowest  frequency  band,  0.01-5  Hz,  figure  Dl,  exhibits  high  power 
levels  and  is  dominated  by  microseisms:  low  frequency  pressure 
disturbances  caused  by  non-linear  interactions  of  ocean  surface  waves 
(Webb,  1984),  and  teleseismic  events.  Temporal  variations  are  as  rapid 
as  a  few  minutes  for  interface  waves  travelling  near  the  sea  floor  from 
teleseismic  and  local  sources,  and  as  long  as  a  few  hours  for  variations 
in  sea  state  and  swell. 

Power  in  the  higher  frequency  band  from  5-100  Hz,  produced  by 
distant  shipping,  is  dependent  on  sound  that  has  travelled  tens  to 
hundreds  of  kilometers  and  is  strongly  influenced  by  wave  guide 
propagation  effects  such  as  sound  channelling  due  to  a  velocity  minimum 
in  the  water  column.  Travelling  great  distances  through  the  ocean 
broadens  the  spectral  peak  and  removes  waves  travelling  outside  a  few 
degrees  of  horizontal.  Receivers  placed  in  a  depression  or  valley  are 
likely  to  be  screened  off  somewhat  from  such  distant  sources  and  have 
lower  noise  levels  (Urick,  1986).  Variations  in  this  band  are  slow,  on 
the  order  of  hours  to  days  and  depend  on  changes  in  ship  traffic  and 
large  scale  temperature  and  salinity  structure  (in  the  sound  channel). 

Transient  sources  also  contribute  to  noise  at  a  given  location. 

These  include  local  ships,  which  are  characterized  by  high  energy  narrow 
band  peaks,  often  including  harmonics  of  a  fundamental  frequency, 
anywhere  from  5  to  200  Hz.  Biological  sources  can  also  contribute  to  the 
noise  field,  with  whales  and  dolphins  capable  of  producing  high 
amplitude,  short  duration  sounds  at  frequencies  from  18-100  Hz  (Watkins, 
1981) . 
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HYDROTHERMAL  SOUND  GENERATION 


No  unambiguous  field  data  exists  characterizing  vent  sound,  but 
theoretical  estimates  on  power  and  frequency,  as  well  as  laboratory 
results  (Little  and  Stolzenbach,  in  prep)  give  a  range  of  values  to  help 
delimit  the  frequency  bands  of  possible  energy.  Sound  produced  by 
turbulent  mixing,  as  is  found  at  the  jet  orifice,  will  be  generated  at  a 
frequency  close  to  V/D  where  V  is  exit  velocity  and  D  is  orifice  diameter 
(Morfey,  1973).  For  typical  vent  velocities  and  geometries,  this 
frequency  range  is  10-100Hz.  Other  source  mechanisms,  such  as  fluid 
pumping  or  cavity  resonance,  may  generate  lower  frequencies,  depending  on 
the  distribution  and  geometry  of  subsurface  fluid  pathways  as  well  as 
hydrodynamic  flow  characteristics  deep  within  the  hydrothermal  system. 
Higher  frequencies  can  be  produced  by  pipe  resonance,  and  are 
proportional  to  c/L  where  c  is  sound  speed  and  L  pipe  length.  In  this 
paper  we  are  confined  to  examining  frequencies  of  1-30  Hz  by  the 
bandwidth  set  for  the  raicroearthquake  study. 


HYDROTHERMAL  SITE  DESCRIPTION 

The  TAG  hydrothermal  area  (Rona  et  al.,  1986;  Thompson,  et  al., 

1988)  is  located  at  26°08.28'N,  44°49.63'W,  on  a  wall  1.5  kilometers  east 
of  the  spreading  axis  of  the  Mid-Atlantic  Ridge  (figure  D2).  This  east 
wall  rises  from  the  valley  floor  at  4000  m  to  a  height  of  2000m  through  a 
series  of  block  faults.  The  hydrothermal  site  is  composed  of  an  inner 
mound  of  sulfides  250m  wide  and  55m  tall  which  sits  atop  an  outer  mound 
580  ra  wide  and  25m  high.  The  base  of  this  structure  is  situated  at  a 
maximum  depth  of  3700  m.  White  smokers  and  shimmering  water,  with 
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temperatures  up  to  226°C,  were  found  around  the  periphery  of  the  inner 
mound.  Tall  chimneys,  at  least  10  m  high,  and  black  smokers  were 
predominantly  located  at  the  very  top  of  the  inner  mound,  emitting  fluid 
at  temperatures  of  350°C  at  velocities  of  about  1-2  m/s.  A  zone  of  low 
temperature  hydrothermal  activity  is  located  between  2400  and  3100  m 
depth  on  the  east  wall,  3.7  km  upslope  to  the  east  of  the  main  field 
(Rona,  1980).  Water  column  anomalies  suggest  a  heat  flux  from  this  10  km 
section  of  ridge  of  5x10®  W  (Klinkhammer ,  et  al.,  1986). 


DATA  DESCRIPTION 

During  the  summer  of  1985,  a  successful  three  week  deployment  of  a 
network  of  six  continuously  recording  ocean  bottom  hydrophones  (OBH) , 
(Koelsch  and  Purdy,  1979;  Kong  et  al. ,  1986)  was  carried  out  within  an 
area  of  19x10  km  on  the  crest  of  the  Mid-Atlantic  Ridge.  This  network 
provided  coverage  of  the  known  low- temperature  hydrothermal  field  on  the 
eastern  wall  (Rona,  1976;  Kong  et  al.,  1986)  and,  coincidentally, 
surrounded  the  subsequently  discovered  high  temperature  vent  field 
(figure  D2) .  The  instruments  were  deployed  within  the  median  valley  (OBH 
1,2, 3,6)  and  on  the  valley  walls  (OBH  5,8).  The  six  instruments  were 
located  around  the  vent  at  ranges  of  approximately  0.75,  3.5,  5.5,  6.5, 
12.0,  and  14.0  km.  Data  for  this  paper  were  obtained  from  a  small  subset 
of  the  earthquakes  and  shots  which  were  recorded  on  all  six  instruments. 
The  data  set  included  thirteen  separate  time  periods  consisting  of 
approximately  8-13  seconds  of  digitized  (at  100  Hz)  ambient  noise 
preceding  each  event.  These  thirteen  sets  of  simultaneously  recorded 
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ambient  noise,  each  12.8  s  long,  spanned  a  period  of  15  days.  The 
reliable  bandwidth  of  the  system  was  1-30  Hz.  Power  spectra  were 
calculated  using  the  Welsh  method,  dividing  each  record  into  256  point 
length  sets,  multiplying  by  a  Hanning  window  and  averaging  the  results. 

Temporal  variations  are  evident  in  the  raw  time  series  data  as  seen 
in  figure  D3  which  depicts  two  records  taken  11  days  apart  by  the  same 
instrument.  The  second  record  shows  high  amplitude  low  frequency 
components  not  seen  in  the  first.  An  examination  in  the  frequency  domain 
reveals  power  variations  over  time  of  three  orders  of  magnitude  on  the 
instrument  nearest  the  vent  (figure  D4a)  and  almost  four  orders  of 
magnitude  at  the  instrument  12  km  away  (figure  D4b) .  Spectra  from  the 
lower  noise  periods  resemble  typical  ambient  ocean  noise  with  high  power 
at  low  frequencies  (0.8-5  Hz),  a  minimum  near  15  Hz,  and  increasing 
levels  out  to  30  Hz  (the  bandwidth  of  this  system)  as  the  effects  of 
distant  shipping  become  important.  Spectral  peaks  at  22  Hz  were  seen 
when  the  R.V.  Knorr  was  within  the  array  and  not  at  other  times. 

Spatial  variability  is  exemplified  by  two  time  series  taken  at  0.75 
and  6.5  km  from  the  vent  (figure  D5) .  A  look  at  the  power  spectra 
reveals  several  orders  of  magnitude  difference  over  the  array  for  a 
number  of  the  events,  figure  D6a,  D6b,  D7a,  D7b.  Plotting  power  versus 
proximity  to  vent  does  not  reveal  any  frequency  at  which  power  decreases 
monotonically  with  distance  for  typical  events  (figures  D8a,  D8b,  D9a, 
D9b)  . 

SOURCE  LOCATIONS 
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In  order  to  determine  whether  the  hydrothermal  vent  site  is 
generating  detectable  sound,  and  if  there  are  any  other  point  sources  of 
sound  near  the  array,  we  use  the  power  spectral  measured  at  each 
instrument.  Our  use  of  power  levels  to  locate  point  sources  is  based  on 
the  assumption  that  such  sources  radiate  acoustic  waves  whose  power 
decays  inversely  with  the  square  of  distance.  To  minimize  contamination 
from  different  sources  we  have  divided  the  frequency  into  3.2  Hz 
bandwidths.  Of  these,  only  ones  with  relatively  high  power  levels  are 
chosen  (0.8-3. 5,  3. 5-6. 6,  6. 6-9. 7,  12.9-16.0,  19.1-22.3,  and  28.5-31.6 
Hz)  in  order  to  reduce  errors  caused  by  the  superposition  of  the 
isotropic  noise  field  in  this  area  and  any  point  source  noise  fields. 
Power  in  each  bandwidth  is  calculated  for  each  of  the  six  instruments. 
These  levels  are  then  compared  to  power  that  would  be  expected  if  a  point 
source  of  amplitude  A,  call  this  a  predictive  source,  was  located  at 
point  x,y  in  a  40  km  x  40  km  grid  surrounding  the  study  area  where: 

P”  — 2~A~2~ 

OT+y) 

To  limit  the  number  of  free  parameters  to  A,  x  and  y,  and  improve  the 

estimate  of  location,  source  depth  was  fixed  at  3750,  the  depth  of  the 

high  temperature  vent  field.  Next,  the  root-mean- square  difference 

between  actual  and  predicted  power  is  calculated  using  the  six 

o 

instruments.  This  value,  the  misfit  in  units  of  Pa  /Hz  for  at  3.1  Hz 
bandwidth,  is  then  plotted  (figure  DlOa  and  DlOb)  for  each  grid  point 
(spacing  x:  2  km,  y:  2  km).  The  minimum  misfit  is  then  the  best  estimate 
for  point  source  location  for  the  given  measured  power  levels  at  the 
specified  bandwidth  and  event.  Perfect  data  from  a  point  source  would 
yield  a  misfit  of  0  at  the  source  location.  The  misfit  value  approaches 
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the  rms  value  of  the  six  instrument  power  levels  as  the  predictive  source 

is  placed  further  and  further  away  from  the  array.  High  misfits  occur 

near  the  instruments  (provided  the  actual  source  is  not  there)  because 

the  calculated  power  becomes  very  high  for  predictive  sources  placed  near 

the  receivers.  Sources  were  sought  for  each  of  the  six  bandwidths  and  13 

events.  Some  events  and  bandwidths  do  not  result  in  isolatable  misfit 

minima  and  hence  have  not  point  source  solutions.  This  would  be  the  case 

for  sound  recorded  from  an  isotropic  noise  field.  Figure  DlOa  shows  a 

contour  plot  of  the  misfit  for  an  event.  The  low  in  the  north-west 

corner  is  the  best  source  location  for  an  event  in  bandwidth  0.8 -3. 5  Hz. 

A  similar  example  is  seen  in  figure  DlOb  for  an  event  one  day  earlier. 

Events  in  the  next  bandwidth,  3. 5 -6. 6  Hz,  have  more  poorly  defined  minima 

located  in  between  the  receivers  for  all  events.  High  misfits  occur  near 

2  2 

the  receivers  since  as  range  -  (x  +y  )  decreases,  P  increases  rapidly. 
Sources  of  error  in  calculating  these  locations  include  errors  in 
instrument  calibration,  non-spherically  spreading  waves,  or  the  existence 
of  multiple  sources. 

The  only  bandwidth  with  consistent  source  locations  from  event  to 
event  was  0.8 -3. 5  Hz  (figure  D2).  This  bandwidth  also  exhibits  the 
highest  overall  power  level  and  hence  may  be  less  affected  by  superposed 
isotropic  noise.  The  other  bandwidths  with  relatively  high  power  levels, 
show  great  scatter  in  source  location  over  time  and  hence  it  is 
impossible  to  attribute  them  to  a  stationary  geophysical  process.  The 
bandwidth  0.8-3.5HZ  sources  cluster  near  the  tip  of  a  small  ridge  on  the 
western  part  of  the  Mid-Atlantic  Ridge  (figure  D2) . 
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Harmonic  tremor  characterized  by  narrow  band  frequencies  in  the 
range  of  1-10  Hz  has  been  associated  with  eruptive  activity  and  magma 
movement  in  volcanoes  (Aki  et  al.,  1981;  Koyanagi  et  al . ,  1987). 
Persistent,  long  tremor  events  which  last  for  hours  and  recur  for  days 
probably  result  from  magmatic  activity  beneath  the  volcano.  In  Hawaii, 
this  shallow  (<5  km)  activity  is  confined  to  areas  as  small  as  5  km  in 
diameter  ce.  rered  on  the  summits  of  the  volcanoes.  We  postulate  that  the 
high  0.8-3. 5  Hz  signal  seen  in  the  TAG  area,  with  a  peak  at  2.5  Hz,  is 
produced  by  volcanic  processes,  such  as  magma  flow,  on  the  western  wall 
of  the  rift  valley. 


CONCLUSIONS 

The  ambient  sound  field  in  the  bandwidth  1-30  Hz  recorded  within  13 
km  of  the  TAG  hydrothermal  area  shows  considerable  spatial  and  temporal 
variation  over  the  course  of  fifteen  days.  The  hydrothermal  vent  did 
not  produce  enough  sound  to  be  detected  by  the  ocean  bottom  hydrophone 
array  located  from  0.75  to  14  km  away. 

Using  a  model  to  find  a  point  source  given  actual  power  levels  at 
each  of  six  OBH  instruments,  under  the  assumption  of  spherical  spreading, 
the  only  stationary  source  uncovered  was  generating  sound  in  a  0.8-3. 5  Hz 
bandwidth  and  located  outside  the  median  valley.  It  appears  to  be 
harmonic  tremor  associated  with  the  tip  of  a  ridge  on  the  western  side 
of  the  spreading  axis  and  may  be  volcanic  in  origin. 
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Microseisms 
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Spectra  of  ambient  noise  in  the  band  1  to  1000  Hz, based  upon 
measurements  in  the  literature  and  present  models  for  noise  spectra  in 
various  frequencies  ranges.  (From  Urick,  1986). 
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Figure  D2 

Topographic  nap  of  TAG  study  area  In  the  Mid-Atlantic  Ridge  (Rona,  et 
al.,  1986),  showing  hydrophone  locations  (large  open  circles), 
earthquake  locations  (snail  open  circles)  from  Kong  et  al.,  1986,  high 
tenperature  vent  location  (double  cross),  and  0.8 -3. 2  ..z  source 
locations  for  events  used  In  this  study  (solid  circles).  Contours  are 
In  100  •  Intervals,  nedlan  valley  runs  fron  north-east  to  south-west. 
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Figure  D4 

a)  Power  epecCra  for  Che  firet  record  (solid  line),  one  hour  later 
(dashed  line) ,  eleven  days  later  (dotted  line)  and  twelve  days 
later  (dot-dash  line),  from  the  lnstrunent  closest  to  the  vent, 
750n,  showing  tenporal  variation.  The  22  Hz  peak  nay  be  due  to 
the  R.V.  Knorr.  The  anomalously  high  levels  on  day  twelve  are 
unexplained  (dot-dash  line). 

b)  Power  spectra  for  the  first  record  (solid  line),  one  hour  later 
(dashed  line),  eleven  days  later  (dotted  line)  and  twelve  days 
later  (dot-dash  line),  fron  the  lnstrunent  12  ka  fron  the  vent 
showing  the  temporal  variation.  The  22  Hz  peak  nay  be  due  to  the 
R.V.  Knorr.  The  anomalously  high  levels  on  day  twelve  are 
unexplained  (dot-dash  line). 

138 


Pressure  Pa 
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Figure  DS 

Tlae  series,  showing  spatial  variability,  froa  Instrument  closest 
to  vent  (2),  and  instrument  6.5  ka  away  (1). 
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TAG  ls,2s,3s,5s,6s,8s 


Figure  D6 

«)  Power  epeccra  recorded  ae  each  of  six  Instruments  for  the  first 
record  ahov  two  orders  of  magnitude  range.  Note  22  Hz  peak  which 
may  be  caused  by  the  ship  or  a  whale . 
b)  Power  spectra  recorded  at  each  of  six  instruments  for  record  taken 
one  hour  after  figure  6a.  Note  22  Hz  peak  which  may  be  caused  by 
the  ship  or  a  whale. 
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TAG  lab,2ab,3ab.5ab,6ab,8ab 


Figure  D7 

a)  Power  spectra  recorded  at  each  of  six  Instruments  for  record  taken 
eleven  days  after  figure  6a.  Ship  la  absent.  Note  high  power  at 
low  frequencies. 

b)  Pewer  epectra  recorded  at  each  of  six  instruments  for  record  taken 
one  dey  after  figure  7a.  Shows  anomalously  high  broad  band  noise 
on  all  Instruments. 
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TAG  2s,3s,5s,ls,6s,8s 


Vent  — -Near  Far 


Figure  D8 

•)  Power  spectra  vs.  distance  from  vent  site  for  first  record  (same 
as  figure  6a).  Note  that  at  no  frequency  does  power  fall  off 
monotonlcally  with  distance.  The  22  Hz  ship  peak  can  be  seen  In 
all  Instruments. 

b)  Power  spectra  vs.  distance  from  vent  site  for  same  record  as 
figure  6b.  Note  that  at  no  frequency  does  power  fall  off 
monotonlcally  with  distance.  The  22  Hz  ship  peak  can  be  seen  In 
all  Instruments. 


TAG  2ab.3ab.5ab,  lab, dab.lab 


Vent  —  Near  Far 


Figure  D9 

•)  Power  spectra  vs.  distance  free  vent  site  for  sane  record  as 
figure  7a.  Ship  Is  absent.  There  Is  high  power  at  lew 
frequencies  but  power  does  not  fall  off  Mnotonlcally  with 
distance. 

b)  Power  spectra  vs.  distance  fron  vent  site  for  sane  record  as 
figure  7b,  showing  anomalously  high  broad  band  noise  on  all 
InatruMnCs  of  relatively  slallar  amplitude. 
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Figure  010 

a)  Contour  of  alsflt  between  actual  powar  va  calculated  power  for 
predictive  eourcee  within  the  atudy  area.  Contours  are  1.95x10'" 
Pa  /Hr  In  a  3.1  Hr  bandwidth  of  the  ras  aiaflt.  Lowest  point  Is 
the  beet  fitting  source  locstlon.  The  low  In  the  north-west 
corner  Is  the  location  of  the  source  that  produces  the  alnlaua 
alsflt  (0.0023  Pa  /Hr)  with  recorded  power  levels  for  the 
bandwidth  0.8-3. 5  Hr  for  the  event  shown  In  ftguras  7b  and  9b. 

The  ras  power  froa  the  six  lnstruaents  for  this  event  was  0.0062 
Pa'/Hr . 

b)  Contour  of  alsflt  betwsen  actual  power  vs  cslculatsd  power  for 
predictive  eources  within  the  study  area.  Contours  are  3.7x10' 
Pa^/Hz  In  a  3.1  Hr  bandwidth  of  the  ras  alsflt.  Lowest  point  Is 
the  best  fitting  source  location.  The  low  In  the  north-west 
corner  Is  the  location  of  the  source  thac  produces  the  alnlaua 
alsflt  with  recorded  power  levels  for  the  bandwidth  0.8-3. 5  Hr  for 
the  records  ehown  in  figures  7a  and  9a.  Contour  values  are 
different  froa  above  because  the  total  ras  power  Is  different 

for  the  two  events.  The, ras  power  froa  the  six  lnstruaents  for 
this  event  was  0.0125  PaVHr.  2.44 
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